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Preface 


Our  original  aspiration  with  regard  to  "Manual  Astronaut  Navi¬ 
gation"  was  to  prepare  a  set  of  charts,  tables,  and  simplified  proce¬ 
dures  the  astronaut  could  use  in  solving  the  orbit  determination 
problem.  But  we  soon  learned  that  a  great  deal  of  background  theory 
needed  development  or  refinement  before  such  a  flight  manual  approach 
could  be  justified.  Our  original  plan  of  attack  was  therefore 
changed  from  preparing  a  flight  manual  to  research  and  presentation  of 
a  variety  of  theoretical  approaches  to  the  manual  space  navigation 
problem.  Our  thesis  more  closely  resembles  a  textbook  approach  rather 
than  the  originally  intended  flight  manual. 

We  have  incorporated  and  presented  the  theoretical  development  of 
a  number  of  topics  applicable  to  the  problem  of  "Manual  Space  Navi¬ 
gation."  We  claim  no  truly  original  ideas  or  significant  "break¬ 
throughs."  In  fact,  most  of  the  concepts  we  have  built  upon  have 
been  known  for  centuries.  Only  recently,  though,  has  man  found  himself 
making  observations  from  the  vehicle  for  which  the  orbit  is  desired. 

In  the  past,  the  observation  point  has  been  rigidly  constrained  to 
the  surface  of  the  earth.  With  the  new  and  exciting  viewpoint,  many 
of  the  classical  orbit  determination  methods  assume  new  and  interesting 
possibilities.  We  hope  our  research  points  out  some  of  the  many 
advantages  the  astronaut  has  in  his  position  as  an  observer  aboard  the 
spacecraft. 
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Abstract 


Manual  Astronaut  Navigation,  onboard  orbit  determination  inde¬ 
pendent  of  primary  system  computers,  must  be  developed  to  insure  the 
safety  of  future  manned  space  flight.  This  report  presents  selected 
methods,  procedures,  and  equipment  that  form  a  basis  from  which  a 
flexible  operational  manual  navigation  system  may  eventually  be 
developed. 

The  Delta-H  method  of  graphically  determining  eccentricity  is 
simple  and  potentially  accurate.  Prepared  charts  are  entered  with 
two  changes  in  altitude,  or  "Delta-H* s,"  and  eccentricity  is  read 
directly.  An  estimate  of  the  size  of  the  orbit,  as  represented  by 
the  major  axis,  must  be  made  to  allow  selection  of  the  correct  chart. 

The  geometric  elements  of  the  orbit  can  be  obtained  by  use  of 
numerical  differentiation  of  a  series  of  range  measurements  to  obtain 
radial  velocity  and  radial  acceleration.  This  procedure  is  limited 
by  the  effects  of  truncation  and  measurement  errors.  Differential 
correction  of  the  orbital  elements  can  be  made  if  a  back-up  computer 
is  available  to  handle  the  computational  load. 

The  three-dimensional  position  fix  opens  the  door  to  many  orbit 
determination  possibilities.  Three  fixes  input  to  the  Gibbs  method 
is  the  most  promising  investigated.  .  This  procedure  may  be  performed 
by  hand  or  with  a  small  back-up  computer. 
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If  the  geometric  elements  are  known,  the  orientation  elements 
may  be  obtained  without  resorting  to  lengthy  iterative  techniques. 
Knowledge  of  angular  positions  at  two  times  coupled  with  the  true 
anomaly  at  one  of  those  times  is  all  that  is  required. 

Accurate  ranging  is  essential  in  most  orbit  determination  schemes. 
A  reference  trajectory  may  be  employed  with  linear  perturbation 
techniques  to  obtain  a  good  estimate  of  spacecraft  range.  Lineariza¬ 
tion  is  valid  for  relatively  large  deviations  from  the  reference 
trajectory.  For  midcourse  trajectories,  this  method  is  vastly 
superior  to  methods  incorporating  angular  diameter  measurements. 

Orbit  determination  procedures  exist  that  do  not  depend  upon 
range  measurements.  The  complete  orbit  may  be  obtained  through 
knowledge  of  three  angular  positions  and  three  true  anomaly  rates. 

To  date,  very  little  has  been  done  in  the  area  of  manual  naviga- 
tion  equipment  design.  "Aids"  to  manual  navigation  such  as  back-up 
inertial  platforms,  hand-held  mechanical  calculators,  and  battery 
operated,  hand-held  computers  are  necessary  to  optimize  the  performance 
of  future  manual  navigation  systems. 

The  manual  navigation  problem  can  be  solved  by  comtining  the 
proper  methods,  procedures,  and  equipment. 
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MANUAL  ASTRONAUT  NAVIGATION 


I.  Introduction 
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Background 

The  dramatic  success  of  the  Apollo  8  lunar  miosion  has  brought, 
this  nation  to  the  threshold  of  an  age  of  interplanetary  space  travel. 

i 

The  earth  orbits  of  projects  Mercury  and  Gemini,  once  considered  with 
awe,  are  now  classified  as  simple  "parking"  orbits  from  which  manned  .< 

*  -  t 

vehicles,  such  as  Apollo,  will  depart  upon  long  midcourse  trajecto¬ 


ries  to  the  moon  and  planets.  As  in  the  flight  of  Apollo  8,  the 


CP 


success  of  these  missions  will  depend  heavily  upon  accurate  and  ° 
reliable  navigation  and  guidance.  In  the  event  of  primary  system0 

•C 

failure,  an  adequate  back-up  system  must  be  available  to  allow^  safe 
return  to  earth.  Manual  navigation  and  guidance  procedures  should "be 

"O  v*- 
o 

provided  to  eliminate  all  dependence  on  primary  system  computers. $ 

Such  procedures  would  involve  use  of  equipment  such  as  the  sp&ce 
stadimeter,  a  space  sextant,  prepared  tables  and  charts,  and  hand-held 
computing  devices.  o  *> 
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The  Problem  ... 

This  report  is  intended  to  be  a  continuation  of  the  work  accora- 

O  ,  -v 

plished  by  Captains  Richard  R.  Schehr  and  Patrick  J.  Smith  in  their 
thesis,  "Manual  Astronaut  Navigation:  Apollo  Mission  Applications," 
(Ref  17) .  The  term  "navigation"  is  defined  as  the  determination  of 
any  six  independent  parameters  which  specify  uniquely  the  geometric 
qualities  of  an  orbit,  or  trajectory,  and  its  orientation  in  inertial 


-'■Vk 


1 


GA/AE/69-1 


space.  "Guidance"  is  defined  as  the  subsequent  determination  of  the 
velocity  correction  vector  necessary  to  insure  the  arrival  of  the 
•  spacecraft  at  a  predetermined  point.  "Manual"  navigation,  the  problem 
attacked  in  this  thesis,  is  defined  as  the  onboard  determination  of 
the  orbital  parameters  using  procedures  and  equipment  totally  inde¬ 
pendent  of  the  primary  system  computer. 

Scope 

y  — 

The  previous  thesis  dealt  with  manual  navigation  in  low  altitude, 
near-circular  earth  and  lunar  orbits.  This  thesis  will  explore  general 
manual  navigation  techniques  applicable  to  midcourse  trajectories  such 
as  the  translunar  phase  of  the  Apollo  mission.  No  attempt  is  made  to 
consolidate  the  procedures  described  in  an  operational  "flight  manual" 
presentation,  as  was  done  by  Captains  Schehr  and  Smith.  The  principal 
objectives  are  to  extend  and  improve  upon  their  methods,  and  to 
investigate  other  promising  manual  navigation  techniques.  Sources  of 
error  are  discussed,  and  some  error  analyses  are  presented.  Numerical 
examples  are  provided,  wherever  possible,  to  enhance  the  theoretical 
development. 

Assumptions 

The  primary  assumption  made  in  this  thesis  is  that  of  restricted 
two-body  motion.  This  implies  that: 

1.  The  mass  of  the  vehicle  is  negligible  when  compared 
to  the  mass  of  the  attracting  body. 

2.  The  gravity  field  of  the  attracting  body  is  ideal, 
inverse-square,  and  central. 
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The  only  other  assumption  made  is  that  all  measurements  indicated  can 
be  performed. 

The  first  assumption  allows  all  orbits  and  trajectories  to  be 
considered  as  conic  sections.  This  assumption  is  valid  for  near  planet 
orbits  and  the  major  part  of  midcourse  trajectories.  The  second 
assumption  is,  of  course,  a  practical  necessity  in  any  theoretical 
discussion  of  manual  navigation. 

The  reader  is  assumed  familiar  with  the  basics  of  astrodynamics 
and  celestial  mechanics. 

Organization 

Each  of  Chapters  II  through  VII  in  this  report  deals  with  an 
independent  method  of  obtaining  one  or  more  of  the  six  orbital  param¬ 
eters  necessary  for  solution  of  the  manual  navigation  problem.  These 
parameters,  or  elements,  are  divided  into  two  sets:  the  geometric 
parameters, 

e  -  eccentricity 
a  -  semimajor  axis 
t  -  time  from  perifocus 
and  the  orientation  parameters, 

i  -  inclination  of  the  orbital  plane 
ft  -  longitude  of  the  ascending  node 
w  -  argument  of  perifocus 

Chapter  VIII  deals  with  an  improved  method  of  optical  ranging,  and 
Chapter  IX  is  a  discussion  of  certain  "aids"  to  manual  navigation 
such  as  battery  operated,  hand-held  computers. 
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Discussion 

As  stated  previously,  this  is  not  a  "flight  manual."  No  "best" 
set  of  procedures  is  determined.  It  will  be  seen  that  this  determina¬ 
tion  depends  upon  the  answers  to  several  questions.  Most  important 
are: 

1.  What  back-up  computational  tools  are  available? 

2.  What  types  of  measurements  are  possible? 

3<  How  much  time  is  allowable  for  the  completion  of  the 
procedures? 

4.  What  is  the  required  accuracy? 

This  thesis  forms  a  basis  from  which  an  operational  flight  manual  may 
eventually  be  developed  once  the  above  questions  have  been  answered. 
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II.  The  Delta-H  Method 

Introduction 

The  cornerstone  of  the  geometric  parameters  determination  for 
low  eccentricity  elliptical  orbits  in  the  thesis,  "Manual  Astronaut 
Navigation:  Apollo  Mission  Applications,"  by  Captains  Schehr  and 
Smith,  is  what  they  labeled  the  "Delta-H"  method  of  obtaining  orbital 
eccentricity.  This  is  the  first  of  three  independent  geometric 
elements  obtained  in  their  approach  to  the  manual  navigation  problem. 

The  other  two  parameters  are  the  time  from  peripoint  and  the  period. 

In  their  scheme,  once  eccentricity  is  determined  it  is  used  to  obtain 
graphically  the  ratio  of  time  from  peripoint  to  period.  The  period 
is  also  determined  graphically  by  a  comparison  with  a  zero  eccentricity, 
or  circular,  orbit. 

As  will  be  seen,  the  Delta-H  method  of  obtaining  orbital  eccen¬ 
tricity,  original  to  the  aforementioned  thesis,  is  relatively  simple 
and  can  be  sufficiently  accurate.  For  this  retason,  and  because  the 
Delta-H  concept  was  not  fully  developed  due  to  lack  of  time,  this 
chapter  is  devoted  primarily  to  an  investigation  of  this  method  of 
obtaining  orbital  eccentricity.  Particularly,  an  extension  of  the 
previous  work  to  highly  eccentric  geocentric  ellipses  such  as  the 
translunar  phase  of  the  Apollo  mission  is  accomplished  as  a  follow- 
on  to  the  original  low  eccentricity  development. 
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Background 

The  Delta-H  method  of  eccentricity  determination  involves  the 
measurement  of  three  consecutive  altitudes  (or  ranges)  separated  by 
equal  time  intervals.  The  second  altitude  (range)  is  subtracted  from 
the  first,  and  the  third  from  the  second,  to  give  two  changes  in 
altitude,  or  "Delta-H* s."  By  definition 

AH1  =  h1  -  h2  (2-1) 

and 

AH2  =  h2  -  h3  (2-2) 

In  their  investigation.  Captains  Schehr  and  Smith  discovered  that 
if  the  changes  in  altitude  were  taken  in  pairs  successively  along  the 
ellipse  and  plotted  as  ordinate  and  abscissa  in  an  orthogonal  axis 
system,  curves  of  constant  eccentricity  resulted  (Ref  17:18).  Figure 
1  is  an  example  of  one  of  their  Delta-H  plots.  To  6btain  eccentricity, 
the  plot  is  entered  with  AH1  and  AH2  and  eccentricity  is  read  directly. 
The  curves  were  generated  for  time  intervals  of  5  and  15  minutes  along 
geocentric  and  selenocentric  orbits.  Only  very  low  eccentricities 
were  considered;  from  .002  to  .016  in  geocentric  orbit,  and  from  .01 
to  .08  in  selenocentric  orbit.  Most  importantly,  for  the  range  of 
major  axes  considered,  i.e.,  from  planet  radius  plus  80  nautical  miles 
to  planet  radius  plus  170  nautical  miles,  the  curves  appeared  to  be 
invariant.  This  has  since  been  confirmed  by  the  authors  of  this  thesis 
as  well  as  by  the  Universal  Technology  Corporation  working  under  Air 
Force  contract.  The  apparent  Insensitivity  to  the  size  of  the  orbit, 
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as  represented  by  the  major 
axis,  means  that  relatively 
few  plots  are  necessary  to 
cover  a  sizable  range  of 
orbits.  This  advantage 
makes  the  Delta-H  method 
more  practical  than  the 
previously  developed  graph¬ 
ical  approaches  (Ref  17:18). 

Operationally,  the  accuracy 
of  the  Delta-H  method  de¬ 
pends  upon  the  accuracies  Fig.  1 

of  timing  and  altitude  (or  Typical  Low  Eccentricity,  Low 

Altitude  Delta-H  Plot 

range) measurement.  These,  (From  Ref  17:126) 

in  turn,  depend  upon  the 

hardware  available,  the  skill  of  the  user,  and  the  inherent  accuracy 
of  the  measurement  procedure. 

Within  the  scope  of  the  problem  defined,  the  Delta-H  method 
appeared  to  Captains  Schehr  and  Smith  to  be  a  simple,  direct,  and 
potentially  accurate  approach  to  obtaining  orbital  eccentricity  man¬ 
ually.  Motivated  by  their  work,  the  rest  of  this  chapter  deals  with 
the  generation  of  Delta-H  plots,  an  investigation  of  their  properties, 
and  a  discussion  of  the  possibilities  of  using  them  in  a  high  eccen¬ 
tricity  manual  navigation  scheme. 


-40  -30  -20-10 

0  10  20  1C  40 

a  nz 

(Ml 

GA/AE/69-1 


Data  Generation 

A  double-precision  computer  program  was  written  to  provide  accurate 
reference  data  for  use  in  investigating  the  Delta-H  method  and  other 
manual  navigation  techniques  to  be  described  in  later  chapters.  The 
computer  graphs  the  Delta-H  data  via  a  tape-driven  digital  plotter. 
Details  are  provided  here  to  minimize  the  amount  of  backtracking 
necessary  in  future  follow-on  work. 

Theory .  As  stated  in  Chapter  I,  a  two-body  approach  is  taken  in 
this  thesis.  The  equations  used  in  the  development  of  the  computer 


program  are 


(2-3) 


(2-4) 


where 


V2  = 

(2-5) 

h  = 

r  -  R 

(2-6) 

M  « 

nt  =  E  -  e  sin  E 

(2-7) 

r  ~ 

a(l  -  e  cos  E) 

(2-8) 

n  =  mean  motion 
p  =  gravitational  constant 
a  =  seraimajor  axis 
t  =  period 
V  *  total  velocity 
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e  =  eccentricity 
r  =  radial  distance 
h  =  altitude 
M  =  mean  anomaly 
t  =  time  from  peripoint 

E  =  eccentric  anomaly  measured  from  peripoint 
R  =  planet  radius 

It  is  desired  to  have  the  computer  solve  for  altitude  as  a  function 
of  semimajor  axis,  eccentricity,  and  time  from  peripoint.  This  involves 
the  solution  of  Eq  (2-7)  for  eccentric  anomaly  as  a  function  of  the 
same  three  variables,  and  then  Eq  (2-8)  for  radial  distance.  Altitude 
is  then  obtained  from  Eq  (2-6) .  Difficulties  arise  because  Eq  (2-7) , 
Kepler's  equation,  is  transcendental  in  eccentric  anomaly.  This 
problem  is  not  serious  for  low  eccentricities  since  many  excellent 
approximate  solutions  have  been  formulated  (Ref  A: 41).  Most  of  these 
solutions  are  in  the  form  of  series  expansions.  Captains  Schehr  and 
Smith  used  the  first  three  terms  of  a  series  expansion  in  powers  of 
eccentricity  about  the  mean  anomaly  (Ref  17:17).  This  series  diverges 
for  values  of  eccentricity  above  .662743  (Ref  11:11-5).  Consequently, 
to  investigate  orbits  with  eccentricities  above  .9,  a  different 
approach  had  to  be  formulated.  Further  research  disclosed  that  all 
series  expansions  for  eccentric  anomaly  diverge  in  this  region,  so  it 
was  decided  to  employ  an  iterative  procedure  based  upon  the  Newton 
method  of  approximations.  Rearranging  Eq  (2-7) 

sin(E)  =  -  (E  -  M)  (2-7a) 

e 
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or 

f (E)  =  sin  E  -  -  (E  -  M)  (2-7b) 

e 

When  Eq  (2-7)  is  satisfied,  f(E)  is  identically  zero.  An  estimate"  of 
eccentric  anomaly  is  made  to  start  the  iteration.  Then,  the  Newton 
method  is  employed  to  obtain  a  better  estimate,  or  second  approxima¬ 
tion. 


E„  =  E,  - 


f(E) 


d  f  (E) 
dE 


or 


-  E,  - 


sin  EX  -  e  (Et  -  M) 
cos  Ei  -  i 


(2-9) 


(2-9a) 


This  procedure  is  repeated  until  f(E)  is  zero.  The  Newton  method  is 
generalized  in  the  text  "Astrodynamics — Advanced  Topics  and  Applica¬ 
tions,"  by  Baker  (Ref  3:26).  Since  the  altitude  computations  start  at 
the  peripoint,  the  initial  value  of  eccentric  anomaly  is  zero. 

Altitudes  are  to  be  computed  for  every  2.5  minutes  thereafter.  This 
short  time  interval  makes  the  final  value  of  eccentric  anomaly  at  one 
time  an  excellent  initial  value  for  the  next.  Also,  once  the  desired 
altitudes  are  obtained,  they  can  be  combined  in  pairs  to  give 
Delta-H's  for  timing  intervals  of  any  integral  multiple  of  2.5  minutes. 
Thus,  the  effect  of  the  timing  interval  can  be  investigated  and  the 
number  of  iterations  for  convergence  is  kept  small. 
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The  Computer  Program.  The  computer  program  mentioned  at  the  begin 
ning  of  this  section  is  reproduced  in  Appendix  A.  The  statements 
required  for  automatic  plotting  are  not  included  since  they  complicate 
the  program  considerably  and  are  not  essential  for  an  understanding  of 
the  logic  involved.  Input  to  the  program  are  semimajor  axis,  eccen¬ 
tricity,  and  time  from  peripoint.  The  output  includes  period,  altitude 
AH1,  AH2,  eccentric  anomaly,  and  the  ratio  of  velocity  to  the  velocity 
at  peripoint. 

High  Eccentricity  Delta-H  Plots 

The  translunar  phase  of  the  Apollo  mission  is  a  geocentric  ellipse 
with  a  typical  eccentricity  of  .9668  and  an  apogee  distance  of  approxi¬ 
mately  60  earth  radii  (Ref  1) .  These  figures  are  used  as  a  point  of 
departure  in  this  section. 

The  computer  was  programmed  to  generate  geocentric  Delta-H  plots 
for  seven  values  of  eccentricity  for  each  of  seven  values  of  semimajor 
axis.  Thus,  the  effect  of  the  semimajor  axis  can  be  isolated  and 
studied.  The  eccentricity  was  initialized  at  .9068  and  incremented  by 
.0100  until  the  .9668  value  was  attained.  The  semimajor  axis  was 
initialized  at  105,000  nautical  miles  and  incremented  by  1,000  nautical 
miles  until  a  final  value  of  111,000  nautical  miles  was  attained.  This 
set  of  values  for  eccentricity  and  semimajor  axis  gives  a  minimum 
apogee  distance  of  approximately  58  earth  radii  and  a  maximum  of 
approximately  63  earth  radii,  since  at  apogee 

r  =  a  (1  +  e)  (2-10) 


from  Eq  (2-8) .  The  timing  interval  was  set  at  15  minutes  with  the 
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first  altitude  determined  at  perigee  and  the  last  3.5  hours  from  perigee. 
From  the  point  of  view  of  the  Apollo  mission,  this  represents  the 
portion  of  the  translunar  coast  phase  prior  to  the  first  scheduled 
midcourse  correction  (Ref  1) .  To  insure  a  sufficient  number  of  data 
points  for  accurate  curve  plotting,  the  Delta-H  pairs  were  computed 
using  2.5  minute  displacements  of  the  three  altitudes,  i.e.,  the  first 
pair  consists  of  altitudes  computed  at  times  0,  15,  and  30  minutes;  the 
second  pair  consists  of  altitudes  computed  at  times  2.5,  17.5,  and 
32.5  minutes,  and  so  on.  Figure  2  is  one  of  the  resulting  Delta-H 
plots.  (To  keep  the  data  points  in  the  first  quadrant,  -AH1  is  plotted 
versus  -AH2.)  The  apparent  overlapping  observed  is  highly  undesirable 
since  it  produces  ambiguity  for  certain  combinations  of  AH1  and  AH2. 
Clearly,  the  curves  of  constant  eccentricity  must  be  distinct  every¬ 
where  for  these  plots  to  have  any  real  value  in  a  manual  navigation 
scheme.  (NOTE:  THE  COMPUTER  SCALES  AUTOMATICALLY  AND  USES  THE 
NOTATION  X  101  TO  MEAN  MULTIPLY  THE  NUMBER  ON  THE  AXIS  BY  TEN.  THIS 
MUST  BE  KEPT  IN  MIND  SINCE  IT  IS  THE  OPPOSITE  OF  WHAT  IS  USUALLY  SEEN 
IN  ENGINEERING  PRACTICE.  THE  LIBRARY  SUBROUTINE  WAS  NOT  ALTERED  DUE 
TO  THE  PROHIBITIVE  COST  OF  RERUNNING  THE  PLOTS  REFERRED  TO  IN  THIS 
CHAPTER. ) 

The  Effect  of  the  Timing  Interval.  Intuitively,  it  was  felt  that 
perhaps  the  timing  interval  was  too  short  to  allow  adequate  "sampling" 
of  the  trajectories  in  the  region  of  ambiguity.  Working  on  this  possi¬ 
bility,  the  program  was  altered  by  increasing  the  timing  interval  from 
15  minutes  to  30  minutes.  Figure  3  is  the  resulting  plot  for  the  same 
semimajor  axis  as  in  Fig.  2.  Some  improvement  is  noted  in  that  the 
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High  Eccentricity  Delta-H  Plot  for 
a  =  111,000  NM  With  At  =  30  Minutes 


Fig.  4 

ty  Delta-H  Plot  for 
With  At  =  60  Minutes 


High  Eccentric! 
a  =  111,000  NM  1 
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curves  are  generally  better  separated  than  those  in  the  prt-  ious  plot. 
The  timing  interval  was  changed  to  60  minutes  with  the  result  shown 
in  Fig.  A.  Complete  separation  is  attained,  which  indicates  that 
regions  of  ambiguity  can  be  eliminated  by  simply  increasing  the 
timing  interval.  The  effect  of  the  timing  interval  on  the  high 
eccentricity  Delta-II  plots,  then,  is  to  increase  or  decrease  the  rela¬ 
tive  curve  separation  as  the  timing  interval  is  made  correspondingly 
larger  or  smaller. 

The  Effect  of  the  Major  Axis.  As  mentioned  previously.  Captains 
Schehr  and  Smith  found  that  their  constant  eccentricity  curves  appeared 
to  be  invariant.  However,  the  high  eccentricity  curves  considered 
here  shift  slowly  as  the  major  axis  is  changed.  This  is  made  apparent 
by  comparison  of  the  plots  in  Figs.  5,  5a,  and  5b  with  Fig.  A.  As 
the  major  axis  is  decreased,  the  curves  shift  slowly  upward  and  to 
the  right.  It  is  probable  that  Captains  Schehr  and  Smith  did  not 
encounter  this  dependence  upon  the  major  axis  because  of  the  limited 
range  of  values  investigated.  In  che  present  case,  there  must  be  a 
change  of  2,000  nautical  miles  in  the  major  axis  to  materially  affect 
the  Delta-H  plots. 

To  help  in  establishing  the  general  effect  of  the  major  axis  on 
these  high  eccentricity  plots,  several  computer  runs  were  made  with 
a  variety  of  values  of  major  axis.  The  conclusion  drawn  upon 
examination  of  the  results  is  that  for  high  eccentricity  ellipses 
(e  above  .9),  a  variation  of  more  than  1  per  cent  in  the  major  axis 
above  or  below  a  particular  value  is  necessary  to  displace  the 
constant  eccentricity  curves  a  significant  amount.  The  Delta-H 
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High  Eccentricity  Delta-H  Plots  for  a  =  108,000  NM 
and  a  =  107,000  NM  With  At  =  60  Minutes 


Fig.  5b 


High  Eccentricity  Delta~H  P.__ts  for  a  =  106,000  NM 
and  a  =  105,000  NM  With  At  =  60  Minutes 
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plots  shown  in  Figs  2  through  5b  are  considered  typical., 

The  Analytic  Approach 

Attempts  have  been  made  to  analyze  and  predict  the  behavior  of 
the  Delta-H  curves  mathematically.  A  general  expression  (valid  for 
any  eccentricity  less  than  1)  relating  the  changes  in  altitude 
directly  to  the  major  axis,  timing  interval,  and  the  eccentricity  such 
as 


AH1  =  f (AH2,a,At,e) 


(2-11) 


or 


AH2  =  f (AHl,a,At,e)  (2-lla) 

has  been  sought,  but  all  attempts  at  formulating  such  an  expression 
have  failed  due  to  the  transcendental  nature  of  Kepler's  equation 
which  links  time  to  altitude  via  Eqs  (2-6)  through  (2-8) .  With  these 
equations  plus  Eqs  (2-1)  and  (2-2),  it  can  be  shown  that,  for  alti¬ 
tudes  1,  2,  and  3 

AH1  =  ae(cos  E?/  -  cos  Ej)  -  AH2  (2-12) 

If  s,  e,  At,  and  AH2  are  specified,  E1  and  E3  can  be  obtained  from 
the  simultaneous  solution  of 


and 


£H2  =  ae(cos  E3  -  cos  E2) 


nAt  =  (E3  -  E2)  -  e(sin  E3  -  sin  E2) 


*  *>\ 


(2-14) 
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which  yields  E2  and  E3,  and  then  the  solution  of  either 

nAt  =  (E2  -  Ej)  -  e(sin  E2  -  sin  E^)  (2-15] 

or 

2nAt  =  (E3  -  Ej)  -  e(sin  Eg  -  sin  Ej)  (2-16] 

which  yields  Ej .  Therefore,  Eq  (2-12)  can  be  solved  for  AH1  given 
AH2,  a,  At,  and  e.  First,  however,  E1 ,  E2,  and  E3  must  be  obtained  by 
iteration.  The  necessity  of  taking  this  intermediate  step  seems  to 
preclude  the  possibility  of  obtaining  an  explicit  relation  such  as 
Eq  (2-11)  or  (2-] la)  for  use  in  a  mathematical  analysis.  An  empirical 
approach  has,  therefore,  been  taken  in  this  chapter.  The  computer 
program  provides  an  implicit  solution  of  Eq  (2-12)  so  the  individual 
effects  of  varying  a,  e,  and  At  can  be  observed. 

General  Characteristics  of  the  Delta-H  Plot 

To  determine  some  general  characteristics  of  the  Delta-H  plot, 

15  computer  runs  were  made  for  intermediate  (.50  to  .75)  and  low 
(.01  to  .10)  values  of  eccentricity.  Both  regions  were  examined  for 
sensitivity  to  changes  in  the  timing  interval  and  the  major  axis. 

Figure  6  is  a  plot  for  eccentricities  ranging  from  .50  to  .56 
in  increments  of  .01.  The  timing  interval  is  15  minutes  and  the 
semimajor  axis  is  106,000  nautical  miles.  Figures  7,  7a,  and  7b  are 
a  set  of  plots  for  the  same  values  of  eccentricity,  but  with  a  60 
minute  timing  interval.  The  major  axis  differs  from  plot  to  plot  by 
2,000  nautical  miles.  These  plots  display  typical  intermediate 
eccentricity  characteristics.  Their  behavior  is  very  similar  to 
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Intermediate  Eccentricity  Delta-H  Plot  for 
a  =  106,000  NM  With  At  =  15  Minutes 


GA/AE/69-1 


Fig.  7a 

Intermediate  Eccentricity  Delta-H  Plots  for  a  =  105,000  NM 
and  a  =  104,000  NM  With  At  =  60  Minutes 
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that  noted  in  the  high  eccentricity  region. 

Figures  8,  8a,  and  8b  are  a  typical  set  of  low  eccentricity  plots. 
Separation  was  again  attained  with  a  60  minute  timing  interval.  A 
notable  difference  in  the  behavior  of  these  plots  lies  in  the  lessened 
effect  of  changes  in  the  major  axis.  A  change  of  more  than  3  per  cent 
appears  necessary  to  displace  the  curves  a  significant  amount. 

Conclusions.  The  general  conclusions  drawn  are: 

1.  Curve  separation  for  all  eccentricities  can  be 
controlled  by  changing  the  timing  interval. 

2.  Sensitivity  to  changes  in  the  major  axis  is  variable. 

It  appears  that  at  low  eccentricities,  a  change  of 
more  than  3  per  cent  is  necessary  for  significant 
displacement,  and  in  the  high  and  intermediate 
regions,  a  change  of  more  than  1  per  cent  is  necessary. 

Conclusion  2  tends  to  confirm  the  supposition  that  Captains  Schehr  and 
Smith  failed  to  encounter  the  dependence  upon  the  map  or  axis  due  to 
the  limited  range  of  values  investigated.  They  considered  low  eccen¬ 
tricities  and  a  maximum  change  of  approximately  2  1/2  per  cent  in  the 
major  axis. 

While  this  has  not  been  a  comprehensive  analysis  of  the  Delta-H 
plot  characteristics  (since  such  an  analysis  is  beyond  the  scope  of 
this  thesis) ,  two  important  trends  have  been  discovered  in  the  effects 
of  the  timing  interval  and  the  major  axis.  These  trends  will  be  us  2d 
as  a  basis  for  fitting  the  Delta-H  method  of  eccentricity  determina¬ 
tion  into  a  high  eccentricity  manual  navigation  scheme. 
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Fig.  8 

Low  Eccentricity  Delta-II  Plots  for  a  =  106,000  NM 
and  a  =  105,000  NM  With  At  =  60  Minutes 
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Low  Eccentricity  Delta-H  Plots  for  a  -  104,000  NM 
and  a  =  103,000  NM  With  At  -  60  Minutes 
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Discussion 

It  has  been  seen  that  the  Delta-H  plots  are  sensitive  to  changes 
in  the  timing  interval  and  the  major  axis.  These  sensitivities  must 
be  taken  into  account  and  used  to  best  advantage  if  the  method  is  to 
be  extended  beyond  near  circular  parking  orbits  to  midcourse 
trajectories  such  as  the  translunar  phase  of  the  Apollo  mission. 

It  would  be  unrealistic  to  i.  'it  the  Apollo  midcourse  problem 
by  assuming  a  small  (±  1  per  cent)  variation  in  the  major  axis  so 
that  only  one  Delta-H  plot  is  necessary  to  cover  each  set  of  possible 
eccentricities.  Au  inflight  emergency  involving  only  moderate 
deviations  from  the  planned  injection  velocity  could  easily  place  the 
actual  major  axis  outside  this  range;  see  Eq  (2-5).  Therefore,  a 
practical  manual  scheme  incorporating  the  Delta-H  method  would  have 
to  include  several  plots  for  each  set  of  eccentricities  to  assure 
adequate  coverage  for  such  a  contingency.  A  set  such  as  the  plots  in 
Figs.  4  through  5b  would  suffice  for  a  12,000  nautical  mile  range  of 
major  axes,  and  additional  plots  could  be  carried  to  increase  this 
coverage.  The  astronaut  would  have  to  know  the  major  axis  to  within 
2,000  nautical  miles  to  choose  the  correct  plot  for  accurate  determi¬ 
nation  of  the  orbital  eccentricity.  Larger  errors  in  the  major  axis 
are  tolerable  if  some  uncertainty  in  the  eccentricity  is  acceptable. 
Careful  inspection  of  Figs.  4  through  5b  shows  that  if  the  major  axis 
is  known  to  within  plus  or  minus  6,000  nautical  miles  (about  ±  3 
per  cent) ,  the  uncertainty  in  the  eccentricity  is  no  greater  than 
approximately  .003.  This  can  be  demonstrated  by  assuming  a  combina¬ 
tion  of  AH1  and  AH2  and  then  entering  each  of  the  seven  plots  for  an 
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estimate  of  the  eccentricity.  (If  the  central,  or  fourth,  plot  is 
considered  to  be  the  correct  plot,  then  this  set  represents  a  ±  3 
per  cent  variation  from  the  actual  major  axis.)  The  difference 
between  the  highest  or  lowest  estimate  and  the  true  eccentricity  is 
in  no  case  greater  than  .003.  Therefore,  this  represents  an  upper 
bound  on  the  uncertainty  in  the  eccentricity  given  that  the  error 
in  the  major  axis  is  within  6,000  nautical  miles. 

The  errors  assumed  in  the  previous  paragraph  are  considered  to 
be  the  result  of  inaccuracies  in  the  method  used  to  determine  the 
major  axis.  This  does  not  include  the  effect  of  measurement  errors. 
Measurement  errors,  of  course,  affect  both  the  determination  of  the 
major  axis  and  the  two  changes  in  altitude.  It  is  essential  that 
the  range  (or  altitude)  measurements  be  as  accurate  as  possible. 
However,  it  is  evident  from  the  relatively  large  changes  in  a] tit  e 
noted  in  Figs.  4  through  5b  that  small  errors  can  be  absorbeu  ch 
little  affect  on  the  eccentricity  determination  as  long  as  the  major 
axis  is  known  to  within  the  limits  prescribed  above.  Also,  an 
advantage  of  the  Delta-H  method  lies  in  the  elimination  of  measure¬ 
ment  errors  due  to  human  bias.  These  errors  tend  to  be  self¬ 
canceling  (Ref  17:27). 

The  major  axis  can  be  obtained  from  knowledge  of.  the  orbital 
period;  see  Eq  (2—14) .  The  method  used  by  Captains  Schehr  and 
Smith  is  not  applicable  in  this  case  since  it  assumes  eccentricities 
close  to  zero.  More  flexible  methods  of  determining  the  major  axis 
are  discussed  later  in  this  thesis.  Once  the  major  axis  and  the 
eccentricity  are  obtained,  they  can  be  used  to  find  the  final 
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geometric  parameter,  the  time  from  peripoint  to  one  of  the  measurements 
from 


where 


_  E  -  e  sin  E 
n 


(2-17) 


E  =  cos-1  (2-18) 

and  r  is  the  range  at  the  time  of  interest.  If  the  cosine  of  E  is 
positive  and  the  range  is  increasing,  then  E  is  less  than  90  degrees. 

If  the  range  is  decreasing,  the  angle  is  greater  than  270  degrees,  and 
so  on. 

The  curve  separating  effect  cf  increasing  the  timing  interval  can 
be  used  to  generate  Delta-H  plots  for  eccentricity  increments  of  .001 
or  smaller,  thereby  improving  the  accuracy  with  which  the  graphs  may 
be  read.  Figure  9  is  a  plot  for  a  90  minute  timing  interval,  and 
Fig.  10  is  a  portion  of  the  same  plot  with  .001  increments. 

Sample  Problem.  The  following  is  a  sample  calculation  of  the 
eccentricity  and  time  from  peripoint  on  a  known  ellipse.  The  calcu¬ 
lations  were  performed  with  a  slide  rule,  standard  trigonometric 
tables,  and  a  scratch  pad.  The  errors  introduced  are  assumed  typical 
for  the  measurements  and  methods  involved. 

On  the  known  ellipse  (exact  values) 


a  =  105,000  NM  (NM  for  all  distances) 


e  «  .9588 

hj  =  7245.584  at  50  minutes  from  perigee 
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h2  =  15540.412  at  110  minutes  from  perigee 
h3  =  22629.467  at  170  minutes  from  perigee 
The  assumed  measurements,  with  errors,  are 

h  =  7300 

h2  =  15650 

h3  =  22750 
At  =  60  minutes 

The  assumed  calculated  major  axis  is  108,000  nautical  miles . 
With  this  input  data 


AH1 

=  hl 

-  h2 

=  -8350 

AH2 

h2 

"  h3 

=  -7100 

Entering  the  AH  plot  for  a  =  108,000  NM  yields 


e  =  .9598 


Then 


E 


1 


-1 

=  coo 


a  -  rl 
ae 


a  -  t1  =  97260 
ae  =  103658 
Ej  =  cos-1  (.9382) 
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Ej  =  20°  15'  =  .3534  radians  (range  increasing) 

sin  Ej  =  .34612 

n  =  .0254 


.3534  -  .3322 
.0254 


minutes 


t1  =  (.8346)  (60)  =  50.07  minutes  from  perigee 


Summary .  The  positive  and  negative  aspects  of  using  the  Delta-H 
method  in  a  high  eccentricity  scheme  may  be  summarized  as  follows: 

Positive  Aspects 

1.  The  computations  required  are  simple. 

2.  An  uncertainty  of  no  more  than  about  .003  in  the 
eccentricity  is  possible  with  a  ±  3  per  cent  error  in 
the  estimate  of  the  major  axis. 

3.  Small  measurement  errors  can  be  absorbed  with  little 
adverse  affect  as  long  as  the  estimate  of  the  major 
axis  is  within  the  limits  prescribed. 

4.  Human  bias  as  a  source  of  error  is  eliminated. 

5.  If  an  increased  timing  interval  is  acceptable,  the 
Delta-H  curves  can  be  separated  enough  to  allow 
eccentricity  increments  of  .001  or  smaller. 

Negative  Aspects 

1.  A  sufficient  number  of  plots  must  be  carried  to  assure 
adequate  coverage  of  possible  orbits. 

2.  The  major  axis  must  be  known  to  within  ±  1  per  cent 
of  the  actual  value  for  accurate  eccentricity 


determination. 
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3.  Measurement  inaccuracies  can  give  totally  erroneous 
estimates  of  the  eccentricity  through  errors  in 

a,  AH1  and  AH2. 

4.  The  timing  interval  must  be  large  enough  to  eliminate 
ambiguity.  Therefore,  it  may  take  two  or  more  hours 
to  obtain  AH1  and  AH2. 

All  of  these  factors  must  be  weighed  before  inclusion  of  the 
Delta-H  method  in  a  high  eccentricity  manual  navigation  scheme. 
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III.  Geometric  Parameters  from  Numerical  Differentiation 

Introduction 

The  following  geometric  properties  of  an  orbit  can  be  obtained 
from  a  set  of  timed  range  (or  altitude)  measurements: 

1.  the  semimajor  axis  -  a 

2.  the  eccentricity  -  e 

3.  the  eccentric  anomaly  E  associated  with  a  position 
and  time. 

The  range  r  is  the  vehicle's  distance  from  the  center  of  the  earth. 

The  computation  method  considered  in  this  chapter  involves  a  numerical 
differentiation  to  obtain  the  radial  components  of  velocity  and 
acceleration  r  and  r.  Once  these  quantities  have  been  obtained,  the 
geometric  elements  may  be  computed  from  the  two-body  orbits..  equations. 
A  block  diagram  outlining  the  general  procedure  is  shown  in  Fig.  11. 

Knowledge  of  the  eccentric  anomaly  E  at  a  specific  time  t  is 
sufficient  to  relate  position  to  time  through  Kepler's  equation 

E  -  e  sin  E  =  M  =  (t  -  t0)  ^  iL.  (3-1) 

where  M  is  the  mean  anomaly,  tQ  the  perifocus  time,  and  y  the  earth's 
gravitational  constant.  In  this  chapter,  the  eccentric  anomaly  at 
its  associated  time  will  be  considered  as  one  of  the  orbital  elements, 
thus  simplifying  portions  of  the  presentation. 
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Fig.  11 

Determination  of  Geometric 
Orbital  Parameters 

Numerical  Differentiation  for  _r  and  i? 

The  radial  velocity  and  radial  acceleration  can  be  obtained  by 
numerically  differentiating  a  set  of  range  measurements.  Development 
and  tabulation  of  numerical  differentiation  formulas  can  be  found  in 
any  text  on  numerical  analysis.  In  general,  numerical  differentiation 
is  achieved  by  using  the  first  few  terms  of  a  series  not  wholly  unlike 
the  familiar  Taylor  series.  Table  I  is  presented  as  a  summary  of 
numerica]  differentiation  formulas  available  for  sets  of  range 
measurements  having  up  to  seven  readings.  The  time  interval  between 
readings  is  denoted  by  t. 
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Orbit  Determination  -  Geometric  Elements 
The  following  notation  is  introduced: 

r  =  jr|  =  vehicle  distance  from  earth  center 

y  =  earth  gravitational  constant 

a  =  semimajor  axis 

E  =  eccentric  anomaly 

e  =  eccentricity 

6  =  true  anomaly 


The  geometric  elements  are  determined  from  r,  r,  and  r  through 
three  relatively  simple  equations.  For  an  elliptic  orbit 


a 


U _ 

(y/r)-  £2  -  r  ir 


(3-2) 


e  cos  E  =  1  -  —  (3-3) 

a 

e  sin  E  =  ■  (3-4) 

/ya 

If  the  trajectory  is  hyperbolic,  a  is  negative.  The  equations  for  the 
geometric  elements  become 


a 


_ y _ 

(y/r)-  f2  -  r  r 


(3-5) 


e  cosh  F  =  1  -  — 

a 


(3-6) 
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esinhF  =  7=4  (3"7) 

The  hyperbolic  functions  result  from  the  substitution  of  iF  for  E  in 
the  elliptic  orbit  equations  (Ret  2:89). 

In  the  unlikely  event  the  orbit  is  exactly  parabolic,  then 

a  =  ®  i.e.  —  -  r2  -  r  r  =  0  (3-8) 

r 

e  =  1  (3-9) 

cos  0  =  ll-L  (3-10) 

U 

A  derivation  of  these  equations  is  presented  in  Appendix  B. 

There  are  no  approximations  involved  in  the  derivations .  The  numerical 
•  •> 

determination  of  r  and  r  involves  some  approximation. 

Example  Problems  -  Geometric  Element  Determination 

Example  Problem  3-1  -  Low  Eccentricity .  Liftoff  from  Cape 
Kennedy  was  normal,  but  before  parking  orbit  insertion,  the  space¬ 
craft  was  besieged  with  a  series  of  emergencies.  The  present  situation 
is  complete  communication  failure  coupled  with  unreliable  operation  of 
the  navigation  computer.  The  immediate  problem  is  determination  of 
the  safety  of  the  orbit. 

Solution: 

The  astronauts  obtain  the  following  altitudes  at  time 
intervals  of  five  minutes: 


42 


GA/AE/69-1 


hx  =  166.8  NM* 

h,  =  158.7  NM 
h3  =  147.8  NM 


The  expressions  for  r2  and  r2  become 


h3  -  h1  _  147.8  -  166.8 


2t 


2(5) 


=  -  1.9  NM/MIN 


r2  = 


h  +  h  -  2h  166.8  +  147. 8  -  2(158.7) 
13  2 


(5) (5) 


0.112  NM/MIN2 


For  the  earth,  p  and  /p"  are 


p  =  225.90258  x  106  NM3/MIN2 

/F  =  15.03005598  x  103  (NM/MIN2//2 


The  radial  distance  from  the  earth's  center  is  obtained  by  adding  the 
earth's  radius  to  the  altitude  reading  h2 


r2  =  h2  +  Re  =  158.7  +  3440.2  =  3599  NM 

The  semimajor  axis  a  is 


(225.9)  106 

(225.9)  106 
(3599) 


(1.9)2  +  (3599) (.112) 


*These  data  are  identical  to  Example  Problem  One  in  Ref  17:65. 
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and  using  a  slide  rule 

a  =  3575  NM  and  /a  =  59.79 


The  period  of  the  orbit  can  now  be  determined,  if  desired,  from 


and 


(4.180413)  10~4  a  /a 


P  =  (4.180)  10“4  (3575) (59.79) 


89.3  min 


The  eccentric  anomaly  at  is  obtained  from 

(3599) (-  1.9) (59.79) 

(15.03)  103  (3575  -  3599) 

228.58°  and 
sin  E2  =  -  0.7499 

The  eccentric  anomaly  is  greater  than  180°  since  the  radial  distance 
is  decreasing.  The  orbital  eccentricity  is  given  by 

r2  r2  _  (3599)  (-  1.9) 


then 


tan  E, 


r2  r2  /a 
/u  (a  -  r9) 


tan  E, 


1.133,  E2  = 


e 


f) T  /a  sin  E 


(15.03)  103  (59.79) (-  .7499) 


.01015 
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With  the  geometric  elements  determined,  the  perigee  altitude  is  given 
by 


hp  =  a(l  -  e)t-  R*  =  3575(1  -  .0101)  -  3440 
hp  =  99  NM 

Therefore,  the  safety  of  the  orbit  is  established  and  the  astronauts 
can  then  direct  their  attention  to  other  aspects  of  the  emergency. 

The  time  from  perigee  of  the  second  range  measurement  is  obtained 
from  the  "easy  solution"  of  Kepler’s  equation 

t  -  t  =  -  [E  -  e  sin  E]  (3-11) 

on 

where  n  =  2ir/P.  Then 

t2  -  tQ  =  — [(228.6)  (.01745)  +  (.0101)  (.7499)] 

and 

t2  -  tQ  =  56.8  rain 

The  iterative  computer  solution  to  this  example  problem  is 
presented  in  Ref  17:67.  Table  II  is  a  comparison  of  the  numerical 
differentiation  and  computer  solutions.  Note  that  identical  input 
data  were  used  and  these  data  were  assumed  perfect.  Therefore,  at 
this  point,  no  claim  can  be  made  regarding  the  ability  to  handle 
altitude  measurements  containing  normal  measurement  error. 
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Table  II 

Comparison  of  Computer  and  Numerical 
Differentiation  Solutions  of 

Example  Problem  3-1 

Orbital 

Parameter 

Numerical 

Differentiation 

Solution 

Computer 

Solution 

a 

3575  NM 

3575.9  NM 

e 

0.0101 

0.0100 

C2  “  Co 

56.8  min 

57.43  min 

Period 

89.3  min 

89 . 40  min 

Example  Problem  3-2.  The  translunar  injection  burn  has  been 
catastrophic.  In  addition  to  premature  engine  shutdown,  the  space¬ 
craft  tumbled  during  the  burn.  The  result  is  that  the  present 
spacecraft  orbit  is  completely  unknown.  Other  emergencies  dictate 
that  the  orbit  must  be  determined  manually. 

Solution; 

The  following  altitudes  are  obtained  at  30  minute 

intervals; 
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hj  -  7232.8  NM 

h2  =  9078.8  NM 

h3  =  10373.2  NM 

These  altitudes  are  known  to  be  slightly  in  error.  The  actual  alti¬ 
tudes  are: 

(h i > ACT  =  7195 * 833  M 

(h2)ACT  =  9031.767  NM 


(h3)ACT  =  10348.165  NM 

but  the  astronauts  must  use  their  only  information — the  measured  alti¬ 
tudes.  The  expressions  for  r  and  r  become 


2t  V 


10373.2  -  7232.8 
2(30) 


52.34  NM/MIN 


hj  +  h3-2h2  =  7232  +  10373  -  2(9078) 

=  -  0.6138  NM/MIN2 


For  the  earth,  p  and  \/\i  are 


p  =  225.90258  x  106  NM3/MIN2 

1 , 

=  15.03005  x  103  (NM3/MIN2)  2 


The  radial  distance  from  the  earth's  center  is  obtained  by  adding 
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the  earth's  radius  R  to  the  altitude  reading  h0 

©  c. 


r„  =  h,  +  R  =  9079  +  3440  =  12,519  NM 
z  © 


The  semimajor  axis  a  is 


(225 >9)  106 


(225.9)  106 
12519 


(52. 34)2  +  (12519) (.6138) 


a 


,(225^)_10i  _  9330  nm 

(22980) 


The  period  of  the  orbit  is 


P  = 


(4.180)  10~4  (9830) (99.15) 
(60) 


6.79  hr 


The  eccentric  anomaly  at  r2  is  obtained  from 


r0  r,  Ja. 

tan  E-  =  - ± - 

A*  (a  -  r2) 


(12519) (52.34) (99.15) 
(15.03)  103  (9830  -  12519) 


E2  =  121-88° 


The  eccentric  anomaly  is  less  than  180°  since  the  radial  distance  is 
increasing.  The  orbital  eccentricity  is  given  by 


r2  r2  (12519) (52.34) 

✓y"  ii  sin  E2  (15.03)  103  (99.15) (.8491) 
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The  actual  orbit,  the  computed  orbit  using  perfect  range  measurements, 
and  the  approximate  orbit  using  the  measured  values  of  range  are  shown 
in  Table  III. 


Table  III 

Comparison  of  Geometric  Parameters 
of  Example  Problem  3-2 

Actual 

Orbit 

Computed  Orbit 

Using  Perfect 

Range  Values 

Computed  Orbit 

Using  Measured 
Values  of  Range 

a 

10,000  NM 

10,017  NM 

9831.5  NM 

e 

0.5000 

0.4998 

0.5177 

E2 

119.629° 

119.364° 

121.88° 

Example  Problem  3-2  was  selected  to  illustrate  two  points: 

1.  Measurement  error  is  always  present. 

2.  Comparison  of  computer  versus  manual  solutions 
where  the  same  input  data  is  used  is  meaningless 
when  measurement  error  is  present.  The  important 
question  is:  How  far  is  the  manual  solution  from 
the  real  orbit? 
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The  adverse  effect  of  measurement  error  will  be  examined  further  in  the 
next  section. 

Summary 

The  geometric  elements  can  be  obtained  quite  easily  by  obtaining 
a  series  of  successive  range  measurements  at  equal  time  intervals. 
Numerical  differentiation  of  these  measurements  yields  the  radial 
velocity  and  radial  acceleration  at  one  of  the  range  readings.  With 
knowledge  of  range,  radial  velocity,  and  radial  acceleration,  the 
geometric  elements  a,  e,  and  E  are  readily  determined  from  fairly 
simple  equations.  A  block  diagram  outlining  the  steps  required  to 
obtain  the  geometric  parameters  is  shown  below. 
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Preliminary  Error  Analysis 

Introduction.  Before  any  orbit  determination  scheme  can  be 
presented  in  "flight  manual"  form,  it  must  be  determined  if  the 
method  will,  in  fact,  yield  the  desired  results.  Example  Problem  3-2 
illustrated  that  measurement  error  cannot  be  ignored. 

This  section  is  a  preliminary  error  analysis  of  the  numerical 
differentiation  method  of  orbit  determination.  The  word  "preliminary" 
is  used  since  a  statistical  approach  is  not  used.  Rather,  the  role 
of  measurement  error  is  studied  by  determining  the  worst  effect  it 
could  produce.  The  general  results  and  conclusions  from  either  type  of 
error  analysis  are  similar. 

The  trajectory  considered  is  that  which  tests  manual  navigation 
capability  to  the  utmost — the  high  eccentricity  ellipse.  Although 
this  will  highlight  the  shortcomings  of  the  numerical  differentiation 
scheme,  readers  should  deduce  that  as  eccentricity  is  lowered,  the 
effect  of  measurement  error  is  somewhat  reduced. 

The  major  error  sources  in  the  numerical  differentiation 
technique  of  orbit  determination  are: 

1.  numerical  differentiation  truncation  error 

2.  range  measurement  error 

3.  .timing  measurement  error. 

Each  of  these  error  sources  will  be  discussed. 

Truncation  Error  Effect.  Truncation  error  results  from  using 
numerical  differentiation  to  obtain  radial  velocity  and  radial 
acceleration.  Using  five  tabulations,  the  truncation  error  for  the 
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first  derivative  is  given  by 

y*  (x0)  -  p'  (x0)  =  ^  y<5>  (?)  (3-12) 


where 

y*  (xQ)  =  actual  derivative 

p’  (xQ)  =  approximate  derivative 

h  =  Xj  -  xQ  (the  tabular  interval) 

?  =  some  value  of  x,  generally 

indeterminate,  within  the  range 
of  values  of  x  under  consideration. 


A  derivation  of  this  equation  may  be  found  in  Ref  16:101.  Applied  to 
the  orbit  determination  problem 


r  -  r  =  r<5)  (?) 

3 (Actual)  3 (Approx)  30 


(3-13) 


for  five  range  measurements.  The  presence  of  t4  in  the  numerator 
suggests  that  t,  the  time  interval  between  readings  of  range,  should  be 
small  if  truncation  error  is  to  be  minimized.  A  similar  conclusion 
holds  for  the  second  derivative.  For  three  range  measurements,  the 
expression  for  truncation  error  for  r2  is 

r  -  r  =  r(4)  (?)  (3-14) 

2 (Actual)  2 (Approx)  12 


Since  determination  of  truncation  error  using  the  formulas  above  is 
only  an  estimate,  the  effect  of  truncation  error  on  the  orbit  determi¬ 
nation  scheme  was  obtained  in  a  slightly  different  manner.  Accurate 


52 


GA/AE/69-1 


values  of  range  at  equal  time  intervals  were  obtained  by  computer  for 
a  representative  high  eccentricity  trajectory.  Then  the  actual  orbital 
parameters  were  compared  with  those  obtained  from  the  formulas  Involving 
numerical  differentiation.  The  results  are  presented  graphically  in 
Figs.  12,  13,  and  14.  Note  that  in  all  cases,  truncation  error,  as 
expected,  increases  as  the  time  interval  between  range  measurements 
is  increased. 

Range  and  Timing  Measurement  Error.  Inaccurate  values  of  range  are 
a  major  source  of  error  in  the  numerical  determination  of  r  and  ir.  In 
addition  to  its  effect  through  the  numerical  differentiation  process, 
an  inaccurate  range  measurement  also  directly  affects  the  values  of  the 
orbital  elements  through  the  value  of  r  itself  (since  the  equations 

for  a,  e,  and  E  each  contain  an  explicit  value  of  range).  Therefore, 

* 

it  is  necessary  to  consider  the  effect  of  errors  in  range  measurement 
on  the  orbit  determination  method  being  considered.  Another  error 
source  is  that  of  timing  between  measurements.  The  following  presenta¬ 
tion  will  begin  with  a  brief  discussion  of  how  range  measurement  error 
affects  the  numerical  differentiation  process.  Then  the  effect  of 
range  and  timing  measurement  error  will  be  investigated  for  each  of 
the  orbital  parameters  being  considered. 

Measurement  Error  Effect  -  Numerical  Differentiation.  A  major 
step  in  the  orbit  determination  method  being  considered  is  the 
numerical  determination  of  r  and  r  from  a  set  of  range  measurements. 

A  brief  discussion  of  how  inaccurate  range  measurements  can  affect 
these  derivatives  is  thus  in  order.  The  cumulative  effect  of 
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Fig.  12.  Effect  of  Truncation  Error  on  Semimajor  Axis 


Effect  of  Truncation  Error  on  Eccentricity 


Calculations  are  for  Ellipse 
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Pig.  14.  Effect  of  Truncation  Error  on  Eccentric  Anomaly 
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measurement  errors  on  the  final  output — the  orbital  parameters— will 
be  treated  in  a  later  section. 

Assume  that  e^,  e 2,  —  e5  represent  the  errors  associated  with 

each  measurement  of  range  and  that  Rj,  R2,  -  R5  are  the  actual 

range  values  that  should  be  observed.  Then 

Rj  =  rj  +  t}  (3-15) 

where  r,  represents  the  observed  range.  Then  using  this  definition 

4 

in  the  five-measurement  differentiation  formula,  the  error  in  output 
due  to  input  inaccuracy  is 

(Rj  -  Rs)  -  8(R2  -  R4)  (rx  -  r5)  -  8(r2  -  r4) 

—  — 

(Sl  -  e5)  -  8(e2  -  e4) 

12t  (3-16) 

If  gj,  .  .  .  eg  do  not  exceed  emax  in  magnitude,  the  output  error  is, 
at  worst 


Maximum  output  error  in  I8e 

•  max 

r3  due  inaccurate  data  =  — 


(3-17) 


A  small  time  interval  t  generally  associated  with  high  accuracy  (i.e., 
small  truncation  error)  thus  magnifies  the  effect  of  measurement 
error — many  times  making  the  trouble  of  minimizing  truncation  error 
pointless. 

By  similar  analysis,  the  maximum  error  for  the  second  derivative 
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(using  five  measurements  again)  is 

Maximum  output  error  in 
r  due  inaccurate  data  =  - 222. 

3  3t2  (3-18) 

This  shot^s  that  higher  derivatives  become  increasingly  inaccurate  as 
the  time  interval  is  reduced.  Note,  however,  the  particular  arrange¬ 
ment  of  measurement  errors  such  that  r  and  r  both  exhibit  their 
maximum  errors  at  the  same  time  cannot  occur.  Nevertheless,  the  adverse 
effect  of  measurement  error  on  numerical  differentiation  is  severe. 

Long  time  intervals  between  measurements  act  to  reduce  this  adverse 
effect. 

The  cumulative  effect  of  measurement  errors  (thus  including 
numerical  differentiation)  on  the  final  output~the  orbital  parameters — 
is  the  topic  of  the  next  section. 

Estimate  of  the  Error  in  Determining  a_,  e_,  and  E.  The  orbit 
determination  scheme,  as  developed  so  far,  involves  the  measurement  of 
several  values  of  range  at  equal  time  intervals.  These  observations 
are  numerically  differentiated  to  yield  r  and  r  for  the  central 
measurement  and  then  the  geometric  orbital  elements  are  obtained  from 
the  three  equations 


tan  E 

e  sin  E 


Jv  (a  -  r) 


r  r 

/T  /a 


(3-2) 

(3-19) 

(3-4) 
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The  error  arising  from  numerical  differentiation  has  been  briefly 
discussed.  It  remains  to  determine  the  overall  effect  of  measure¬ 
ment  errors  on  the  values  of  a,  e,  and  E  obtained  through  the  use  of 
these  approximate  data  in  the  three  equations  listed  above. 

In  general,  the  maximum  error  Eu  in  the  value  of  a  function  is 
given  by 


E 


u 


3f 

V  4- 

3f 

3x 

Ea  + 

3y 

x  =  a  x  =  a 

y  =  b  y  =  b 


(3-20) 


where 


E 


a 


Max  error  in  function  argument  a 


Eb 


Max  error  in  function  argument  b 


A  derivation  of  this  expression  is  shown  in  Appendix  C.  This 
expression  will  be  used  to  examine  the  possible  errors  in  orbital 
parameters  due  to  measurement  errors. 

Error  in  Major  Axis  Determination.  For  the  following  analysis, 
it  is  assumed  that  truncation  error  in  the  determination  of  r  and  r 
is  negligible  (a  valid  assumption  if  a  sufficient  number  of  range 
measurements  are  used) .  Then  the  maximum  error  in  a,  due  to  range 
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meisurement  and  timing  error,  is 


E 


a 


= 

3a 

E_  + 

3a 

E„  +  3a  1  E  + 

3a 

3ri 

rl 

3r2 

r2  3r3!  r3 

3r4 

+ 

3a 

E„  + 

3a 

3r5 

r5 

3t 

(3-21) 


for  five  measurements  of  range.  Assuming  Er  =  E„  .  .  .  E„  then 

rl  r2  r5 


K 

3a 

+ 

3a 

! 

+ 

3a 

4* 

3a 

+ 

3a 

3rl 

3r2 

3r3 

3r4 

3r5 

+ 


3a 

3t 


(3-22) 


The  equation  for  a  is 


a 


(3-2) 


and  for  five  measurements. 


3a  = 
3x 


+ .. 

+  r3  3x  /  (3-23) 
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where  x  is  taken  to  successively  represent  r  ,  r  ,  .  .  .  r  ,  and  t. 

The  partial  derivatives  of  r  and  r  may  be  obtained  by  differentiating 
the  approximate  differentiation  formulas.  For  convenience,  let  Ca  r 
represent  the  coefficient  of  Er  and  Cajt  the  coefficient  of  Et.  Then 

Ea  "  ca,r  Er  +  Ca,t  Et  (3-24) 

where 


3a 

+ 

3a 

+ 

3a 

+ 

3a 

9rl 

3r2 

3r3 

3r4 

9r5 

(3-25) 


and  Ca  t  =  1 3 a/3 1 1 .  Then  by  examining  the  behavior  of  Ca^r  and 

C  ^  the  steps  required  to  minimize  the  effect  of  range  measurement 
a>  l 

error  can  be  determined. 

A  computer  program  was  written  to  evaluate  the  coefficients  C 

a,r 

and  Ca>t  for  various  values  of  t  and  for  a  high  eccentricity  trajectory 
similar  to  the  initial  translunar  coast  portion  of  a  moon  mission.  The 
results  are  presented  graphically  in  Figs.  15  and  16.  It  should  be 
noted  that  Figs.  15  and  16  yield  the  behavior  of  the  maximum  error — not 
the  actual  or  rms  error. 

Error  in  Eccentric  Anomaly  Determination.  Eccentric  anomaly  error 
can  be  investigated  using  the  same  procedure  as  was  followed  in  examining 
the  error  in  major  axis  determination.  The  equation  for  eccentric 
anomaly  is 

E  =  tan  _1  .  (3-19) 

AT  (a  -  r) 


\ 
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and  substituting  for  a,  this  equation  becomes 


E  =  tan' 


.  (V  ♦  2  ..  \ 

-1  r\r  -  r  -  r  r) 

(r2  +  r  r) 


/2 


(3-26) 


For  simplication,  let 


u  =  r 


/  2 


(3-27) 


w 


=  (r2  +  r  r) 


then 


E  =  tan-1  uv 
w 


(3-28) 


Taking  partial  derivatives  as  before. 


9E  _ 
3x 


1  + 


u2  v2 
w2 


w  fu  +  v  -  U  V 
\  3x  3 x  / 


3w 

3x 


w 


(3-29) 


(the  independent  variable  x  is  again  taken  to  successively  represent 
ri»  r2»  .  .  .  r^,  t) .  The  maximum  error  in  E  (Eg)  due  to  range 
measurement  and  timing  errors  is  then 


3E 

E_  + 

3E 

+  •  •  v  -f* 

3E 

E,.  + 

3E 

9rl 

rl 

3r2 

r2 

9rn 

rn 

3t 

(3-30) 
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where  n  denotes  the  number  of  range  measurements  taken.  Again  assuming 

E*.  —  Er  88  •  •  .  E.. 

1  2  rn 


(3-31) 


eE  =  cE,r  Er  +  CE,t  Et  (3"32) 

where  Cg  r  and  Cgjt  represent  the  coefficients  of  Er  and  Et  respec¬ 
tively.  Again,  Cg  r  and  Cg  t  can  be  examined  and  the  steps  required 
to  minimize  the  effect  of  measurement  errors  can  be  determined. 

Figures  17  and  18  are  the  graphical  results  of  a  computer 
program  written  to  evaluate  Cg  r  and  Cg>t  for  a  typical  trajectory. 

Error  in  Eccentricity  Determination.  The  effect  of  measurement 
error  on  the  calculation  of  eccentricity  can  be  examined  in  the  same 
manner  as  for  the  major  axis  and  eccentric  anomaly.  The  expression 
for  eccentricity  is 


e 


a  -  r 
a  cos  E 


Differentiating  as  before 


(3-33) 


9e 

9x 


(a  cos  E)|jf  -  |f[-  (a  -  r)[(-  a  sin  E)  ||  +  (cos  E)  || 


a2  cos2  E 


(3-34) 
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Fig.  17.  Effect  of  Range  Measurement  Error  on  Eccentric  Anomaly 
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Fig.  18.  Effect  of  Timing  Error  on  Eccentric  Anomaly 
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where  the  independent  variable  x  is  again  taken  successively  as  r^ 
.  .  .  rn,  and  t.  The  maximum  error  in  e  is  then  approximately 


E 


e 


De 

+ 

3e 

+  . . .  + 

3e 

Er  + 

3e 

3rl 

3r2 

3rn 

I 

3t 

(3-35) 


or 


Ee  “  Ce,r  Er  +  ce,t  Et 


(3-36) 


where  Ee  is  the  maximum  error  in  eccentricity  and  Ce>r  and  Ce>t 
represent  the  coefficients  of  range  measurement  error  Er  and  timing 
error  Et  respectively. 

Figures  19  and  20  are  the  graphical  results  of  a  computer  program 
written  to  evaluate  Ce>r  and  Ce>t  for  a  typical  trajectory. 

Accuracy  of  Range  Measurement.  Measurement  of  the  earth's  angular 
diameter  is  assumed  the  primary  means  of  determining  distance  from  the 
earth.  The  geometry  of  the  measurement  is  shown  in  Fig.  21. 


Trajectory 


Fig.  21 

Range  Measurement  Geometry 
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on  Eccentricity 


Fig.  20.  Effect  of  Timing  Error  on  Eccentricity 
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The  distance  r  is  given  by  the  equation 


r 


(3-37) 


where  R  is  the  radius  of  the  earth  and  A  is  the  earth’s  subtended 
angle  as  measured  from  the  spacecraft.  Then 


-  Rffl  cos  A/ 2 

dr  - - - - dA  (3-38) 

2  sin2  A/2 

As  distance  from  the  earth  increases,  the  angle  A  becomes  smaller — 
the  net  result  for  a  given  measurement  error  dA  is  a  decrease  in 
accuracy  as  distance  from  the  earth  increases.  Figure  22,  from 
"Space  Position  Fixing  Techniques,"  (Ref  22:318),  is  a  graph  showing 
the  error  in  distance  per  minute  error  in  A  as  a  function  of  distance 
from  the  earth.  Note  that  statute,  not  nautical,  miles  are  presented 
on  the  graph. 

With  proper  training,  practice,  and  instrument  calibration, 
sighting  measurement  error  can  be  held  appreciably  below  one  arc 
minute  (Ref  17:5).  Therefore,  Fig.  22  can  be  considered  to  represent 
the  upper  bound  on  range  measurement  error. 

When  only  a  portion  of  the  earth’s  horizon  is  visible,  the 
stadimeter,  a  device  which  measures  the  curvature  of  the  horizon,  can 
be  used  to  determine  altitude  (Ref  17:12).  The  stadimetar  is  designed 
to  measure  altitudes  below  2500  NM.  The  trajectories  considered  in 
this  thesis  are  generally  above  this  altitude. 
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Another  unfortunate  situation  is  also  apparent.  As  distance  from 
the  earth  increases,  the  accuracy  of  range  measurement  decreases.  At 
the  same  time,  the  effect  of  measurement  error  on  some  orbital  param¬ 
eters  becomes  more  severe.  It  is  essential,  therefore,  to  begin  orbit 
determination  as  near  the  earth  as  possible. 

It  can  be  safely  assumed  that  timing  error  can  be  held  to  less 
than  five  seconds.  For  the  high  eccentricity  orbits  examined  and  a 
one  hour  time  interval  between  readings,  five  seconds  is  capable  of 
causing  a  maximum  error  in  major  axis  of  less  than  100  NM.  Timing 
e. rc  5  thus  not  a  major  drawback  in  orbit  determination  using  this 
method. 

Conclusions.  The  conclusions  derived  from  the  preliminary  error 
analysis  are: 

1.  Decreasing  truncation  error  by  taking  five  range 
measurements  instead  of  three  is  pointless  unless 
extremely  accurate  range  measurements  are  available. 

2.  For  high  eccentricity  lunar  trajectories,  a  minimum 
time  interval  of  five  hours  between  measurements  is 
required  to  determine  the  approximate  orbit  (thus, 
ten  hours  are  required  for  orbit  determination) . 

3.  Every  effort  must  be  directed  toward  the  highest 
attainable  accuracy  in  range  measurements. 

4.  Timing  error  of  five  seconds  or  less  is  insignificant. 

5.  Orbit  determination  must  commence  as  near  the  earth 
as  possible. 
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6.  A  supplementary  correction  scheme  may  be  required  to 
further  improve  the  accuracy  of  the  orbital  elements. 

7.  Alternate,  more  accurate,  means  of  obtaining  range, 
range  rate,  and  radial  acceleration  should  be  explored. 

8.  The  use  of  graphs,  slide  rule,  and/or  uninterpolated 
tabular  data  is  entirely  consistent  with  the  degree  of 
accuracy  attainable  using  this  orbit  determination 
method. 
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IV.  Differential  Correction  of  Orbital  Elements 


Introduction 

The  orbital  elements  obtained  by  any  means  or  technique  will  con¬ 
tain  a  certain  amount  of  error.  For  example,  in  the  orbit  determina¬ 
tion  method  of  Chapter  III  (orbit  determination  using  range,  range  rate, 
and  radial  acceleration) ,  the  orbital  elements  were  found  to  be 
increasingly  in  error  as  the  time  interval  between  range  measurements 
was  reduced.  This  section  is  a  brief  discussion  of  a  means  to 
correct  the  geometric  elements  of  a  preliminary  orbit  to  obtain  more 
precise  values  of  a,  e,  and  M5.  The  method,  however,  requires  the  use 
of  a  digital  or  mechanical  analog  computer  because  of  the  complexity 
of  the  computations.  The  inclusion  of  such  a  method  in  a  thesis 
titled  "Manual  Astronaut  Navigation"  can  be  justified  by  supposing 
the  eventual  development  of  a  hand-held  back-up  computer. 

Differential  Correction  of  Geometric  Orbit  Parameters 

Differential  correction  of  orbits  is  discussed  in  most  books 
covering  methods  of  orbit  determination.  For  example,  see  Ref  6:233 
or  Ref  3:77.  Most  applications  of  differential  corrections  in  the 
field  of  celestial  mechanics,  however,  require  the  simultaneous 
solution  of  six  equations  for  the  corrections  to  the  six  independent 
elements  characterizing  an  orbit.  The  computation  complexity  can 
be  reduced  somewhat  by  the  technique  of  dividing  the  six  parameters 
into  two  groups — three  parameters  (a,  e,  M0)  describing  the  size 
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and  shape  of  the  orbit,  and  the  remaining  three  (i,  Q,  w)  describing 
the  orbit  orientation.  Thus,  instead  of  six  equations  with  six  unknown 
corrections,  there  are  two  groups  of  three  equations  in  three  unknowns. 

Differential  Correction  Equations  -  Geometric  Elements 
The  conic  section  equation  is 

r  =  a(l  -  e  cos  E)  (4-1) 

The  total  differential  dr  is  then 

dr  =  If  da  +  lf  de  +  H  dE  (4-2) 

Then,  approximately, 

ir  -  If  4a  +  If  le  +  ft  4E  <4-3> 

The  quantity  Ar  can  be  obtained  by  comparing  the  measured  value  of  r 
at  a  certain  time  with  the  corresponding  preliminary  reference  orbit 
value  at  the  same  time.  That  ia, 

Ar  =  robserved  "  preliminary  orbit  (4-4) 

The  partial  derivatives  are  obtained  from  the  conic  section  equation 
and  from  Kepler’s  equation.  Three  values  of  Ar,  obtained  at  different 
times,  yield  three  equations  which  can  be  solved  for  the  corrections 
to  be  applied  to  the  preliminary  orbital  elements.  The  partial 
derivatives  with  respect  to  a,  e,  and  E  are,  from  the  conic  section 
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The  differential  change  in  mean  anomaly  is  then 

dM  =  dM0  +  dn  (t  -  t0) 

Substituting  for  dn, 

dM  =  dM0  -  |  n  |2.  (t  -  t0) 

Since 


n(t  -  t0)  =  M  -  M0 

dM  =  dM0  -  (m  -  M0) 

From  Kepler's  equation,  M  =  E  -  e  sin  E 

dM  =  dE  -  e  cos  E  dE  -  sin  E  de 


or 


dM  =  (1  -  e  cos  E)  dE  -  (sin  E)  de 

and  substituting 

—  =  1  -  e  cos  E, 

a 

dM  =  L  dE  -  (sin  E)  de 
a 


and  solving  for  dE 

dE 


fl? 


dM  +  (sin  E)  de 


0 


(4-12) 


(4-13) 


(4-14) 

(4-15) 


(4-16) 


(4-17) 


(4-1) 

(4-18) 


(4-19) 
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Now  substituting  for  AM 

dE  .=  ^  |~dM0  -  -|~  ^  (M  -  Ho)]  h  ^  (sin  E)  de  (4-20) 

Approximating  the  differentials  with  the  incremental  corrections  Aa, 

Ae,  and  AM0 

Ar  -  —  Aa  -  a(cos  E)  Ae  +  a  e(sin  E)  AE  (4-21) 

a 

Substituting  for  AE  and  simplifying  results  in 

Ar  =  [l  -  3(M  -  M0)  a  e  sin  E~]  Aa 

La  2r  J 

+  sin2  E  -  a  cos  eJ  Ae 

+  sin  eJ  AMd  (4-22) 

The  bracketed  coefficients  of  Aa,  Ae,  and  AMQ  are  to  be  evaluated  on 
the  preliminary  orbit  at  the  times  corresponding  to  the  range  measure¬ 
ments. 

Discussion 

Three  measurements  of  Ar  obtained  at  different  times  provides 
three  equations  in  the  unknowns  Aa,  Ae,  and  AMQ.  The  equations  can 
be  inverted  to  yield  the  values  Aa,  Ae,  and  AM0.  When  the  number  of 
readings  exceeds  three,  the  number  of  unknowns,  the  most  probable 
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values  of  the  corrections  can  be  found  by  least  squares — but  this 
greatly  increases  the  computational  load. 

Assuming  the  availability  of  a  back-up  computer,  the  geometric 
parameters  can  be  obtained  directly  from  Eq  (4-22)  by  estimating 
values  of  a,  e,  and  M0  and  then  correcting  these  assumed  values  in  an 
iterative  process  until  the  orbital  elements  approach  the  desired 
accuracy.  However,  such  a  procedure  is  of  questionable  value  since 
any  of  a  number  of  iterative  techniques  could  be  developed  to  evaluate 
the  geometric  parameters.  Also,  a  number  of  iterative  techniques 
are  available  which  yield  all  six  orbital  elements — so  why  iterate  to 
obtain  only  the  geometric  elements? 

A  computer  program  was  written  to  test  the  effectiveness  of  the 
differential  correction  procedure.  Data  generated  from  previous 
computer  programs  were  used  to  improve  an  orbit  known  to  be  in  error. 
Only  one  correction  was  applied  to  the  erroneous  orbit,  hence,  the 
final  orbit  is  still  not  precise — but  further  corrections  could  be 
applied  to  improve  the  final  orbit.  Figure  23  is  a  block  diagram 
presentation  of  the  results. 

Conclusion 

The  computer  makes  simple  work  of  the  computations  involved  in 
this  differential  correction  scheme.  It  is  out  of  the  question  to 
consider  accomplishing  the  calculations  manually.  If  a  small 
back-up  computer  is  available,  the  differential  correction  equations 
may  be  of  great  value  because: 

1.  The  geometric  parameters  can  be  refined  to  more 


accurate  values 
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2.  The  capability  exists  to  handle  redundant  data. 

The  differential  correction  method  will  not  be  considered  further 
in  this  thesis.  However,  the  recommendation  section  will  contain  a 
suggestion  for  continued  study  in  this  arc  a. 
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V.  Orbit  Determination  Using  Position  Fixes 


Introduction 

Two  complete  three-dimensional  fixes  in  space  and  knowledge  of  the 
time  between  the  two  fixes  is  sufficient  for  orbit  determination  using 
Lambert's  theorem  (Ref  4:70).  Lambert's  theorem  requires  iterative 
solution  and  its  use  is  impractical  unless  a  back-up  computer  is  avail¬ 
able.  By  taking  another  position  fix,  thus  providing  three  position 
fixes  in  space,  any  of  a  number  of  orbit  determination  methods  may  be 
considered  (see  Ref  3:Ch  1).  This  section  will  discuss  first  a  simpli¬ 
fied  and  then  a  general  method  of  obtaining  three-dimensional  fixes  in 
space  using  the  hand-held  space  sextant.  These  fixes  will  be  used  as 
input  to  orbit  determination  methods  which  will  yield  both  the  geometric 
and  orientation  orbital  elements.  Unless  otherwise  specified,  the 
geocentric  equatorial  coordinate  system  is  used  in  this  chapter  (and 
succeeding  chapters) . 

Position  Vector  From  Polaris,  Earth,  2nd  Star 

Introduction.  The  navigation  fix  in  space  is  discussed  in 
Ref  4:221.  In  general,  a  three-dimensional  fix  requires  two  star-earth 
horizon  measurements  coupled  with  a  means  to  determine  radial  distance 
from  the  earth.  The  equations  are  nonlinear  and  impractical  for  manual 
computations.  Considerable  simplification  results,  however,  if  one  of 
the  selected  stars  is  Polaris,  the  North  Star.  The  development  to 
follow  is  due,  in  part,  to  Hr.  L.  C.  Ragland,  of  TRW  Systems  (Ref  14). 
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where  Rffi  is  the  earth's  radius. 

Right  Ascension  Determination.  Figure  25  shows  the  geometrical 
arrangement  for  finding  the  right  ascension  of  the  vehicle,  ot^,. 


(f)  =  Star-Earth-Vehicle  Angle 
<5V  =  Vehicle  declination 

Aa  =  Difference  in  Right  Ascensions  of  Star  and  Vehicle 
otg> Ss  =  Right  Ascension  and  Declination  of  Star 


Fig.  25 

Geometry  for  Vehicle  Right  Ascension 


! 


! 


From  the  spherical  triangle  star-Polaris-vehicle, 
cos  v  =  cos  (90  -  6g)  cos  (90  -  Sv) 

+  sin  (90  -  6S )  sin  (90  -  6V)  cos  Aa 


or 


cos  iJj  =  sin  6S  sin  6V  +  cos  6S  cos  6V  cos  Aa 


(5-3) 


(5-4) 
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then 


cos  Ao 


cos  4*  “  sin  6S  sin  6V 
cos  <Ss  cos  6V 


(5 


The  terms  are  defined  on  Fig.  25.  Let  G  be  a  unit  vector  normal  to 
the  star-earth-vehicle  plane.  An  observer  looking  down  from  the  tip 
of  this  vector  toward  the  earth  would  see  the  picture  shown  in 
Fig.  26.  Again,  since  the  lines  of  sight  to  the  star  are  parallel, 

tl>  +  Yg  +  A/2  =  180°  (5 


or 


\p  =  180°  -  (ys  +  A/2) 


(5 


Fig.  26 

Geometry  of  Star-Earth-Vehicle  Plane 

OO 
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Substituting  ip  into  the  expression  for  Act  results  in 

-  cos  (ys  +  A/2)  -  sin  6S  sin  6y 

cos  Aa  =  -  (5-8) 

cos  6s  cos  dv 

from  which  Aa  can  be  determined.  Note  the  resulting  simplification 
if  a  star  can  be  chosen  with  6S  =  0.  With  Aa  known  the  right 
ascension  of  the  vehicle  is  given  by 

av  =  as  ±  Aa  (5-9) 

where  a  simple  "right-hand-rule"  determines  the  sign.  If  the  origin, 
for  the  moment,  is  considered  at  the  vehicle,  then  the  sign  is  plus 
if  the  general  direction  of  istar  x  Ip0iaris  is  toward  the  earth. 

Istar  and  ^Polaris  represent  unit  vectors  in  the  directions  of  the 
star  and  Polaris  respectively. 

When  a  star  with  zero  declination  can  be  chosen  Eq  (5-8) 
becomes 

cos  to  -  -  £°S-lls+Am.  (5-3.0) 

cos  6V 

A  number  of  stars  are  very  close  to  the  celestial  equator.  (Most  of 
these  are  closer  to  the  equator  than  Polaris  is  to  the  pole  position.) 
Table  VI  is  a  tabulation  of  some  stars  with  near  zero  declination. 

An  easily  identifiable  star  is  6  Ori  (Mintak^)  which  is  in  the  belt 
of  Orion.  With  r,  cty,  and  ov  known,  the  coordinates  of  the  vehicle 
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Table  VI 

Selected  Stars  Near  Celestial  Equator 

Star 

IAU 

Right 

Name 

Abbreviation 

Constellation 

Ascension 

Declination 

Kin taka 

S  Ori 

Orion 

84.249° 

-  0.321° 

— 

6  Mon 

Monoceros 

107.557° 

-  0.438° 

— 

6  Vir 

Virgo 

203.265° 

-  0.433° 

Sadalmelik 

a  Aqr 

Aquarius 

331.035“ 

-  0.476° 

are  given  by 


x  =  r  cos  Sv  cos  av 


y  =  r  cos  5V  sin 


(5-11) 


z  =  r  sin  6, 


Conclusion.  The  simplified  three-dimensional  fix  presented  in 
this  section  has  as  its  basis  the  fact  that  Polaris  is  very  near  the 
north  celestial  pole.  Of  course,  Polaris  is  not  exactly  at  the  pole 
position  and  this  discrepancy  must  be  examined.  But  even  if  a 
detailed  error  analysis  proves  the  method  unsatisfactory,  the  sim¬ 
plified  three-dimensional  fix  can  be  of  great  value  to  the  astronaut. 
Op  the  earth,  Polaris  is  extremely  useful  as  an  instant  cross-check 
of  heading  and  latitude.  Its  analogous  role  in  space  navigation  is 
as  a  source  of  vehicle  declination.  Using  another  star  on  the 
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celestial  equator,  the  astronaut  can  quickly  determine  the  right 
ascension  of  the  vehicle.  Therefore,  even  if  eventually  deemed 
unsatisfactory  for  manual  navigation,  the  Polaris  position  fix  will 
serve  as  a  quick  cross-check  of  more  sophisticated  fixing  means. 

General  Solution  for  Position  Vector 

Introduction.  The  Polaris-earth-star  measurement  scheme  is  the 
simplest  method  known  to  the  authors  of.  obtaining  a  three-dimensional 
fix.  This  section  will  attack  the  more  difficult  general  problem 
when,  for  some  reason,  Polaris  is  not  observed. 

The  general  formulation  of  the  three-dimensional  fix  is  given 
in  Ref  4:Ch  7.  In  general,  the  three-dimensional  fix,  using  the 
minimum  number  of  measurements,  requires  the  simultaneous  solution  of 
the  following  nonlinear  equations 

=  -  r  cos  (y1  +  A/2) 

=  -  r  cos  (y2  +  A/2)  (5-12) 

=  r 

__  A  A 

where  r  is  the  position  vector  of  the  vehicle,  Ij  and  I2  are  unit 
vectors  in  the  directions  of  two  selected  stars,  Yj  2»  t*ie  star- 
vehicle-earth  horizon  angles,  and  A,  the  earth's  subtended  angle. 

As  mentioned  previously,  the  manual  solution  of  these  equations  for 
the  components  of  r  is  impractical.  Some  simplification  results, 
however,  if  an  additional  star  sighting  is  made.  • 


91 


GA/AE/69-1 


The  Three  Star  Fix.  The  unit  position  vector  of  the  vehicle  is 


obtained  fairly  simply  by  making  three  star-earth  horizon  measure¬ 
ments.  The  set  of  linear  equations 


r  •  I1  =  -  cos  (Yj  +  A/2) 


r  •  I2  =  -  cos  (y2  +  A/2) 


(5-13) 


r  •  Ig  =  -  cos  (y3  +  A/2) 


must  be  solved  for  the  components  of  the  unit  vector  r.  The  solution 

A 

for  rx,  the  x  component  of  the  unit  vector  r,  is 


rx  =  -  cos  (Yj  +  A/2) 


^Xv2  1zZ  ~  Xv3 


,  (Iv,  I,,  -  Iv,  I?1) 

+  cos  (y0  +  A/2)  — ^2 - - - - — 

2  A 


(5-14) 


(I  I  -  I  I  ) 
-  cos  (y3  +  A/2)  -^1-2-2 - Z2_JLLl 


where  A  represents  the  determinant  of  coefficients 


Ixl  *yl  Izl 


A  “  Ix2  Xy2  “z2 


(5-14a) 


*X3  Iy3  *z3 
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The  solutions  for  ty  and  rz  are,  of  course,  similar.  In  general,  the 
solutions  take  the  form 

rx  =  Axl  cos  (Yj  +  A/2)  +  Bx2  coj;  (y2  +  A/.2) 

+  Cx3  cos  + 

ry  =  Ayi  cos  (yx  +  a/2)  +  By2  cos  (y2  +  A/2) 

+  Cy3  cos  (y3  +  A/2) 

rz  =  Azl  cos  (Yj  +  A/2)  +  Bg2  cos  (y2  +  A/2) 

+  Cz3  cos  (y3  +  A/2)  (5-15). 

with  the  terms  defined  as  follows: 


A  .  = 

xl 

(Iys 

Iz2 

Iy2 

V'4 

y  ‘ 

^Xl 

Iz2 

"  Ix2 

v'4- 

11 

PQ 

(Iyi 

Xz3  - 

Iy3 

V'4 

Azl 

^x3 

\z 

"  Ix2 

V/4 

x° 

OJ 

II 

«y2 

Izl 

V 

Lz2i  ^ 

BZ2  * 

^xl 

Ty3 

Xx3 

Ayi 

V* 

*z3  “ 

Tx3 

!«)/« 

Cz3  - 

<Tx2 

Tyl 

-Txl 

v/4 

By2  " 

0*3 

Jzl  - 

*xl 

W'4 

(5-16) 

A  ~ 

I  f 

’i 

I 

-  I 

I  1  -  I 

r  I 

I 

-  I  „ 

I  *1 

xl  [ 

_  y2 

z3 

ya 

z2_|  X2 

_yi 

z3 

y3 

zlj 

+  **3  [V  ^  '  V  ^l] 
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It  is  a  simple  computer  problem  to  tabulate  the  above  constants  for  a 
considerable  number  of  three-star  sets.  The  astronaut  is  thus 
provided  some  semblance  of  general  capability  to  manually  compute  his 
position  in  space.  Once  the  unit  vector  r  is  obtained,  the  vehicle 
position  vector  is  given  by 

f  =  r  r  (5-17) 

where  r  =  |r|  and  is  obtained  by  semi-diameter  or  other  means. 

Manual  calculation  of  position  using  this  procedure  will  require 
a  publication  listing  the  constants  AX1  through  Cz3  for  selected 
three-star  sets.*  Table  VII  shows  a  typical  format  for  such  a 
publication.  Although  the  requirement  exists  for  a  third  star 
sighting  (compared  to  the  nonlinear  three-dimensional  fix  equations) , 
the  additional  measurement  removes  the  nonlinearity  plus  the  possible 
ambiguity  in  the  direction  of  r  (Ref  4i221).  Following  Table  VII  is 
a  sample  problem  illustrating  the  manual  position  fix. 


*  If  25  stars  are  established  as  "space  navigation  stars,"  there  are 
2,300  possible  three-star  sets.  Assuming  20  sets  per  page,  the 
publication  would  require  only  115  pages  to  cover  all  possibilities. 
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Table  VII 

Example  of  Three-Star  Position  Fix  Data 

- 

^1 

Bx2 

Cx3 

Star  Set 

V 

By2 

cy3 

Azl 

Bz2 

Cz3 

1  Dubhe 

-.90368847 

-.29882173 

-1.68227200 

2  Betelgeuse 

.23797789 

-.99759279 

.16392111 

3  Fomalhaut 

-1.63183470 

-.02032258 

-.89050598 

1  Pollux 

3.16269770 

-5.84253850 

-3.38404430 

2  Betelgeuse 

-.19226972 

-.28394138 

.58373726 

3  Rasalague 

.78242022 

-4.28971420 

-3.74551760 

1  Capeila 

-.32702313 

.28773730 

-.98883642 

2  Acrux 

-1.67469260 

-1.15633370 

-.35108396 

3  Fomalhaut 

.25868869 

1.04071820 

.52196385 

1  Acrux 

.22619851 

.84496636 

-.51440839 

2  Arcturus 

.28928391 

.35811369 

.90990474 

3  Altair 

.99216655 

-.44923096 

.21120275 

1  Dubhe 

.25166101 

3.02868280 

2.96921080 

2  Sirus 

-.24731313 

-.00109432 

1.00028020 

3  Rasalague 

-.97028150 

1.56380340 

1.39883810 

1  Capeila 

1.76155570 

-1.74660870 

2.98233530 

2  Canopus 

-.93319008 

-.62009385 

-.19967779 

3  Acrux 

-.84378145 

.92195900 

-.3816431* 

1  Dubhe 

-1.69239350 

-.52695925 

-2.16196100 

2  Canopus 

-2.39505190 

-1.75921190 

-1.43748310 

3  Fomalhaut 

-I. 68600160 

-.03619068 

-.92345019 

1  Acrux 

.77487909 

.90722594 

-.64576203 

2  Arcturus 

-.68124277 

.24798659 

1.14224790 

3  Shaula 

.76689250 

-.47479314 

.26513315 
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Example  Problem  5-1  -  Manual  Three-Dimensional  Fix.  Communication 
failure  occurred  during  the  translunar  injection  burn.  Although  the 
navigation  computer  seems  to  be  working  properly,  the  astronauts  suspect 
the  cause  of  communication  failure  may  have  also  affected  the  computer. 
The  astronauts  decide  to  obtain  a  manual  fix  and  cross-check  the 
manual  position  with  the  computer  read-out.  If  the  positions  cross-check 
reasonably  well,  the  computer  will  be  relied  upon  for  further  naviga¬ 
tion.  Looking  toward  the  earth,  the  astronauts  can  identify  the 
following  navigation  stars: 

Acrux 

Arcturus 

Altair 

The  following  angular  measurements  are  obtained  at  TLI  (Trans  Lunar 


Injection)  plus  three  hours: 

Acrux-vehicle-earth  horizon  -  43.970° 

Arcturus-vehicle-earth  horizon  -  48.011° 

Altair-vehicle-earth  horizon  -  56.503° 

Earth’s  subtended  angle  -  14.018° 

The  s Car-vehicle-earth  center  angles  are  then 

Acrux-vehicle-earth  -  50.979° 

Arcturus-vehicle-earth  -  55.020° 

Altair-vehicle-earth  -  63.512° 


Using  data  from  Table  VII,  the  components  of  the  vehicle’s  unit 
position  vector  are 

r  =  (.22619)  cos  50.979°  +  (.84496)  cos  55.020° 

X 

-  (.51440)  cos  63.512° 


l 
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ry  =  (.28928)  cos  50.979°  +  (.35811)  cos  55.020° 

+  (.90990)  cos  63.512° 

rz  =  (.99216)  cos  50.979°  -  (.44293)  cos  55.020° 

+  (.21120)  cos  63.512° 

or 

rv  =  (.22619) (.62960)  +  (.84496) (.57329)  -  (.51440) (.44600) 

Ty  =  (.28928) (.62960)  +  (. 35811) (.57329)  +  (.90990) (.44600) 

rz  =  (.99216) (.62960)  -  (.44923) (.57329)  +  (.21120) (.44600) 

and  carrying  out  the  operations 

rx  =  .39739 

vy  =  . 79324 

rz  =  .46132 

The  radial  distance  from  the  earth  is 


R@  _  3443.93 

14.018°  ~  0.12201 

sxn  - o - 


28,226  NM 


so  the  coordinates  of  the  vehicle  are  (geocentric  equatorial  system) 

x  =  (.39739) (28226)  =  11,216  NM 


y  =  (.79324) (28226)  =  22,389  NM 
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The  declination  of  the  vehicle  is  simply 

5v  =  (Yns  +  A/2)  -  90°  =  117.467°  -  90.000°  -  27.467° 
The  expression  for  Act  is 


cos  Act 


-  cos  (vs  +  A/2)  -  sin  6S  sin  <$v 
cos  6  s  cos  6V 


and  substituting  values 


cos  Act  =  (-.445937)  -  (.152398)  (.46 1232') 

(.988273) (.887279) 

Aa  =  ±  126.08° 

The  astronauts  determine  the  proper  sign  as  "plus"  by  looking  out  the 

A  A 

window  and  establishing  the  direction  of  I...  .  x  I„  ,  as 

Altair  Polaris 

generally  toward  the  earth.  The  right  ascension  of  the  vehicle  is 
then 


oig  *f  Act 

The  right  ascension  of  Altair  is  obtained  from  the  ephemiris  as 
297.305°.  Then 


av  =  297.305°  +  126.08°  =  63.38° 

The  polar  coordinates  of  the  vehicle  are  thus 

otv  =  63.38° 

6V  =  27.47° 
r  =  28,226  NM 
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The  rectangular  coordinates  of  the  unit  position  vector  are  given  by 

rx  =  cos  6V  cos  av 

ry  =  cos  <SV  sin  av  (5-18) 

rz  =  sin  6V 

Substituting  values 

rx  =  (.887279) (.448019)  =  0.39752 

ry  =  (.887279) (.894024)  =  0.79325 

rz  =  (.461232)  =  0.46123 

The  coordinates  of  the  vehicle  are 

x  =  rrx  =  (28226) (.39752)  =  11,220  NM 

y  =  rry  =  (28226) (. 79325)  =  22,390  NM 

z  =  rrz  =  (28226) (.46123)  =  13,019  NM 

which  agree  with  the  values  obtained  using  three-star  fix  proce¬ 
dures.  Note,  however,  that  Polaris  was  assumed  at  the  pole  position. 

Position  Fix  Simultaneous  Measurement  Requirement ♦  The 
position  fix  techniques  presented  require,  Lit  the  ideal  case, 
simultaneous  measurements.  Obviously,  if  only  one  sextant  is 
available,  it  is  impossible  to  measure  all  angles  simultaneously. 

A  number  of  techniques  resolving  the  difficulty  are  available. 
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Linear  interpolation  between  an  early  and  late  sighting  could 
be  used  to  solve  the  simultaneous  measurement  requirement.  If 
linear  interpolation  is  insufficient,  more  measurements  can  be 
obtained  for  use  with  a  higher  order  interpolation  formula  (Ref 
24:  Ch  10).  This  procedure  is  illustrated  in  the  following 
example: 

r 

Example  Problem  5-3  -  Interpolation  for  Measurement.  The 


angular  measurement  of  Arcturus-vehicle-earth  is  desired  at  time 
12:25.  The  following  measurements  are  obtained: 

Reading 

Number  Time  Y 


12:10 


41.265' 


12:20 


39.271 


12:30 


37.370 


12:40 


35.568 


Bessel's  interpolation  formula  (Ref  24:107)  reduces  to 


=  Y2  +  Y3  Yi  -  Y2  -  V3  + 


12:25 


(5-19) 


_  /6.641  0.192  qq  onoo 

Y12:25 - 2 - 16 - 


For  comparison,  linear  interpolation  for  the  angle  at  12:25  yields 


12:25 


•Y  ■?- *  — 3-  =  38.320' 
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A  number  of  other  interpolation  formulas  are  available — see  any  text  on 
numerical  analysis. 

The  extrapolation  of  measurements  to  an  earlier  or  later  time  can 
also  be  achieved  using  an  interpolation  formula.  Another  interesting 
related  technique,  which  has  been  used  by  astronomers  for  hundreds  of 
years,  is  also  available.  The  procedure  will  be  illustrated  by  an 
example. 

Example  Problem  5-4  -  Early  and  Late  Measurement.  The  angular 
measurements  of  Arcturis-vehicle-earth  are  desired  at  12:00  and  12:50. 
The  measurements  of  Example  Problem  5-3  are  to  be  used. 

To  obtain  the  early  and  late  measurement,  the  set  of  measure¬ 
ments  is  used  to  construct  a  difference  table  (Table  VIII).  First, 
the  difference  column  of  the  measurements  is  obtained  by  subtracting 
from  each  reading  the  reading  immediately  preceding  it.  Then  the 
differences  of  the  differences  are  tabulated.  The  table  is  continued 
as  far  as  possible.  Table  VIII  is  the  completed  difference  array 
for  the  set  of  measurements  used  in  this  example. 

To  extend  the  table,  the  last  difference  is  assumed  constant 
and  the  last  column  of  the  difference  array  is  extended.  Starting 
with  the  extended  values  in  the  last  column,  the  table  is  built 
"backward"  to  finally  arrive  at  the  desired  early  and  late  values. 

Table  IX  illustrates  the  procedure.  The  underlined  values  in 
Table  IX  represent  the  extension  of  the  difference  array. 

The  desired  values  at  12:00  and  12:50  are  then: 

12:00  -  43.346° 

12:50  -  33.871° 


102 


GA/AE/69-1 


Table  VIII 
Difference  Array 


1st 

2nd 

3rd 

Time 

Angle 

Diff. 

Diff. 

Diff. 

12:10 

41.265° 

-1.994 

12:20 

39.271° 

-1.901 

.093 

.006 

12:30 

37.370° 

1 

j 

-1.802 

.099 

12:40 

35.568° 

Table  IX 

Extended  Difference  Array 


Time 


12:00 


12:10 

12:20 

12:30 

12:40 

12:50 


Angle 


43.346 


41.265 

39.271 

37.370 

35.568 

33.871 


3rd 

Diff. 


-2.081 


-1.994 


-1.901 


-1.802 


-1.697 
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The  differencing  technique  is  more  fully  treated  in: 

1.  Ref  2:135 

2.  Ref  24:98 

3.  Ref  6:134 
and  many  others. 

The  simultaneous  measurement  problem  is  also  considered  in 
Ref  8:18.  The  technique  employed  is  to  convert  measurements  to  a 
common  time  by  the  use  of  pre-calculated  data  on  the  rate  of  change 
of  angle  with  time  along  a  reference  trajectory.  The  method  is  fully 
treated  in  Ref  8  and  thus  will  not  be  repeated  here. 

Another  possibility  is  the  use  of  a  photograph  to  obtain  the 
angular  sightings — with  a  photograph,  the  measurements  can  indeed 
be  simultaneous.*  In  any  event,  the  simultaneous  measurement 
requirement  represents  nothing  more  than  a  slight  inconvenience  to  the 
astronaut. 

Graphical  Solution  -  Unit  Position  Vector 

The  graphical  solution  for  the  unit  position  vector  is  straight¬ 
forward  when  the  temporary  view  is  taken  that  the  earth  is  orbiting 
the  spacecraft.  The  graphical  solution  solves  the  equations 

A  A 

r  •  I1  =  -  cos  (y i  +  A/2) 

A  A 

r  •  I2  =  -  cos  (y2  +  A/2)  (5-21) 

|  r  1  *  1 

*  Additionally,  photographs  of  the  stars  are  in  gnomonic  projection 
(Ref  15:81),  thus  perhaps  aiding  in  the  graphical  solution  for 
Cty  and  6V  discussed  later. 
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and  yields  the  right  ascension  and  declination  of  the  vehicle.  A  third 
precision  measurement  is  unnecessary  and  the  only  calculation  involved 
is  addition  or  subtraction.  The  graphical  solution  is  illustrated  in 
fig.  27.  The  step-by-step  procedure  is  as  follows: 

1.  Determine  the  star-vehicle-earth  center  angle 
y  +  A/2  for  two  stars. 

2.  Find  the  same  two  stars  on  a  star  chart.  The  chart 
projection  must  be  such  that  a  straight  line  represents 
a  great  circle  (gnomonic  projection). 

3.  Using  Star  1  as  a  center,  lay  off  an  arc  of  radius 
Yj  +  A/2.  The  angular  scale  on  the  equator  or  a 
meridian  must  be  used  to  set  the  radius. 

4.  Using  Star  2  as  a  center,  lay  off  an  arc  of  radius 
Y2  +  A/2. 

5.  The  arcs  intersect  at  two  points.  Select  the  proper 
point  by  an  eyeball  sighting  of  a  third  star — an  arc 
drawn  using  the  third  star  will  fall  close  to  the 
proper  point..  (Or  observe  Polaris  and  roughly  determine 
the  earth’s  declination  using  60  =  -6V.) 

6.  The  plotted  point  represents  the  apparent  position  of 
the  earth  (the  temporary  view  has  been  taken  that  the 

t  earth  is  orbiting  the  spacecraft) . 

I 

7.  Determine  the  vehicle  coordinates  using 

°v  =  a®  ±  180° 
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Fig.  27 

Graphical  Solution  For  Unit  Position  Vector 


Example  Problem  5-5  -  Graphical  Solution  For  ov  and  6V.  The 
situation  described  in  Example  Problem  5-1  will  be  used  to  illustrate 
the  graphical  determination  of  vehicle  right  ascension  and  declination. 
The  angular  scale  of  the  3tar  chart  (Fig.  28)  has  been  constructed  to 
fit  the  thesis  page — accuracy  is  therefore  sacrificed. 

At  TLI  plus  three  hours,  the  astronauts  determine  the  angles 
Arcturus-vehicle-earth  center  =  55.020° 

Rasalhague-vehicle-earth  center  =  44.491° 
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The  graphical  solution  is  completed  on  Fig.  28.  Using  the  star 
Arcturus  as  center,  an  arc  of  radius  55.0°  is  drawn.  Another  arc 
of  radius  44.5°  is  drawn  using  P.asalhague  as  center  point.  By  an 
eyeball  sighting  of  Polaris,  the  vehicle  declination  is  around 

6V  =  110°  -  90°  =  +20° 

so  the  declination  of  the  earth  as  viewed  from  the  spacecraft  is 
around  -20°.  With  this  knowledge,  the  proper  arc  intersection  point 
is  read  as 

a®  =  244° 

5®  -  -27 0 

The  coordinates  of  the  vehicle  are  then 

av  =  244°  -  180°  -  64° 

6V  -  -  (-27°)  =  +27° 

Limitations .  The  surprising  simplicity  of  the  graphical  solution 
is  marred  by  the  fact  that  it  yields  only  an  approximate  angular 
position.  The  angular  scale  of  the  gnomonic  projection  varies  over 
the  chart,  therefore,  the  intersecting  arcs  are  not  truly  arcs  of 
circles.  However,  within  a  limited  range  of  the  chart  (the  central 
region  of  Fig.  28),  the  scale  is  fairly  uniform;  positions  obtained 
within  this  region  will  be  close  to  the  actual  spacecraft  coordinates. 

Since  the  chart  can  only  yield  an  approximate  position,  the 
graphical  solution  is  limited  to  being  a  cross-check  of  more 
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28.  Example  of  Graphical  Position  Fix 


GA/AE/69-1 


sophisticated  fixing  means.  For  example,  after  obtaining  av  and  <SV  by 
other  means,  the  astronauts  can  quickly  check  their  calculations  by 
determining  the  vehicle  angular  position  graphically.  Gross  errors  in 
computation  can  thus  be  detected  without  actually  repeating  the  cal¬ 
culations  . 

Velocity  Determination 

An  intermediate  step  before  some  orbit  computation  methods  is 
determination  of  the  vehicle's  velocity  vector.  In  general,  this  is 
an  iterative  computer  problem.  Several  approximate  methods,  however, 
are  available. 

Numerical  Differentiation .  Any  of  the  numerical  differentiation 
formulas  introduced  in  Chapter  II  can  be  used  to  obtain  velocity. 

For  three  position  fixes  separated  by  time  t 


r2 


(5-22) 


and  for  four  fixes 


3(2ra  -  r,)  -  (2?,  +  r„) 

6t 


(5-23) 


Each  vector  equation  represents  three  scalar  equations.  Thus,  Eq 
(5-22)  represents  three  equations,  the  first  of  which  is 


x3  -  x 

X,  =  ~ - - 

2  2t 


and  similar  equations  for  y^  and  z2< 


(5-24) 
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Lambert's  Theorem.  Two  vector  positions  and  their  time  separation 
are  sufficient  for  velocity  determination  using  Lambert's  theorem 
(Ref  4:80).  Solution  of  Lambert's  theorem  involves  a  double  iteration 
process  and  is  thus  impractical  unless  a  back-up  computer  is  available. 

Orbit  Determination  From  Position  and  Velocity 

Knowledge  of  position  and  velocity  at  a  specific  time  completely 
defines  an  orbit.  Methods  for  computing  the  orbital  elements  from 
position  and  velocity  are  presented  in  most  texts  on  astrodynamics . 

One  of  the  orbit  computation  methods  discussed  in  Ref  3:Ch  1  is 
briefly  outlined  step-by-step  in  Table  X.  The  items  in  Table  X 
labeled  "intermediate  steps"  are  fully  derived  and  discussed  in 
Ref  3:Ch  1. 

Figure  29  shows  the  orientation  unit  vectors  P,  Q,  and  W  used  in 
the  procedure  as  well  as  a  pictorial  definition  of  the  descriptive 
orientation  parameters  i,  ft,  and  w.  Two  of  the  unit  vectors  establish 
the  orbit  orientation — the  th?'‘d  unit  vector  as  well  as  i,  ft,  and  to 
are  not  really  required. 

t 

The  computation  method  of  Table  X  can  be  extended  to  include 
parabolic  and  hyperbolic  trajectories  (Ref  3:Ch  1). 

Orbit  Determination  From  Three  Position  Fixes:  Gibbsian  Method 

An  approximate  orbit  can  always  be  obtained  from  three  fixes  by 
using  numerical  differentiation  to  determine  velocity  at  one  of  the 
positions.  But  such  methods  are  only  approximate.  The  Gibbsian 
method,  however,  yields  the  orbit  with  no  need  for  approximations  or 
iterations  of  any  kind.  The  computed  orbit  using  the  Gibbsian  method 
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Table  X 

Orbital  Elements  From  Position  and  Velocity 


.  [il 


r  =  position  vector 
r  =  velocity  vector 


Compute  |r|  and  S2 
r2  =  (x2  +  y2  +  z2) 
S2  =  x2  +  y2  +  z2 


S2  =  I #1 2 


Compute  a 


(2u/r)  -  S2 


Yields  semimajor  axis  a 


Compute  D  d  =  intermediate  step 

D  =  rr/Zp"  =  (xx  +  yy  +  zz)  .  ± 

vy  r  =  radial  velocity  f  |r| 


Compute  D 


D  =  Zu 


’  1  _  if 

r  a 


D  =  intermediate  step 


T,  LU 

e  cos  Efl  =  —  =  1  -  • 


Yields  e  and  E, 


ax  =  (Dx  -  xD)//5T 
ay  =  (Dy  -  yD) / Zp" 
az  =  (Dz  -  zD)/ZF 


SL^>y>z  =  intermediate  steps 

(ax,y,z  does  not  refer  to 
semimajor  axis  a) 


3 


<  s 


GA/AE/69-1 


Fig.  29 


AAA 

Orientation  Parameters  and  Vectors  P,  Q,  W 


is  as  accurate  as  the  position  fix  input  data. 

The  Gibbsian  method  is  treated  in  Ref  3:51.  A  step-by-step 
outline  of  the  computation  procedure  is  presented  in  Table  XI. 
Note  that  the  time  interval  between  fixes  is  not  required.  The 
computation  procedure  can  be  extended  to  include  hyperbolic 
trajectories  and  variations  exist  that  use  redundant  data  (Ref 
3:Ch  1). 
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Table  XI 


Gibbsian  Method:  Orbit  Determination  From  Three  Position  Fixes 


<  is 
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GA/AE/69-1 


measured  with  the  sextant.  The  angle  6,  which  corresponds  to  the  true 
anomaly  of  the  orbit  plane  star,  is  obtained  from 

I  •  P  =  -  cos  3  (5-31) 

If  only  the  difference  in  true  anomaly  between  two  positions  is  needed, 
then 


02-6!  -  (v  +  A/Z>2  -  (y  +  A/2)! 

and  there  is  no  need  to  determine  0. 

A  star  in  the  orbit  plane  can  be  found  by  establishing  a  and  6 
of  two  known  orbital  points,  plotting  these  points  on  a  properly 
constructed  star  chart,  and  joining  them  with  a  straight  line.  Any 
star  on  the  line  is  "in"  the  orbit  plane.  The  two  known  points  may 

A  A 

consist  of  two  unit  position  voctois  or  the  unit  vectors  P  and  Q. 
Knowledge  of  1  and  ft  also  establishes  the  line. 

The  type  of  star  chart  is  important  since  the  straight  line 
drawn  between  the  two  known  points  must  represent  a  great  circle  on 
the  celestial  sphere.  Projections  having  this  characteristic  are: 

1.  equatorial  gnomonic 

2.  oblique  gnomonic 

3.  polar  gnomonic 

True  Anomaly  from  any  Star.  If  a  star  in  the  orbit  plane  cannot 
be  found,  the  situation  is  more  complicated.  However,  it  is  possible 
to  sight  any  star  and  derive  a  measure  of  true  anomaly.  Preliminary 
data  needed  are  two  angular  positions — right  ascensions  and 
declinations  or  the  unit  position  vectors.  Complete  derivations  of 
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the  equations  to  follow  may  be  found  in  Appendix  D. 

An  illustration  of  the  orbital  trace  across  the  celestial  sphere 


is  shown  in  Fig.  31. 
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where 

Ys  +  (A/2)  =  the  measured  star-vehicle-earth  center  angle. 

ij)0  =  the  minimum  possible  star-earth-vehicle  angle.* 

0  =  true  anomaly  of  vehicle. 

0S  =  true  anomaly  of  star's  unit  vector  as  projected 
into  orbit  plane. 

A  =  earth's  subtended  angle. 

Note  that  ij>0  represents  the  minimum  possible  star-earth-vehicle  angle. 
Therefore,  if  the  star  and  the  earth  can  be  observed  simultaneously, 
the  vehicle  must  travel  close  to  180°  before  \pQ  will  actually  occur. 
The  angle  t|»Q  can  be  determined  from 

-  cos  (ys  +  A/2) „ 

cos  ip0  = - 0  <  <  90°  (5-34] 

cos  (0S  -  02)  “  0  ~ 

where  (ys  +  A/2)2  represents  the  star-vehicle-earth  center  angle  as 
measured  at  position  two.  The  angle  (0S  -  02)  can  be  determined 
from 

tan  (0g  -  02)  =  cot  (02  -  9i) 

cos  (ys  +  A/2) x 

cos  (ys  +  A/2)2  sin  (02  -  0^  (5-35] 
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and  (02  -  0j)  from 


The  angle  tf>0  for  a  given  star  is  constant.  Therefore,  once 

is  established,  the  true  anomaly  difference  between  two  positions 
o 

can  be  determined  through 


cos  (0-  -  6)  = 


-  cos  (ys  +  A/2) 
cos  t|>0 


(5-37) 


For  the  special  case  when  the  star  is  in  the  orbit  plane  cos  \pQ  =  1 
and  Eq  (5-37)  reduces  to  the  simpler  expression  treated  in  the  first 
part  of  this  section. 

The  angle  ipQ  can  also  be  determined  by  direct  measurement 
(Ref  13:8).  The  astronaut  simply  observes  a  star  passing  close  by 
the  earth  and  records  the  minimum  star-vehicle-earth  angle.  The 
minimum  star-vehicle-earth  angle  corresponds  to  i|»Q.  Figure  32  is 
an  illustration  of  this  concept. 


Conclusion 

The  manual  determination  of  the  three-dimensional  fix  is  within 
the  realm  of  possibility.  With  three  fixes,  the  orbit  can  be 
obtained  using  straightforward  computations.  Although  direct  manual 
computation  of  the  orbit  will  be  tedious  and  time  consuming,  the 
capability  is  there.  Simplifying  techniques  and  procedures  will 
become  apparent  if  the  method  is  adopted  and  analyzed  step-by-step. 
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Star-Vehicle-Earth  (SVE)  Plane  at  Instant 
of  ip0  Measurement  (Orbit  Plane 
and  SVE  Plane  Perpendicular) 

If  a  back-up  computer  is  available,  several  avenues  are  open: 

1.  Program  the  approximate  orbit — i.e.,  obtain  velocity 
using  numerical  differentiation. 

2.  Solve  Lambert's  theorem. 

3.  Program  the  Gibbsian  method. 

Programming  an  approximate  orbit  is  undesirable  if  better  methods  are 
available.  Lambert's  theorem  involves  a  double  iteration  which  would 
probably  exceed  a  back-up  computer's  capability.  Obviously,  the 
Gibbsian  method,  or  a  variation,  should  be  programmed.  Chapter  IX  will 
present  such  a  program.  The  output  of  this  program  will  be  used  as 
input  to  a  second  program  that  "projects"  position  and  velocity  as  a 
function  of  time  along  the  trajectory.  This  information  could  be  used, 
in  turn,  as  input  to  a  back-up  guidance  scheme. 
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VI .  Orientation  Parameters  -  Geometric 
Elements  Known 

Summary  of  Present  Manual  System 

References  12,  17,  and  22  present  a  manual  means  for  determining 
the  orientation  parameters  i  (inclination) ,  ft  (longitude  of  ascending 
node),  and  w  (argument  of  perifocus).  The  method  consists  of  determin¬ 
ing  minimum  coaltitudes  for  a  selected  star  pair  (minimum  coaltitude 
is  the  minimum  star-earth-vehicle  angle  for  a  given  star) .  Given  two 
minimum  coaltitudes,  and  (ip 0)2,  the  following  nonlinear  set  of 

A 

equations*  is  solved  for  the  components  of  W 

W  •  Ij  =  sin  (tl^h 

W  •  I2  =  sin  (4>0)2  (6-1) 

|W|  =  1 

A  A 

where  the  unit  vector  W  is  normal  to  the  orbit  plane  (Fig.  29). 

a 

and  I2  are  unit  vectors  toward  star  one  and  star  two  respectively. 

A 

The  orientation  elements  are  then  determined  from  knowledge  of  W 
and  presented  graphically  as  a  function  of  (^0)j  and  (ij/0)2. 
system  is  ingenious  and  works,  but  there  are  several  disadvantages 
in  determining  the  orientation  parameters  using  this  method. 


*  Notice  the  similarity  of  this  set  of  equations  to  the  nonlinear 
three-dimensional  fix  (Eq  (5-12) )  . 
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1.  The  geometric  parameters  are  used  in  determining 

and  ')0)2'  Therefore,  any  errors  in  the 
geometric  parameters  will  be  carried  over  to  the 
computation  of  each  orientation  element. 

2.  The  system  is  dependent  upon  pre-computations  using 
pre-selected  star  pairs — star  pairs  are  chosen  with 
regard  to  the  expected  orbit  and  an  established  range 
of  probable  values  of  i,  Q,  and  w.  Therefore,  the 
method  is  not  entirely  flexible. 

3.  The  solution  of  the  nonlinear  set  of  equations, 

Eq  (6-1) ,  is  a  disagreeable  task  even  for  a  ground 
based  computer.  Programming  a  small  back-up  computer 
for  inflight  solution  of  Eq  (6-1)  is  of  questionable 
practicality. 

Alternate  methods  for  determining  the  orientation  parameters  are 
available.  The  unit  vector  W  can  be  obtained  from  any  two  vectors 
known  to  be  in  the  orbit  plane.  For  example,  if  Tj  and  ?2  represent 
the  vehicle  position  vector  at  time  one  and  later  time  two,  then 


A 

w 


ri  X  r2 
ri  X  *2 


(6-2) 


The  position  vectors  r1  and  ?2  can  be  determined  by  the  methods 
discussed  in  Chapter  V. 

The  object  in  this  discussion  is  not  blind  criticism  of  a  manual 
method  which  works  (and  quite  well  -  see  Ref  12) ,  but  is  rather  an 


GA/AE/69-1 


attempt  to  start  development  of  a  general  method  of  orientation  element 
computation  applicable  to  any  orbit  whatsoever — with  no  requirement 
for  prior  knowledge  of  what  that  orbit  might  be.  The  following 
paragraphs  treat  the  problem  of  orientation  element  computation  from 
partial  knowledge  of  two  positions  in  space. 


Orientation  Elements  From  Position  Fixes 

If  the  geometric  elements  are  obtained  by  other  means,  only  two 
partial  fixes  are  needed  to  specify  the  orientation  elements  i,  S2, 
and  a)  (defined  in  Fig.  29).  In  the  development  to  follow,  it  is 
assumed  that  two  such  fixes  have  been  obtained  using  one  of  the 
methods  of  Chapter  V. 

At  position  one  and  a  later  position  two,  the  unit  position 
vectors  are 


A 


Xl 


i  + 


ryi  J 


+  r 


zl 


r  =  r  i+r  j+r  k 
2  x2  y2  z2 


(6-3) 


A 

Then  the  unit  vector  W  (see  Fig.  29)  is  given  by 


from  which 


(6-4) 


cos  i  =  Wz  (0  i  <_  ir) 


(6-5) 
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and 


tan  ft 


The  ambiguity  for  ft  is  resolved  by  noting  that 


W 

sin  ft  = 

sin  1 


(6-6) 


(6-7) 


Thus,  the  orbital  inclination  i  and  the  longitude  of  ascending 
node  ft  are  determined  with  no  requirement  for  prior  knowledge  of 
the  geometric  parameters. 

The  solution  for  w,  the  argument  of  perifocus,  requires  knowledge 
of  the  true  anomaly  at  some  position  r.  This  can  be  position  one, 
two,  or  any  other  known  position.  Therefore,  the  geometric  elements 
must  be  determined  before  w  can  be  calculated.  The  argument  of 
perifocus  is  given  by 


where  0  is  the  true  anomaly  at  the  vehicle  position  r.  The  ambiguity 
for  w  can  be  resolved  using 


co 


cos 


-1 


x 


cos  ft 


r2  tan 
tan  i 


(6-9) 


Example  Problem  6-1  -  Orientation  Elements .  Several  spacecraft 
emergencies  have  compounded  to  the  point  that  manual  determination 
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of  the  orientation  elements  is  necessary.  Using  the  three-star  fix 
procedures  of  Chapter  V,  the  astronauts  obtain  the  following  unit 
position  vectors  and  corresponding  radial  distances: 


t  = 


t  =  tj  +  3  1/2  hr 


rxl  =  • 50837 
ryl  «  .74113 

rzl  =  .43850 

r  =  21098  NM 


rx2  *  ‘ 25442 

ry2  =  .83992 

rz2  =  .47937 

r  =  unknown — but 

‘  greater  than 


The  cross-product  rx 

A  A 

r2Xr2  = 


x  r2  is 


i  (r 


yi 


given  by 


z2 


-  r 


zl 


A 

+  j 


Zl 


r 

xl 


z2 


) 

) 


+  k  (r  ro-r  r  ) 
xl  y2  yl  x2 


(6-10) 


and  substituting  values 


i1  x  r2  =  i  (-.01303)  +  j  (-.13213)  +  k  (.23843) 


x  r0 


=  0.27290 


A 

The  unit  vector  W  is  then 


W  =  i  (-.04774)  +  j  (-.48417)  +  k  (.87369) 
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from  which 


cos  i  =  0.87369 


and 


i  =  29.11° 


The  longitude  of  the  ascending  node,  ft,  is  given  by 


tan  ft  «  _ —  =  (.04774)  _ 


W„ 


(-.48417) 


=  -0.09860 


from  which 


ft  = 


174.36° 

or 

354.36° 


Since  sin  ft  -  Wx/sin  i  is  negative,  the  correct  value  is 


ft  =  354.36°  =  -5.64° 


The  true  anomaly  at  one  of  the  positions  must  be  determined  before 
the  argument  of  perifocus,  u>,  can  be  calculated.  The  geometric 
elements  (assumed  previously  determined)  yield  the  true  anomaly 
position  one  through  the  equation 


at 


cos  0 


-  a  (1  -  e2)  -  r, 
1  - 1 


er. 


(6-11) 
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Using  the  values 


a  =  132,875  NM 


e  =  .9737 


r  =  21,098  NM 


the  expression  for  cos  6j  is 


cos  0, 


-  (132875) (1  -  .9481)  -  21098 

(.9737) (21098) 


=  -  0.69132 


01  =  133.73° 


The  true  anomaly  is  less  than  180°  since  the  radial  distance  is 
increasing.  The  argument  of  perifocus  is  given  by 

.  _i  rzl  .  4-1  .43850  .-o  -oo 

(d  =  sin  1  — — -  -6.  =  sin  1  —,-zr<rr  ~  133.73° 

sm  i  1  .48651 


64.33° 

or 

115.67° 


-  69.40 


-  133.73°  = 


-  18.06 


The  proper  angle  is  determined  by  a  rough  calculation  or 


rvi  r_.  tan  SI 

a)  =  cos  1  — +  _JLi - - —  -  e 

cos  S2  tan  x 
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which  becomes 


0) 


cos 


-1 


(.51)  _  (.44) (.1) 
(.99)  (.56) 


-  134° 


Then  roughly 


ov 

o 

_ l 

0 

O 

0)  - 

or 

-  134°  = 

or 

-64° 

-198° 

The  angle  common  to  both  calculations  is  -69.40°  so  the  argument  of 
perifocus  is 


a  =  -  69.40° 

Orientation  Elements  from  av  and  6y 

Knowledge  of  vehicle  right  ascension  and  declination  at  two 
points  coupled  with  the  true  anomaly  at  either  point  is  sufficient  to 
define  the  orientation  elements.  Conversion  from  polar  to  rectangular 
coordinates  is  unnecessary.  Figure  33,  an  illustration  of  the 
vehicle's  path  on  the  celestial  sphere,  will  be  used  for  some  of  the 
derivations-  to  follow. 

From  the  previous  section 

cos  i  =  Wz  =  (rxl  ry2  -  ryl  rx2)/ 


rl  x  r2 


(6-12) 


t 

I 


I 
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The  argument  of  perifocus  requires  knowledge  of  the  true  anomaly 


at  one  of  the  positions.  Assume  ©2  is  known.  Then,  from  the 
spherical  triangle  acd 


sin  6, 

sin  (<o  +  8,)  =  -  (6-19) 

*  sin  i 


The  ambiguity  can  be  resolved  using 


sin 

sin  (w  +  0.)  =  — — (6-20) 
1  sin  i 


Example  Calculation  6-2  (i,  0,  and  io  from  a_  and  6) .  Using  the 
Polaris  fix  technique  (Chapter  V)  the  astronauts  establish  the 
following  angular  data: 

t  =  tx  t2  =  tj  +  2  hr 

<*!  =  55.55°  ot2  =  68.29° 


6j  =  26.01° 


S2  =  28.15° 


In  addition,  the  true  anomaly  at  position  one  is  known  to  be 

0j  =  133.73° 

The  orientation  of  the  orbit  is  desired. 
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Solution.  The  true  anomaly  difference  between  position  one  and 


two  is 


cos  A0  =  sin  61  sin  <$2  +  cos  cos  62  cos  (a2  -  (6-15) 


and  substituting  values 


cos  A 9  -  (.43851) (.47178)  +  (.89873) (.88171) (.97541) 

=  .97981 


from  which 


A8  =  11.53e 


sin  A0  =  .19991 


The  inclination  is  given  by 


cos  ^  cos  62  sin  (a2  ~ 
sin  A0 


(6-16) 


and  substituting  values 


cos  i  =  (.89873) (.88171) (.22042)  =  >g7372 

(.19991) 


i  =  29.11° 
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The  longitude  of  the  ascending  node  is  found  from 


sin  (a2  -  ft) 


tan  6 2  (.53508)  = 

tan  i  '  (.55673) 


.96110 


“2  -  ft 


73.97° 

or 

106.03° 


■  68.29°  -  ft 


-5.68° 

0  =  or 

-37.74° 


The  ambiguity  is  resolved  by  a  rough,  calculation  of  ft  using  the 
other  position.  Accordingly 

tan  6 ,  aq 

sin  (a.  -  ft)  =  - - -i.  -  =  .87 

1  tan  i  .56 


-5° 

ft  =  or 
-65° 


The  angle  common  to  both  calculations  is  near  -5°  so 


ft  =  -5.68° 


The  argument  of  perifocus  is  determined  from 


sin  (to  +  0j) 


<j)  +  0, 


sin  6} 
sin  i 


64.35° 

or 

115.65° 


(.43851) 

(.48643) 


=  .90149 


=  a)  +  133.73° 


i 
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The  common  value  is  selected,  i.e., 

o)  =  -69.38° 

Graphical  Solution  for  jL ,  ,  and  w. 

A  graphical  solution  for  the  orientation  parameters  is  illus¬ 
trated  in  Fig.  34.  The  graphical  solution  for  i,  J2,  and  oi  consists 
of  the  following  steps: 

1.  Plot  two  known  unit  position  vectors  on  a  properly 
constructed  star  chart,  i.e.,  plot  (avl,  <5vl) 

and  (ov2>  <5v2)  where  position  two  is  later  than 
position  one. 

2.  Draw  a  straight  line  from  position  two  through 
position  one  and  inte_sect  the  equator. 
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Celestial  Sphere-Equatorial  Gnomonic  Projection  j 


Fig.  34 

Graphical  Solution  -  Orientation  Elements 


3.  The  right  ascension  of  the  intersection  point  is  ft.* 
The  angle  between  the  orbit  plane  and  the  equatorial 
plane  is  the  inclination — but  the  angle  i  cannot  be 
measured  directly  from  the  chart.  Instead,  the 
inclination  must  be  calculated  using  the  right 


The  intersection  will  yield  the  descending  node  if  the  line  is 
drawn  from  position  one  to  two  to  the  equator. 
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spherical  trigonometry  relationship 


tan  i 


tan  (5V 
sin  (av  -  fl) 


(6-22) 


4.  Assume  the  true  anomaly  of  position  one  is  known. 

Plot  the  point  representing  perifocus  by  measuring 
backward  on  the  constructed  line  (i.e.,  in  direction 
of  point  two  to  point  one)  from  point  one  an  angular 
distance  equal  to  the  true  anomaly  (the  angular 
scale  on  the  equator  or  a  meridian  must  be  used) . 

The  angular  distance  from  the  equator  to  perifocus, 
measured  on  the  constructed  line,  is  to. 

The  star  chart  must  be  constructed  such  that  a  straight  line  repre¬ 
sents  q  portion  of  a  great  circle.  Projections  having  this 
characteristic  are: 

1.  equatorial  gnomonic 

2.  oblique  gnomonic 

3.  polar  gnomonic 

The  scale  distortion  of  the  gnomonic  projection  does  not 
appreciably  affect  the  accuracy  of  the  graphical  determination  of  £2 
and  i.  However,  the  distorted  scale  must  be  used  to  determine 
0) — therefore,  w  can  only  be  determined  approximately. 

Example  6-3  (Graphical  Solution  for  i  and  £2) .  Using  the 
Polaris  fix  techniques  (Chapter  V),  the  astronauts  establish  the 
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following  angular  data: 


c 

t2  =  tj^  +  2  hr 

o  =  39.12° 

a  =  63.39° 

2 

fi  =  21.40° 

6  =  27.47° 

2 

The  inclination  of  the  orbital  plane  and  the  longitude  of  the 
ascending  node  is  desired. 

Solution.  Figure  35  represents  an  equatorial  gnomonic  projection 
of  the  celestial  sphere.  The  two  angular  positions  are  plotted  and 
a  straight  line  is  drawn  from  position  two  through  position  one.  The 
line  intersects  the  equator  at  a  =  -5.5°,  therefore 

fl  =  -5.5° 


The  inclination  is  given  by 


tan  i 


tan  6  2 
sin  (o2  -  0) 


tan  27.47° 
sin  (63.39°  +  5.5°) 


tan  i 


.51983 

.93285 


.55725 


from  which 


i  =  29.1° 

The  solution  for  u)  is  not  shown  since,  for  the  example  orbit 
selected,  the  construction  lies  outside  the  range  of  the  chart. 
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Orientation  Elements  From  Minimum  Coaltitudes 

This  section  is  an  outline  of  an  alternate  method  for  determining 

A 

the  unit  vector.  W  (Fig.  29).  The  scheme  is  related  to  the  present 
minimum  coaltitude  method  used  in  Refs  12,  17,  and  23.  The  present 
method  was  briefly  summarized  at  the  opening  of  this  chapter. 

A 

Three  minimum  coaltitudes  can  yield  the  unit  vector  W  through 
the  equations 

W  •  £x  =  sin  (i|)0)  j 

A  A 

W  •  I2  =  sin  (t|»0)2 

A  A  ♦ 

W  •  I3  =  sin  (^0)g  (6-23) 

Eq  (6-23)  is  the  same  form  as  the  three-dimensional  fix  equations  of 
Chapter  V,  F.q  (5-13) .  Solution  for  the,  components  of  the  unit  vector 

A 

W  leids  to 

Wx  =  “  Axl  siTl  ^o>i  “  Bx2  sin  ^+0)2  "  cx3  sin  ^0)3 

wy  ■  -  Ayi  sin  (l|)G)1  -  By2  sin  0|>o)2  -  Cy3  sin  0|>o)3 

(*o>,  -  “z2  «o>2  -  Cz3  sin  (*0)3  ;6_24) 

The  constants  A^  through  CZ3  are  the  same  quantities  defined  by 

Eq  (5-16)  and  used  in  the  three  star  fix  procedure  of  Chapter  V. 

The  same  publication  tabulating  the  constants  can  therefore  be  used 
either  for  three  star  fix  procedures  or  for  directly  determining  the 

A 

unit  vector  W.  The  orientation  elements  can  be  derived  from  knowledge 
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of  W. 

The  minimum  coaltitudes  can  be  dete'-mined  geometrically  using 
spherical  trigonometry  if  the  geometric  elements  are  known  (Ref 
17:45).  Minimum  coaltitude  can  also  be  observed  directly — the 

\ 

vgeometry  is  shown  in  Fig.  32  in  Chapter  V. 

Other  than  eliminating  the  simultaneous  measurement  requirement 
of  the  three-dimensional  fix,  the  minimum  coaltitude  method  offers 
little  advantage  over  the  position  fix  technique.  Indeed,  since 
similar  equations  must  be  solved  for  either  method,  the  position  fix 
is  preferred.  A  position  fix  in  space  yields  more  information  the 
astronaut  can  use  directly  in  solving  the  navigation  and  guidance 
problem. 

Summary  -  Orientation  Elements 

Compared  to  the  geom  trie  elements,  the  solution  for  the 
orientation  parameters  is  simple.  The  exact  solution  for  the 
geometric  elements  from  a  set  of  range  measurements  involves  a  set 
of  nonlinear  equations.  The  orientation  elements,  with  the  exception 
of  w  which  depends  on  true  anomaly,  can  be  obtained  using  simpler 
equations.  There  is  no  need  to  toil  over  nonlinear  interative 
techniques  to  obtain  the  orientation  elements.  Notice  that  the 
solutions  for  i  and  fl  do  not  require  knowledge  of  the  vehicle's 
distance  from  the  earth — a  plus  feature  since  the  range  measurement 
is  the  most  difficult  to  obtain  accurately. 
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VII .  Orbit  Determination  -  No  Range  Measurements 

Most  of  the  orbit  determination  methods  presented  in  this  thesis 
require  one  or  more  measurements  of  range  (range  is  radial  distance 
from  the  center  of  the  earth).  In  general,  the  range  measurement  is 
the  most  difficult  to  obtain  accurately.  For  this  reason,  it  seems 
logical  to  seek  an  orbit  determination  scheme  that  does  not  require 
any  range  measurement  whatsoever.  One  such  method  will  be  briefly 
described.  Time  permits  the  presentation  of  only  the  bare  details — 
the  recommendation  section  of  this  thesis  will  suggest  further 
analysis  of  this  orbit  determination  scheme. 

Orbit  Determination  -  N£  Range  Measurements 

One  possible  method  of  computing  an  orbit  without  knowledge  of 
radial  distance  is  briefly  presented. 

1.  Determine  the  orbit  plane — two  angular  positions 
are  required  from  which  i  and  ft  can  be  determined. 

2.  Find  a  star  in  the  orbit  plane.  This  can  be  done 

using  procedures  described  in  Chapter  V. 

•  • 

3.  Obtain  three  values  of  6.  0  might  be  determined 

using  numerical  differentiation.  For  three 
measurements 


02  = 


2t 


(7-1) 
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Note  that  if  a  suitable  measurement  schedule  is 
followed j  a  number  of  terms  can  be  pre-calculated. 

5.  Compute  l  (semilatus  rectum)  from 

— 1 1/3 

2  -  na  +  e  cos  6,)1*  ,  a(1  .  e2)  (7.8) 

L  ^  J 

Summary 

Only  the  bare  essentials  of  a  "no  range"  orbit  determination  method 
have  been  presented.  Other  orbit  determination  schemes  that  involve 
indirect  measurements  of  distance  are  also  available — the  classical 
Laplacian  method  (Ref  3:36)  is  one  example.  Intuitively,  it  seems 
that  an  orbit  determination  scheme  of  this  type  would  be  highly 
advantageous  because: 

1.  To  date,  the  accuracy  of  range  determination  is 
marginal. 

2.  Angular  measurements  can  be  obtained  more  directly 
and  probably  with  greater  accuracy. 

3.  The  semi-diameter  method  of  obtaining  range  is 
limited  to  close  earth  orbits.  Range  in  deep  space 
requires  the  solution  of  a  more  complex  equation 
(Ref  4:222). 

4.  Range  determination  methods  using  a  reference  trajectory 
are  useless  if  the  reference  trajectory  is  invalidated 
(this  topic  is  treated  in  Chapter  VIII) . 
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VIII.  An  Improved  Optical  Method  of 
Range  Measurement 


Introduction 


As  has  been  seen  in  the  previous  chapters,  accurate  range  or 
altitude  measurement  is  essential  in  most  manual  navigation  schemes. 
Optical  measurement  using  the  planet  angular  diameter  is  subject  to 
large  errors  as  the  distance  from  the  measured  body  increases. 
Therefore,  other  optical  ranging  techniques  have  been  investigated. 
Harold  A.  Hamer,  of  the  NASA  Langley  Research  Center,  has  proposed  a 
method  based  upon  linear  perturbations  from  a  pre-computed  reference 
trajectory  (Ref  8).  This  chapter  deals  with  Mr.  Hamer's  approach  to 
accurate  optical  range  measurement. 


The  Method 

The  geometry  of  the  measurement  procedure  is  depicted  in  Fig.  36. 
The  angle  D,  which  can  vary  between  0  degrees  and  180  degrees,  lies  in 
the  instantaneous  earth-moon-vehicle  plane.  If  Fig.  36  represents  the 
vehicle  on  a  reference  trajectory  at  a  given  time,  then  the  incremental 
range  Ar^  (actual  range  minus  reference  range)  at  that  time  is  obtained 
as  follows: 


em 

sin  A 


sin  C  = 


sin  (A  +  B) 

sin  A 


(8-1) 


rev  *  ^em  cos  B  +  rem  sin  B  cot  A 


(8-la) 
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(From  Ref  8; 6) 
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As  shown  in  Eq  (8-la) ,  the  variables  which  define  rgv  at  a  given  time 
are  the  angles  A  and  B;  therefore 


dr. 


ev 


3rev  3rev 
ev  dA  +  — —  dB 


9A 


9B 


(8-2) 


From  Eq  (8-la) 


9r. 


ev 


9A 


-  r  „  sin  B  esc2  A 
em 


(8-3) 


but 


rmv 

sin  B  =  ~ —  sin  A 
rem 


so  that 


9A 


r 

mv 

sin  A 


(8-4) 


(8-5) 


where 


cos  A 


ev 


+  r  2  -  r  2 


mv 


em 


2r  .  r 
ev  mv 


(8-6) 


Further,  from  Eq  (8-la) 

9rev 

--- —  =  r  (cos  B  cot  A  -  sin  B)  (8-7) 

o  B  em 

where  angle  B  is  given  by  Eq  (8-4)  and  angle  A  is  given  by  Eq  (8-6) . 

The  values  of  rev,  rmv,  and  rem  are  known  for  the  reference 


trajectory. 
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The  variables  A  and  B  in  Eq  (8-2)  refer  to  angular  measurements 
which  must  be  made  onboard  the  vehicle.  Angle  B  cannot  be  measured 
from  the  vehicle.  Therefore,  a  related  angle  which  can  be  measured 
must  be  substituted.  As  shown  in  Fig.  36,  the  angle  F,  measured 
between  a  star  and  the  earth’s  center,  can  be  used  in  the  following 
manner.  For  the  case  in  which  the  star  is  in  the  instantaneous 
earth-moon-vehicle  plane,  dD/dB  =  ±  1,  (this  is  true  because  rem 
and  the  line-of-sight  to  the  star  are  fixed  in  direction  for  the 
instant  being  considered),  where  the  sign  is  determined  by  the 
relative  directions  of  the  earth  and  the  star.  Equation  (8-2) , 
therefore,  can  be  written 


3r 


dr, 


ev 


3r 


ev 


3A 


dA  +  k 


ev 


3B 


dD 


(8-8) 


where  dD  represents  the  change  in  the  star-to-earth  angle  and  k  is 
plus  or  minus  one.  (The  equation  for  determining  k  is  derived  in 
Appendix  E.)  If  the  star  is  not  in  the  instantaneous  plane,  as 
shown  in  Fig.  36,  dD  can  be  determined  from  the  right  spherical 
trigonometric  relationship 

cos  F  =  cos  D  cos  E  (8—9) 


or 


D  = 


(8-9 a) 
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Thus 


or 


dD 


sin2  F 
E  -  cos 


2 


1/2 


dF 


(8-10) 


(8-11) 


Therefore,  in  terms  of  two  angles  which  can  be  measured  from  the 
vehicle  (earth  center  to  moon  center  and  star-to-earth  center) , 
Eq  (8-8)  can  be  written 


dr 


3r 


ev 


ev 


3A 


dA  +  k 


1/2 


1  -  cos2  F 

2  E  -  cos2  F; 


dF 


(8-12) 


or,  for  the  region  in  which  a  change  in  rev  is  linear  with  changes 
in  A  and  B 


Ar 


ev 


AA  +  k  ffev 

3A  3B 


(8-13) 


where 


A  A  s  A  M  A 

actual  reference 

(8-14) 

AF  =  F  ,  -  F 

actual  reference 

(8-15) 

The  incremental  range  Arev,  then,  is  determined  by  the  two  measure¬ 
ments  A(actual)  and  F (actual) ,  and  by  the  pre-calculated  values  of 
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A(reference) ,  F(reference) ,  3rgv/3A,  3rgv/3B,  and  k  [(1  -  cob2 
(cos2  E  -  cos2  F)]1/2.  Equations  for  calculating  F,  E,  and  k,  which 
pertain  to  the  reference  trajectory,  are  derived  in  Appendix  E. 

Linearity  Characteristics .  Mr.  Hamer's  investigation  of  the 
variations  of  Arev  with  AA  and  AB  for  a  typical  lunar  trajectory 
showed  that  at  a  given  time  the  partials  3rev/3A  and  3rev/3B  can 
be  considered  to  be  essentially  constant  over  a  wide  region,  as 
indicated  in  Fig.  37.  (Note  that  the  kilometer  is  used  as  the  unit  of 
distance  in  this  chapter.)  The  reference  trajectory  used  was  selected 
from  Ref  7.  In  Fig.  37,  the  curve  for  AB  was  determined  from  the 
relations 

sin  C  =  sin  A  (8-16) 

r 

em 


and 


B  =  180°  -  (A  +  C)  (8-17) 

where  A  is  held  constant  at  the  reference  value  and  is  given  by 
Eq  (8-6).  The  curve  for  AA  was  determined  from  the  relation 

sin  A  =  —  sin  B  (8-18) 

r 

mv 


where 

rmv  =  (rev2  +  rem2  "  2  rem  rev  cos  B>l/2  <8" 1 9> 

and  B  is  held  constant  at  the  reference  value,  and  is  determined  by 
Eqs  (8-16)  and  (8-17).  The  curves  shown  in  Fig.  37  are  for  a 
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AA  or  AB,  Degrees 


Fig.  37 

Linearity  Characteristics  of  the  Hamer  Method  for  a  Reference 
Distance  of  162,000  KM  From  Earth  (Ref  8:10) 
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reference  time  from  injection  (perigee)  of  16.125  hours.  (Here,  rgv 
is  approximately  162,000  kilometers.)  At  this  time,  a  two  hour 
difference  at  injection  from  the  reference  injection  time  would 
represent  a  difference  in  range  from  the  earth  of  about  13,000  kilo¬ 
meters.  The  data  in  Fig.  37  show  that  this  difference  can  be 
considered  to  be  within  the  linear  region.  Several  reference 
trajectories,  therefore,  could  suffice  for  range  determination  for  a 
given  launch  window 

The  linearity  characteristics  of  3rev/9A  and  9rev/9B  illustrated 
in  Fig.  37  are  representative  of  those  over  most  of  the  earth-moon 
distance.  Actually,  the  linear  approximation  for  9rev/9A  improves  as 
the  moon  is  approached,  as  shown  by  the  derivative  of  Eq  (8-5) >  which 
is 


92i 


ev 


•3  A2 


=  r, 


mv 


cos  A 
sin2  A 


(8-20) 


The  angle  A  approaches  90  degrees  as  the  moon  is  approached.  From 
the  derivative  of  Eq  (8-7) 


32rev 

9B2 


r  (cot  A  sin  B  +  cos  B) 
em  ■ 


(8-21) 


it  can  be  determined  that  the  linear  approximation  for  9r  /9B  will 

ev 

also  improve  as  the  moon  is  approached. 

Effect  of  Proximity  to  the  Earth-Moon  Line.  The  values  of 

9r  /3A  and  3r  /9B  which  are  pre-calculated  and  supplied  to  the 
ev  ev 
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astronaut  in  the  form  of  charts  or  tables  are  shown  in  Fig.  38  for 
vehicle  positions  along  the  translunar  portion  of  the  reference 
trajectory.  The  partials  3rmv/3A  and  Sr^/SC  required  to  determine 
range  from  the  moon  are  also  shown.  These  values  will  approach  those 
shown  for  3rev/3A  and  3rev/3B  at  distances  near  the  earth.  It  is  to 
be  noted  that  for  a  position  fix  near  the  moon,  higher  accuracy  is 
obtained  by  measuring  range  from  the  earth  rather  than  from  the  moon, 
inasmuch  as  the  geometry  gives  lowei  values  for  3rev/3A  and  3rev/3B. 
It  will  be  seen  that  the  error  in  determining  the  incremental  range, 
Arev,  is  proportional  to  the  magnitudes  of  the  partials. 

In  Fig.  38,  the  large  values  shown  for  3r£v/3A  art}  3rev/3B  near 
the  earth  occur  because  of  the  close  proximity  to  the.  earth-moon  line 
(angle  A  approaches  180  degrees).  Thus,  range  determination  using 
Eq  (8-13)  would  be  inaccurate  for  about  the  first  five  hours  (or 
75,000  KM)  from  injection.  These  results  would  be  typical  for  any 
method  based  upon  earth-moon  measurements.  The  accuracy  of  angular 
diameter  measurements  is  much  better  for  the  first  five  hours,  as 
shown  in  Fig.  39.  The  dashed  curve  represents  the  error  in  the  range 
for  an  angular  diameter  measurement  cr  error  of  10  arc-seconds. 

For  times  along  the  reference  trajectory  except  for  the  first 
five  hours,  the  error  in  determining  range  by  Eq  (8-13)  is  relatively 
small,  as  shown  by  the  solid  curve  in  Fig.  39.  This  curve  represents 
the  error  in  the  incremental  range  along  the  translunar  portion  of 
the  reference  trajectory  where  the  angular  measurements  A  and  F  are 
considered  to  have  random  uncorrelated  errors,  each  with  a  standard 
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deviation  a  of  10  arc-seconds,  such  that 


where 


•3F 


k 


F 

2 


3r 


ev 


3B 


(8-22) 


(8-22a) 


In  Fig.  39,  the  star  which  is  used  to  measure  F  is  assumed  to  lie  in 
(or  near)  the  instantaneous  earth-moon-vehicle  plane  so  that  the 
value  of  the  term  [(1  -  cos2  ^(cos2  E  -  cos2  F)]1/2  in  Eq  (8-13)  is  one, 
and  3rev/3B  =  ±  3rev/3F. 

Effect  of  the  Star  Selected  for  Measurement.  The  value  of  the 
quantity  [(1  -  cos2  E^cos2  E  -  cos2  F)]1/2  can  vary  from  one  to  infinity, 
depending  on  the  reference  trajectory  and  the  star  selected  for  the 
measurement  of  F.  The  effect  of  star  position  on  the  error  in  AreV 
is  shown  in  Fig.  40.  The  error  in  Arev  is  a  minimum  when  the 
quantity  [(1  -  cos2  F)/(eos2  E  -  cos2  F)]1/2  is  a  minimum.  This  quantity 
has  its  minimum  with  respect  to  F  at  F  =  90  degrees  and  rises  rapidly 
as  the  extreme  values,  F  =  E  and  F  =  180  degrees  minus  E,  are 
approached.  The  extreme  values  of  F  occur  when  the  star  position  is 
in  the  plane  containing  the  vehicle-earth  line  and  the  vertical  to 
the  earth-moon-vehicle  plane.  The  desired  stars,  therefore,  are  those 
that  lie  near  the  earth-moon-vehicle  plane  (E  within  about  plus  or 
minus  30  degrees)  with  projections  in  that  plane  away  from  the 
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vehicle-earth  line. 

Corrections  for  Nonsimultaneous  Measurements .  In  the  equation  for 
range  determination,  Eq  (8-13) ,  both  measurements  are  assumed  to  have 
been  made  simultaneously.  Since  simultaneous  measurements  are  not 
usually  possible,  errors  would  be  introduced  in  the  results  unless  the 
measurements  were  corrected  or  converted  to  some  common  time.  This 
time  could  be  either  before  or  after  the  actual  time  of  the  measurements. 
If  the  time  corresponding  to  one  of  the  measurements  is  made  the  common 
time,  only  one  measurement  would  have  to  be  converted. 

The  measurements  can  be  converted  to  a  common  time  (small  time 
increments  only)  by  use  of  pre-calculated  data  on  the  rate  of  change 
of  angle  with  time  along  the  reference  trajectory.  These  data  would 
be  practically  the  same  for  any  other  trajectory  reasonably  near  the 
reference,  so  that  they  would  apply  for  making  corrections  when  off 
the  reference  trajectory.  Equations  for  calculating  dA/dt  and  dF/dt 
along  a  reference  trajectory  are  given  in  Appendix  B  of  Ref  8. 

In  converting  a  measurement  to  a  common  time,  the  corresponding 
value  of  dA/dt  or  dF/dt  at  the  given  range  multiplied  by  the  time 
increment  would  be  sufficient  to  calculate  the  incremental  angle  to 
be  added  algebraically  to  the  measured  angle.  The  range  would  have 
to  be  known  only  approximately.  For  trajectories  launched  near 
nominal  injection  time,  the  nominal  range  at  the  time  of  the 
measurements  would  be  adequate,  whereas  for  trajectories  not  launched 
near  the  nominal  injection  time,  the  approximate  range  could  be 
determined  according  to  the  time  from  actual  injection  time. 
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Measurements  of  the  Center  of  £  Body 

Throughout  this  discussion,  it  has  been  assumed  that  optical 
measurements  to  the  center  of  the  earth  and  moon  are  possible.  The 
USAF  Hand-held  Space  Sextant  has  been  qualified  for  star-horizon 
measurements  only.  It  is  felt  by  the  authors  of  this  thesis  that 
this  sextant  could  be  adapted  to  measure  star-earth  center  and 
moon-earth  center  angles.  The  incorporation  of  an  attachment  such  as 
a  filter  enscribed  with  concentric  circles  and  centering  cross-hairs 
has  been  suggested  to  the  Air  Force  Avionics  Laboratory,  sponsors  of 
the  sextant's  development  and  testing.  Their  response  has  been 
favorable,  and  this  idea  may  be  investigated  with  possible  operational 
testing  in  the  D-009  Apollo  experiment. 

Discussion 

Mr.  Hamer  points  out  that  the  range  measurement  method  he  has 
proposed  could  be  used  just  as  easily  with  the  Sun  and  a  planet  as 
with  the  earth  and  the  moon.  Also,  he  has  determined  that  after  the 
first  five  hours  from  injection  on  a  trcnslunar  trajectory,  the 
uncertainty  in  range  is  relatively  constant  at  about  35  kilometers  as 
long  as  the  star  used  is  within  plus  or  minus  30  degrees  of  the 
earth-moon-vehicle  plane. 

These  results  are  very  encouraging.  With  range  measurements  of 
this  accuracy,  the  manual  navigation  techniques  described  in  this 
thesis  are  of  much  greater  utility  in  cislunar  and,  indeed,  in 
interplanetary  space  than  they  would  have  been  if  angular  diameter 
measurements  were  used  exclusively  to  determine  range.  The  linearity 
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characteristics  of  the  method  allow  large  deviations  from  the  reference 
trajectory.  Thus,  barring  a  catastrophic  failure  such  as  an  uncon¬ 
trollable  propulsion  system,  this  method  is  flexible  enough  to  be 
incorporated  in  a  manual  navigation  scheme.  The  only  equipment 
necessary  would  be  the  space  sextant  with  attachment  for  body  center 
measurements,  a  time  history  of  the  reference  trajectory  (or 
trajectories,  if  several  sets  of  data  are  carried  for  greater 
flexibility),  and  a  computational  tool  for  use  in  the  solution  of 
Eq  (8-13). 


Sums  .ary 

The  precomputations  necessary  far  the  Hamer  method  of  range 
determination  are: 


1*  rey,  rmv,  rem,  A,  B,  and  F  as  functions  of  time  on  the 


reference  trajectory 


2.  3r 


ev 


‘•mv 


3A 


sin  A 


3.  3r 


ev 


3B 


rem  (cos  B  cot  A  “  sin  B) 


h  =  rev  x  rmy  (See  Appendix  E) 


:  i?  — 

'  •  OXU  JU  — 


H  hx  +  m  hy  +  n  hz 


/«  .  - 


voee  appennxx 


6‘  5  =  *ev  x  H 


(See  Appendix  E) 
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7.  k 


£  b„  +  m  b  +  n  b 
x  y  z 

£bx  +  mby  +  nbz 


(See  Appendix  E) 


All  of  these  values  must  be  precomputed  before  flight.  Inflight,  the 
astronaut  is  required  to  solve  only 


Ar  =  C,  AA  +  C„  AF 
ev  £  2 


(8-13a) 


where  Cj  and  C2  are  constants  which  are  precomputed  for  the  reference 
trajectory,  i.e., 


C 


1 


^rev 

OB 


(8-23) 


and 


Or 

k  — 1 — 
OB 


1/2 


(8-24) 


Sample  Problem.  The  following  is  a  sample  calculation  of  distance 
from  the  earth  using  the  Hamer  method.  The  reference  values  were 
taken  from  Ref  8. 

At  31.125  hours  from  injection  (about  halfway  to  the  moon), 
the  reference  trajectory  values  are: 


x  =  -  9233.1921 
ev 

y=„  =  211785.18 

CV 

z  =  116836.74 
ev 

rev  =  242051.73 
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x  =  -  113475.61 

mv 

y  =  -  101136.23 

mv 

z  =  -  49059.589 

mv 

r  =  159725.06 

mv 

r  =  369198.98 

em 

A  =  132.4675  degrees 

B  =  18.0877  degrees 

F  =  82.4794  degrees  (the  star  chosen  is 
Algenib) 

The  precomputation  yields 


and 


Cj  =  -  212234.26 


C  -  -  450855.29 
2 


An  actual  position  was  assumed  at 

x  =  -  15000 

ev 

y  =  220000 

ev 

z  =  122000 

ev 
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So  that 


ev 


252010 


which  is  approximately  10000  kilometers  off  the  reference  value .  (This 
position  was  chosen  to  test  the  linearity  characteristics  of  the  Hamer 
method.)  At  this  position,  the  actual  angular  measurements  would  be, 
assuming  10  arc-seconds  error  in  each: 

^actual  =  127.2446  degrees 


^actual  =  83  * 6540  degrees 


Applying  Eq  (8-13a) 


Ar  =  C.  AA  +  C0  AF 
ev  1  2 


or 


r  =  r  +  C,  AA  +  C„AF 

actu?l  ref  1  2 


yields 

r  ,  =  252155  kilometers 

actual _ • 

This  value  is  in  error  by  145  kilometers,  since  the  precise  assumed 
value  is  252010  kilometers.  This  represents  a  .06  per  cent  error. 
An  angular  diameter  measurement  would  have  been  in  error  by  more 
than  500  kilometers.  It  is  probable  that  the  error  in  the  Hamer 
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estimate  of  the  actual  range  is  primarily  due  to  nonlinearity  at  such 
a  large  distance  off  the  reference  trajectory.  The  coefficients 
C1  and  C2  are  fairly  large  also,  and  this  magnifies  the  measurement 

error  introduced.  In  any  case,  the  Hamer  method  is  a  vast  improvement 
over  the  angular  diameter  method. 
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IX ,  Aids  to  Manual  Navigation 

As  progress  is  made  in  the  exploration  of  space,  operational 
manual  navigation  and  guidance  systems  will  become  a  necessity.  The 
authors  of  this  thesis  believe  that  as  this  requirement  makes  itself 
felt  certain  "aids"  will  be  developed  to  enhance  the  performance  of 
these  systems. 

To  date,  no  concerted  effort  has  been  made  to  develop  a  general 
manual  system.  Those  who  have  worked  on  the  idea  have  been  limited 
by  the  constraints  of  the  hardware  available  at  the  present  time. 

(The  sum  total  c*  this  hardware  consists  of  the  USAF  Hand-held 
Space  Sextant  and  the  USAF  Space  Stadimeter.)  The  manual  navigation 
techniques  developed  thus  far  are  inflexible  due  to  these  constraints 
Many  of  the  methods  described  in  this  thesis  are  impractical  when 
viewed  from  the  standpoint  of  current  equipment.  The  authors  feel 
that  no  promising  method  or  approach  should  be  discarded  on  this 
basis.  Indeed,  the  equipment  must  be  designed  to  fit  the  best 
approach,  and  not  vice  versa. 

It  is,  therefore,  the  purpose  of  this  chapter  to  explore  the 
utility  of  three  categories  of  "aids"  to  manual  navigation;  the 
back-up  inertial  platform,  the  hand-held  mechanical  calculator,  and 
the  hand-held,  battery  operated,  digital  computer. 
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The  Back-up  Inertial  Platform 

A  back-up  inertial  platform  would  be  invaluable  in  determining 
spacecraft  position.  An  instrument  such  as  a  theodolite  mounted  on 
an  inertial  platform  is  capable  of  measuring  directly  the  right 
ascension  and  declination  of  the  vehicle  by  simply  sighting  the  earth, 
sun,  or  any  other  body  whose  coordinates  are  known.  No  computations 
would  be  necessary  to  obtain  this  information.  The  space  sextant 
could  be  used  to  align  the  platform  by  sighting  two  known  stars. 

Since  the  lines-of-sight  define  a  plane,  the  platform  can  be  oriented 
so  as  to  be  parallel  to  any  desired  reference  plane.  The  power  source 
could  be  electrical  (batteries)  or  mechanical  (hand-crank) . 

The  Hand-held  Mechanical  Calculator 

Many  of  the  manual  navigatj  n  techniques  described  in  this  report 
involve  several  computation  steps.  It  would  be  difficult  to  achieve 
consistent  accuracy  with  a  slide  rule.  Reducing  the  equations  to 
graphical  form  is  not  always  practical  because  of  the  volume  of  graphs 
necessary  for  flexible  operation.  A  better  approach  may  lie  in  the 
use  of  a  small,  highly  accurate  mechanical  calculator.  Such  a  device 
is  manufactured  by  the  Curta  Corporation,  of  Van  Nuys,  California. 
Figure  41  is  a  photograph  of  the  Curta  Calculator.  This  machine  is 
capable  of  performing  addition,  subtraction,  multiplication,  and 
division.  It  also  squares,  cubes,  and  extracts  square  roots.  The 
accuracy  of  the  Curta  Calculator  is  an  amazing  eight  significant 
figures!  (This  is  based  upon  the  limit  imposed  by  division. 
Multiplication  can  be  carried  to  fifteen  figures.)  The  diameter  of 
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Fig.  41 

The  Curta  Calculator 

the  machine  is  2  9/16  inch,  and  the  height  is  3  5/8  inch.  Most 

importantly,  the  weight  of  the  Curta  Calculator  is  but  12  1/2  ounces. 
This  hand-held  mechanical  calculator  combines  the  accuracy  of  a  desk 
calculator  with  the  portability  of  a  slide  rule.  The  authors  feel 
that  such  a  computational  tool  is  ideally  suited  as  an  aid  to  manual 
navigation. 
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The  Hand-held ,  Battery  Operated,  Digital  Computer 

The  hand-held,  battery  operated,  digital  computer  is,  in  the 
opinion  of  the  authors,  the  most  important  tool  in  the  future  of  man¬ 
ual  navigation  in  space.  With  the  proper  measuring  equipment,  such  a 
device  could  be  capable  of  solving  explicitly  the  manual  navigation 
problem.  The  most  promising  methods  explored  in  this  thesis  are 
ideally  suited  to  being  programmed  into  a  small  computer.  The  proper 
use  of  such  a  tool  is  paramount  in  determining  its  true  utility. 

Memory  space  is  the  prime  limitation.  It  would  be  unreasonable  to 
attempt  "brute  force"  techniques  that  require  excessive  iterations. 

The  hand-held  computer  should  not  be  used  to  duplicate  the  functions  of 
the  primary  system  computers.  It  should  not  be  pre-programmed  with 
inflexible  routines  that  cannot  be  altered  as  needed.  It  should  be 
designed  so  that  the  astronaut,  who  is  obviously  familiar  with  astro- 
dynamics  and  celestial  mechanics,  can  program  the  equations  he  needs 
under  any  circumstances.  Several  small  computers  are  in  existence 
today.  The  most  promising  investigated  is  the  Hewlett-Packard  9100A. 
This  computer  is  capable  of  performing  at  least  as  well  as  the  IBM  1620. 
It  can  be  programmed  with  magnetic  "credit  cards"  to  solve  problems 
up  to  and  including  numerical  integration.  The  present  weight  of  the 
9100A  is  40  pounds.  There  is  no  reason  to  believe  that  this  couldn't 
be  reduced  for  manual  navigation  and  guidance  applications. 

To  demonstrate  the  utility  of  a  small,  hand-held  computer  in 
manual  navigation,  the  Gibbs  method  of  orbit  determination  (see 
Chapter  V)  was  programmed  in  Fortran  II  and  run  on  the  IBM  1620. 

This  program  would  easily  fit  into  the  Hewlett-Packard  9100A. 
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Figure  42  is  a  reproduction  of  the  Gibbs  program  with  output.  This 
program  is  totally  flexible  in  that  it  may  be  applied  to  elliptic  and 
hyperbolic  trajectories.  The  output  includes  the  type  of  trajectory, 
the  eccentricity,  the  semimajor  axis,  the  time  past  perifocus  to  the 

A  A  A 

third  fix,  the  P,  Q,  and  W  orientation  vectors,  and  the  velocity 
components  at  the  third  fix.  The  three  input  fixes  were  obtained  from 
a  precision  trajectory  in  Ref  8.  .  1  of  the  output  values  are  within 

.01  per  cent  of  the  quoted  precision  (translunar)  trajectory  values. 
Since  the  program  is  written  in  general  form,  it  is  applicable  to  any 
inverse-square  field,  as  long  as  the  gravitational  constant  for  the 
attracting  body  is  known.  This  constant  must  be  input  with  the  three 
fixes.  The  output,  then,  is  in  the  units  of  the  gravitational  constant 
and  the  position  fixes.  For  example,  in  Fig.  42,  the  gravitational 
constant  and  fixes  are  in  the  kilometer-second  system  of  units.  There¬ 
fore,  the  output  distances  are  in  kilometers  and  the  velocities  are  in 
KM/SEC. 

Since  position  and  velocity  are  known  at  fix  three,  the  position 
and  velocity  at  any  time  thereafter  can  be  projected.  Figure  43  is 
a  reproduction  of  a  program  which  uses  the  universal  formulas  developed 
by  Battin  in  Ref  4.  The  output  is  position  and  velocity  at  hourly 
intervals  past  the  third  fix.  Again,  a  general  form  is  used  for 
flexibility,  and  the  units  of  the  input  determine  the  units  of  the 
output.  (Since  the  output  of  the  Gibbs  program  is  input  to  this 
program,  its  output  is  in  kilometers  and  KM/SEC.)  This  information 
could  be  used  to  initiate  a  guidance  scheme. 
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C  C  GIBBS  METHOO  OF  ORBIT  DETERMINATION— HORRIGAN  AND  WAlSH 
PI=3. 1415927 

READtXl » Y1 »Z1 »X2 » Y2  »Z2  »X3*Y3»Z3  »EMU 

A=SQRT( (Y1*Z3-Y3*Z1)**2+(Z1*X3-X1*Z3)**2+(X1*Y3-Y1*X3)**2) 
B»SQRT( <Y2*Z  3-Y3*Z2)*#2+(Z2*X3-X2*Z3)**2+(X2*Y3-Y2*X3)**2) 
C»SQRT( (Y2*Z1-Y1*Z2)**2+(Z2*X1-X2*Z1)**2+(X2#Y1-Y2*X1)**2) 
C1«B/A 
C3*C/A  . 

R1=SQRT(X1**2+Y1##2+Z1**2) 

R2*J 2RT<  X2**2+Y2*«2+Z2**2 ) 

R3-S  3RT < X3**2+Y3**2+Z3*#2 ) 

P«(C1*R1+C3*R3-R2)/(C1+C3-1.) 

WX=(Y1*Z3-Y3*Z1)/A 

WY*(Z1*X3-X1*Z3)/A 

WZ«(X1*Y3-Y1.*X3)/A 

EQX= ( X3* ( P-R 1 ) -X 1* ( P-R3 ) ) /A 

EQY=(Y3*(P-R1)“Y1*(P-R3) ) /A 

EQZ* (Z3* ( P-Rl ) -Zl# ( P-R 3 ) ) /A 

E35QRT(EQX*#2+EQY**2+EQZ**2 ) 

QX*EQX/E 

QY=EQY/E 

QZ-EQZ/E 

AXIS=P/(1.-E**2) 

PX3QY*WZ-QZ*WY 
PY=QZ*WX-QX*WZ 
PZ3QX*WY-WX*QY 
YW33QX*X3+QY*Y3+QZ*Z3 
IF( AXIS) 1»1»2 

2  COSE33 ( AXIS— R3 }/ { AXIS*E) 

SINE3=YW3/<AXIS*SQRT(1.-E**2n 
GAMMA=ATAN( ( E/ (SORT  < 1»-E**2 ) ) )*SINE3 ) 

PUNCH  100 

100  FORMAT ( 26HTHE  TRAJECTORY  IS  ELLIPTIC*//) 

IF (C0SE3 ) 5  *5  *3 

5  EANOM=ATAN (SINE3/C0SE3 )+PI 
GO  fO  11 

3  EANOM*AT  AN ( SI NE3 /C0SE3 ) 
ll  EN-EMU/(AXIS**1*5) 

DELTT=UEANOM-E*SINE3)/EN)/3600. 

GO  TO  15 

1  SINHF=-YW3/(AXIS*SQRT(E**2-1.) > 

GAMMA=SQRT( (P**2)/(R3*(Z.*P+R3*(E**2-1.) > ) ) 

PUNCH  101' 

101  FORMAT ( 28HTHE  TRAJECTORY  IS  HYPERBOLIC*//) 

F=LOG( S INHF+5QRT  ( I SINHF*#2 )  +  l ♦ ) > 

EN-»EMU/(  (—AXIS  )  -**1*  5  ) 

DELTT«=  <  ( E*SI  NHF-F )  /EN)  /3600» 

15  RTDOT= ( EMU/R3) *SQRT <  AXI S*ABS ( l.-E**2 ) ) 

ROOT  =RTDOT*S I N ( GAMMA ) /COS ( GAMMA ) 

VX= ( 1D0T*X3+RTD0T* ( WY*Z3-WZ*Y3 ) ) /R3 
VY=  ( RD.OT#Y3+RTDOT*  ( WZ*X3-WX*Z3 )  )  /R3 
VZ= ( RD0T*Z3+RTD0T* ( WX*Y3-WY*X3 ) ) /R3 
PUNCH  102»E 

102  FORMAT ( 20HTHE  ECCENTRICITY  IS  #F7.4*//) 

PUNCH  103 » AX  IS 

103  FORMAT( 23HTHE  SEMI-MAJOR  AXIS  IS  »F12.4*//) 

PUNCH  104,DELTT 


Fig.  42 

Typical  Gibbs  Method  Program 
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104  FORMAT I40HTHE  TIME  PAST  PERIFOCUS  TO  FIX  THREE  IS  »F12.4*//) 

PUNCH  105 

105  FORMAT! 16HPX.PY .AND  PZ  ARE) 

PUNCH  106.PX.PY.PZ 

106  FORMAT! 3F12. 4*//) 

PUNCH  107 

107  FORMAT! 16HQX.QY. AND  QZ  ARE) 

PUNCH  1Q6.QX  *QY *QZ 

PUNCH  108 

108  FORMAT ( 16HWX *WY. AND  W Z  ARE) 

PUNCH  106*WX*WY»WZ 

PUNCH  109 

109  FORMAT! 54HTHE  X.Y.AND  Z  VELOCITY  COMPONENTS  AT  THE  THIRD  FIX  ARE) 
PUNCH  106»VX*VY  »VZ 

STOP 

END 

17o00«141  13826.706  8591.5583  20284.239  66965.325  38219.273 

15797.195  106257.19  59747.541  631.35093 


THE  TRAJECTORY  IS  ELLIPTIC 

THE  ECCENTRICITY  IS  .9739 

THE  SEMI-MAJOR  AXIS  IS  247285.1800 

THE  TIME  PAST  PERIFOCUS  TO  FIX  THREE  IS  10.5060 


PX.PY.AND  PZ 

ARE 

.2694 

-.8485 

-.4555 

CX.QY.AND  QZ 

ARE 

.9619 

.2136 

.1709 

WX.WYtAND  WZ 

ARE 

-.0478 

-.4842 

.8737 

THe  X*Y*AND  Z  VELOCITY  COMPONENTS  AT  THE  THIRD  FIX  ARE 
-.3010  1.9199  1.0475 


Fig.  42  (Continued) 
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C  C  PROJECTED  POSITION  AND  VELOCI TY—HORRIGAN  AND  WALSH 
READ »X»Y#Z#VX»VY»VZ» EMU 
TT*3600. 

B*l. 

4  XX*1. 

DI ST-SQRTI X*X+Y*Y+Z*Z I 
VSQR »VX*VX+VY*VY+VZ*VZ 
ALPh A*2 • /D I ST-VSQR/ ( EMU**2 ) 

1  IF(ALPKA)8»7»2 

8  D=SQRT(-ALPHA*(XX**2) ) 

S= ( t  EXP{ D)-EXP(-D) ) /2.-D) / ( D**3 ) 

C=( (EXP{D)+EXP(-D))/2«-l.)/(D**2> 

GO  TO  3 
7  S=l./6. 

C=.5 
GO  TO  3 

2  E=SQRT ( ALPHA*  <  XX**2 } ) 

S= ( E— SIN { E ) J / { E#*3 ) 

C= ( l.-COS( E) ) / (E**2 ) 

3  Tl°<  ( X*VX+Y*VY+Z*V2 ) * { XX**2  )*C)/( EMU**2 ) 

T2={ ( l.—ALPHA*DIST ) *  <  XX** 3 I+S I/EMU+ ( DI ST*XX  > /EMU 
T=T1+T2 

IF (ABS(TT— T )-»01)6*6.5 

5  Q=(X*VX+Y*VY+Z*VZ)/(EMU**2) 

DER I V=Q* ( XX-ALPHA* ( XX**3 ) *S } + ( ( 1 . -D ( ST*ALPHA ) * ( XX**  2 ) *C+D I ST ) /EMU 
XX=XX— ( T—TT ) /DERI V 
GO  TO  1 

6  RX=< l.~< <XX**2)/DIST)*C>*X+CTf~<XX**3)*S/EMU)*VX 
RY=<l.-< (XX**2)/DIST)*Ci*Y+<TT-(XX**3)*S/EMU)*VY 
RZ= ( 1.- ( ( XX**2 ) /DI ST ) *C ) *Z+<  TT- ( XX**3 ) *S/EMU ) *VZ 
R=SQRT l RX*RX+RY*RY+RZ*RZ ) 

A=EMU 

VXX= ( ( 1 . / <  R*D I ST ) } * ( ALPHA*  C  XX**3 ) *S-XX ) ) *X*A+ ( 1 . - ( ( XX**2 ) /R ) *C ) *VX 
VYY= ( ( 1 */ ( R*DI ST ))  *  ( ALPHA* < XX**3 ) *S-XX ) ) *Y*A+ ( 1 .- ( ( XX**2 )/R)*C) *VY 

VZZ* ( ( 1. / ( R*DIST ) )  *  < ALPHA* ( XX**3 )*S-XX ) ) *Z*A+ ( l.-C  C XX**2 ) /R ) *C ) *VZ 
TTT*B 

PUNCH  100»  TTT 

100  FORMAT ( 38HPOSI T ION  AND  VELOCITY  COMPONENTS  AT  T=»F7.4»//) 

PUNCH  101.RX.RY.RZ 

101  FORMAT (3F12.4.// J 
PUNCH  101 »VXX*VYY. VZZ 
B=R+1. 

X=I<X 
Y»RY 
Z«RZ 
VX*V  XX 
VY=V YY 
VZ=VZZ 
GO  TO  4 
END 

15797.195  106257.19  59747.541  -.3010  1.9199  1^0475  631,35093 


Fig.  43 

Projected  Position  and  Velocity  Program 
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C  C  PROJECTED  POSITION  AND  VELOCITY 


POSITION  AND 

VELOCITY  COMPONENTS  AT  T* 

1.0000 

14693.4090 

113027.0500 

63438.8700 

t 

-.3118 

1.8426 

1.0041 

POSITION  AND 

VELOCITY  COMPONENTS  AT  T* 

2.0000 

13555.2520 

119533*5600 

66982.3800 

-.3202 

1.7733 

•  9652 

POSITION  AND 

VELOCITY  COMPONENTS  AT  T* 

3.0000 

12390.14 40 

125803.2800 

70393.1990 

-.3268 

1.7108 

•  9302 

POSITION  AND 

VELOCITY  COMPONENTS  AT  T« 

4.0000 

11203.8160 

131858.2800 

73683.8680 

-.3320 

1.6539 

.8984 

FDsITION  AND 

VELOCITY  COMPONENTS  AT  T* 

5.0000 

10000.7620 

137717.1100 

76864.9120 

-.3362 

1.6017 

•  8693 

POSITION  AND 

VELOCITY  COMPONENTS  AT  T« 

6.0000 

8784.5643 

143395.5600 

79945.2750 

-.3394 

1.5536 

.8424 

POSITION  AND 

VELOCITY  COMPONENTS  AT  T* 

7.0000 

7558.1158 

148907.2500 

82932.6660 

-.3419 

1.5090 

.8176 

POSITION  AND 

VELOCITY  COMPONENTS  AT  T= 

8.0000 

6323.7722 

154263.9600 

85833.7390 

-.3438 

1.4674 

.7944 

i 

Fig.  43  (Continued) 
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C  C  ORBIT  GEOMETRY  FROM  ANGULAR  RATE  DATA—HORRIGAN  AND  WALSH 

READ*DT12*DT13*TDl»TD2»TD3  *EMU 
STD1=SQRTCTD1) 

STD2aSQRT(TD2 ) 

STD3aSQRT(TD3) 

Q=STD2*SIN<DT13)-STD3*SIN(DT12)+STD1*(SIN(DT12)-SIN<DT13)) 
Ra-(STD3-STD1)*C0S<DT12)+(STD2-STD1)*C0S(DT13)-(STD2-STD3) 
SaSTDl*SIN (DT13-DT12 )+STD3*SI N( DT12 ) -STD2*SIN (DT13) 

EaSQRTI <R/S)**2+!Q/S>**2) 

TlaATAN( R/Q) +3 .1415927 

P= << EMU* << 1. +E*COS< T1 )) **4 ) )/<TDl**2>  >**.33333333 
AaP/ABS<l.-E**2) 

RDOTaE*SQRT ( EMU/ ( A*ABS { 1 .-E**2 > ) >  *S I N ( T 1 ) 

Rl=<  <EMU*P)**.25)/STD1 
TlaT 1*57 *295 7795 
RTD0T»R1*TD1 
IF( A) 1*1 *2 

1  PUNCH  100 

100  FORMAT ( 28HTHE  TRAJECTORY  IS  HYPERBOLIC*//) 

GO  TO  3 

2  PUNCH  101 

101  FORMAT l 26HTHE  TRAJECTORY  IS  ELLIPTIC*//) 

3  PUNCH  1C2*E 

102  FORMAT ( 20HTHE  ECCENTRICITY  IS  *F7.4»//) 

PUNCH  10 3.  A 

103  FORMAT! 23HTHE  SEMI-MAJOR  AXIS  IS  *F12.4*//) 

PUNCH  104.T1 

104  FORK  AT  <  39HTHE  1RUE  ANOMALY  AT  ANGULAR  FIX  ONE  IS  *F12.4»//) 

PUNCH  105.ROOT 

105  FORMAT ( 42HTHE  RADIAL  VELOCITY  AT  ANGULAR  FIX  ONE  IS  »F12.4,//) 
PUNCH  106.RTD0T 

106  FORMAT (45HTHE  TRANSVERSE  VELOCITY  AT  ANGULAR  FIX  ONE  IS.F12.4*//) 
PUNCH  107»R1 

107  FORMATI42HTHE  RADIAL  DISTANCE  AT  ANGULAR  FIX  ONE  IS  *F12.4*//) 
STOP 

END 

,138862  .2271076  .000010214764  .0000059230245  .0000040938722  62750.717 


THE  TRAJECTORY  IS  ELLIPTIC 

THE  ECCENTRICITY  IS  .9068 

THE  SEMI-MAJOR  AXIS  IS  104948.6000 

THE'  TRUE  ANOMALY  AT  ANGULAR  FIX  ONE  IS  138.3645 

THE  RADIAL  VELOCITY  AT  ANGULAR  FIX  ONE  IS  1.1051 

THE  TRANSVERSE  VELOCITY  AT  ANGULAR  FIX  ONE  IS  .5911 

THE  RADIAL  DISTANCE  AT  ANGULAR  FIX  ONE  IS  57870.6980 


Fig.  44 

Orbit  Determination  Program  Using  Angular  Rate  Data 
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If  a  complete  position  fix  is  not  available,  the  trajectory  may 
be  determined  by  knowledge  of  the  angular  position  of  the  vehicle  and 
three  true  anomaly  rates  (see  Chapter  VII).  Figure  44  is  a  reproduction 
of  a  program  which  utilizes  angular  rate  data.  In  this  case,  only  the 
true  anomaly  rates  and  differences  in  true  anomaly  are  known.  The 
orbit  geometry  can,  therefore,  be  determined;  but  not  the  inertial 
orientation.  The  program  can  be  altered  to  provide  this  information 
if  angular  position  data  is  available.  The  units  of  the  input  and 
output  are  nautical  miles  and  NM/SEC. 

Discussion 

This  chapter  has  been  a  discussion  of  some  important  aids  to 
manual  navigation.  This,  by  no  means,  exhausts  the  list  of  possi¬ 
bilities.  Total  commitment  to  the  development  of  manual  navigation  and 
guidance  systems  will  bring  forth  many  more  new  ideas.  ;  Flexibility 
and  accuracy  are  the. bywords  for  the  future.  Equipment  must  be  designed 
with  this  in  mind.  Methods  cannot  be  limited  just  because  hardware  is 
not  available.  It  must  be  made  available. 
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X.  Conclusions  and  Recommendations 


Introduction 

This  chapter  completes  the  theoretical  development  of  this  report 
with  general  conclusions  and  recommendations  for  follow-on  research. 
Many  interesting  possibilities  remain  to  be  examined  and  some  of  these 
topics  are  listed  under  "Recommendations."  The  "Conclusions"  para¬ 
graph  lists  some  of  the  important  determinations  resulting  from 
research  associated  with  this  thesis. 

Conclusions 

The  following  conclusions  are  offered  as  the  result  of  experience 
gained  in  preparing  this  report: 

1.  The  AH  method,  a  graphical  method  of  obtaining 
eccentricity,  is  feasible  for  high  eccentricity 
elliptical  trajectories.  A  major  drawback  is  the 
need  for  prior  knowledge  of  the  semimajor  axis. 

2.  The  numerical  differentiation  approach  for  computing 
the  geometric  elements  is  sensitive  to  measurement 
accuracy.  But  the  method  is  simple,  explicit,  and 
entirely  adequate  if  accurate  values  of  range  are 
available.  The  method  appears  to  be  especially  well 
suited  for  near  earth  orbits. 

3.  Differential  correction  procedures  are  impractical 
unless  a  back-up  computer  is  available.  If  such  a 
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computer  is  provided,  the  differential  correction  iw  thod 
should  be  exploited. 

4.  Position  fixes  can  he  obtained  usir.to  straightforward 
procedures.  The  simultaneous  measurement  requirement, 
frequently  cited  as  a  major  problem  in  obtaining  a 
three-dimensional  fix,  can  be  resolved  using  procedures 
similar  to  those  presently  used  in  astronomy  and 
celestial  aircraft  navigation.  With  position  fixing 
capability,  the  astronaut  possesses  a  powerful  and 
flexible  tool  that  can  be  used  for: 

(a)  orbit  determination 

(b)  cross-check- of  the  primary  navigation 
computer 

(c)  manual  guidance  procedures  (yet  to  be 
developed) 

5.  The  orientation  elements  are  much  easier  to  obtain 
than  the  geometric  elements.  Two  angular  positions 
and  knowledge  of  true  anomaly  at  some  point  of  the 
orbit  are  al ..  that  is  required  to  completely  define 
the  orientation  parameters. 

6.  Range  can  be  determined  optically  with  much  better 
precision  than  previously  thought.  The  accuracy 

of  the  Hamer  method  of  measuring  range  (Chapter  VIII) 
is  far  better  than  the  semi-diameter  method. 

7.  Manual  navigation  aids  such  as  small  mechanical 
calculators  and  back-up  computers  that  can  be 
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programmed  in  fligh  .  art  practical  necessities. 

Recommendations  for  Follow-on  Theses 

The  following  areas  are  suggested  as  worthy  endeavors  for  future 
research : 

1.  Further  study  of  graphical  methods  for  obtaining 

,  orbital  parameters .  Extend  the  AH  method  to  hyperbolic 
trajectories.  Study  the  behavior  of  the  AH  plots 
around  an  entire  high-eccentricity  ellipse  (perigee  to 
perigee) .  Examine  other  graphical  solutions  to  the 
orbit  determination  problem. 

2 .  Extended  sttidy  of  methods  for  determining  range .  The 
range  measurement  accuracy  is  considerably  improved 
using  the  methods  of  Chapter  VIII.  However,  if  a 
suitable  reference  trajectory  is  not  available,  other 
range  measurement  procedures  must  be  provided.  Several 
alternate  methods  are  discussed  in  Ref  4:223. 

3.  Continued  study  and  development  of  orbit  determination 
methods.  Many  of  the  classical  orbit  determination 
methods  are  simplified  by  the  fact  that  the  observer  is 
on  the  spacecraft  and  not  on  the  surface  of  the  earth. 
Several  excellent  orbit  determination  schemes  are 
suggested  in  Ref  20:683,  Ref  13,  and  Ref  3:Ch  1.  All 
of  the  published  methods  should  be  investigated.  The 
goal  should  be  to  assemble  new  orbit  computation 
methods  using  the  good  points  of  previous  schemes. 
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Special  attention  should  be  given  to  methods  that 
require  no  direct  measurement  of  distance  (such  as 
the  method  of  Chapter  Vti) . 

4.  Error  Analyses.  A  number  of  methods  discussed  in 
this  thesis  require  comprehensive  error  analysis. 

For  example,  the  Polaris  fixing  procedure  must  be 
investigated  to  see  what  corrections,  if  any,  are 
required  to  compensate  for  the  fact  that  Polaris 
deviates  from  the  pole  position.  Other  topics  of 
the  thesis  require  error  analysis  before  final 
adoption  as  manual  navigation  procedures. 

5.  Back-up  Computer  Techniques.  A  small  back-up  computer 
will  greatly  enhance  manual  navigation  capability. 

But  the  capability  of  such  a  computer  can  be  easily 
wasted  on- needless  iterative  techniques.  Several 
short  programs  were  presented  in  this  thesis.  Other 
programs  and  procedures  should  be  devised — the  purpose 
being  the  most  efficient  use  of  a  limited  memory 
computer. 

6.  Application  of  Differential  Correction.  The  differential 
correction  equations  are  complex  but  extremely  useful 

if  a  back-up  computer  is  available.  It  appears  that 
a  navigation  scheme  incorporating  redundant  data  could 
be  constructed  using  the  differential  correction  concept. 
The  use  of  redundant  data  is  desirable  since  the  effect 
of  measurement  error  is  reduced  considerably. 
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7.  Systematized  Navigation  Procedures.  After  the  suitability 
of  a  navigation  scheme  is  established,  the  astronaut  must 
be  provided  systematic  procedures  in  flight  manual  form. 
The  flight  r. anual  should  consist  of  charts,  worksheets, 
and  operational  procedures  arranged  to  render  the 
solution  of  the  navigation  problem  to  step-by-step  form. 
Compilation  of  such  a  flight  manual  is  suggested  as 

a  follow-on  thesis  project. 

8.  Manual  Guidance.  After  the  orbit  is  known,  the  astronaut 
must  have  the  capability  to  determine  what  corrections 
are  required  to  make  the  spacecraft  go  where  he  wants  to 
go.  The  problem  of  a  manual  abort  after  translunar 
injection  is  of  prime  concern.  At  first  glance,  it 
appears  the  Gibbsian  procedure  (Chapter  V)  could  be 
molded  into  an  explicit  guidance  scheme.  In  any  event, 
the  guidance  problem  is  complex  and  could  well  support 
several  thesis  topics. 
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XI .  Concluding  Remarks 

As  stated  in  Chapter  I,  the  purpose  of  this  thesis  is  co  lay  a 
foundation  for  the  development  of  an  operational  manual  navigation 
system.  Several  independent  methods  of  obtaining  the  necessary  orbital 
parameters  have  been  presented.  The  combination  of  procedures  selected 
for  use  in  an  operational  system  depends  upon  several  factors,  and 
most  importantly  upon  the  answers  to  the  four  questions  posed  in  the 
Introduction  (Chapter  1:4).  The  authors  believe  that  these  questions 
should  be  answered  as  follows: 

1.  What  computational  tools  are  available? 

Answer  -  A  hand-held,  battery  operated,  computer 
capable  of  flexible  operation,  and  a 
hand-held  mechanical  calculator  for 
incidental  computations. 

2.  What  types  of  measurements  are  possible? 

Answer  -  Equipment  and  procedures  will  be  made 
available  to  allow  complete  position 
fixing  using  star-body  and  range  measure¬ 
ments. 

3.  How  much  time  is  allowable  for  the  completion  of  the 
procedures? 

Answer  -  At  least  five  hours  for  highly  eccentric 

trajectories,  and  correspondingly  less  for 
lower  eccentricity  orbits. 
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4.  What  is  the  required  accuracy? 

Answer  -  This  depends  upon  the  use  to  which  the 
orbital  parameters  are  put.  Different 
requirements  will  exist  for  re-entry 
navigation  and  guidance  than  for  mid- 
course  navigation  and  guidance.  The 
best  accuracy  attainable  should  be  a 
goal,  however,  since  the  same  manual 
system  will  probably  be  used  throughout 
the  mission. 


If  the  above  answers  are  acceptable,  the  manual  navigation  problem 
is  solved.  An  operational  system  consisting  of  position  fixes  input  to 
the  Gibbs  method  of  orbit  determination  programmed  into  a  back-up 
computer  is,  in  the  author's  opinion,  the  most  promising  combination 
devised  to  date.  It  is  completely  flexible;  it  involves  no  approxi- 
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Appendix  A 

The  Computer  Program 

This  Appendix  consists  of  a  reproduction  of  the  double-precision 
computer  program  developed  to  investigate  the  Delta-H  method  of 
eccentricity  determination  and  to  provide  accurate  data  for  use 
throughout  this  thesis.  A  listing  of  the  variables  is  included  to  aid 
the  reader  in  following  the  logic  involved. 
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EMU  =  (gravitational  constant) 

RPLAN  =  earth's  radius  in  NM 

AXSS  =  semimajor  axis  of  ellipse  in  NM 
AXI  =  (AXSS) 1  * 5 

EN  =  mean  motion 

TAU  =  orbital  period  in  hours 

ETAU  =  number  of  2.5  minute  increments 

E>."C  =  eccentricity 

TEE  =  time  in  seconds  from  perigee 
ECAN  =  eccentric  anomaly  in  radians 
RADI  =  perigee  radius 
VEL  =  velocity  magnitude  in  NM/SEC 
ALT  =  altitude 
EMAN  =  mean  anomaly 
RAD  =  radial  distance 
CONT  =  time  in.  minutes  from  perigee 
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DELTB  =  AH2 

ECA  =  eccentric  anomaly  in  degrees 

VELO  =  ratio  of  velocity  to  perigee  velocity 
The  program  that  follows  generates  data  for  high  eccentricity 
ellipses  with  a  five  minute  timing  interval  for  AH1  and  AH2.  A  portion 
of  the  resulting  output  is  also  included.  Total  ninning  time  on  the 
IBM  7094  digital  computer  is  approximately  two  minutes. 


I  JOB  0*3*5Q00  68-514  WALSH»RC  afit-se 

$  I  ft JOB 
SIrFTC  main 

C  GEOCENTRIC  5  MINUTE  DELTA-H  FOR  HIGH  ECCENTRICITY 

DOUBLE  PRECISION  ECAN* ALT  » VEL  *EMU*T  AU» ECC* RADI *RIGHT  »CONT ♦ 
1VEL0*RPLAN»AXI *ETAU * EMAN * RAD* DELTA #ECA»AXSS»EN*T EE *ALEFT  *DELTB 
DIMENSION  ECAN (85) »ALT(85) * VEL ( 83 ) 

EMU=DSQRT( 62750.717 ) 

RPLAN-3440.1728 
AXSS=105000.00 
DO  12  1=1*7 

AXI =DSQRT ( AXSS*AXSS*AXSS ) 

EN^EMU/AXI 

TAU=6.28318530/!EN*3600.) 

ETAU=210./2.5 

N=ETAU 

ECC=0«9068 

DO  11  J=l*7 

TEE=150.0 

WRITE(6*1) 

1  FORMAT ( 1H1*9X*4HAXIS* 13X* 12HECCEN1 RICI TY*7X»6HPERIOD*// ) 

WRITE! 6 *2)  AXSS»ECC*TAU 

2  FORMAT !1X*3F17*4»//) 

ECAN ! 1 ) =0* 

RAi)l=AXSS-AXSS*ECC 

VEL i 1 ) =EMU*DSQRT ( ( 2 . /RADI )  -  ( 1 . /AXSS } ) 

ALT(1)=RAD1-RPLAN 
DO  13  K*2»N 
EMAMEN*TEE 
ECAN ( K) =ECAN t  K-l ) 

14  ALEFT=DSIN(ECAN!K)) 

RIGHT= ( 1 . /ECC) *( ECAN ( K ) -EMAN > 

IF !ABS!ALEFT-RIGHT) -*0000001) 6*6*9 
9  ECAN (K)=  ECAN  <  K ) - <  ALEFT-R I GHT ) / < DCOS !  ECAN  l  K.  >  )  - 1 .  / ECC) 

GO  TO  14 

6  RAD=AXSS-AXSS*ECC*DCOS ( ECAN ( X  > ) 

ALT { K ) =RAD-RPLAN 

VEL!K)=EMU*DSQRT!  !2./RAD)-U*/AXSS)  ) 

10  TEE=TEE+150.0 

WRI TE( 6*3 ) 

3  FORMAT ( 2X»4HTIME *5X»8HALT ITUDE»5X*8HDELTA  Hi *5X»8HDELTA  H2*3X* 
111HECC  ANOMALY »7X .4HV/VP* // } 

M-N-4 
CONT=0»0 
DO  13  L=1»M 
DELTA=ALT ( L ) -ALT  !  L+2  > 

DELTB=ALT ( L+2 ) -ALT ( L+4 ) 

ECA=ECAN(L)*57.2957795D0 


4 

13 

U 

12 

30 

SEqF 


VELO*VEL ( L ) /VEL ( 1 ) 

WRITE! 6*4)  CONT*ALT(L) *DELTA*DELTB*ECA»VELO 
FORMAT !”6*1»5F13*3) 

CONT=CONT+2.5 
ECC=ECC+0*  01 


AXSS=AXSS+1000.00 

STOP 

END 
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Axis 

ECCENTRICITY 

PERIOD 

105000.0000 

0.9068 

237.0566 

TIME 

ALTITUDE 

DELTA  HI 

DELTA  H2 

ECC  ANOMOLY 

V/VP 

0. 

6345.827 

“26.688 

-79.475 

0. 

1.000 

2.5 

6352.508 

“53.228 

-105.291 

0.679 

1.000 

5.0 

6372.515 

"79.475 

-130.547 

1.357 

0.999 

7.5 

6405.737 

-105.291 

-155.128 

2.033 

0.997 

10.0 

6451.991 

-130.547 

-178.933 

2.70  * 

0.994 

12.5 

6511.028 

-155.128 

-201.875 

3.376 

0.991 

15.0 

6582.538 

-178.933 

-223.884 

4.041 

0.988 

17.5 

6666.156 

-201.875 

-244.906 

4.701 

0.983 

20.0 

6761.471 

-223.884 

-264.902 

5.356 

0.978 

22.5 

6868.031 

-244.906  ‘ 

-283.850 

6.004 

0.973 

25.0 

6985.355 

-264.902 

-301.739 

6.644 

0.967 

ff.5 

7112.937 

-283.850 

-318.570 

7.278 

0.961 

30.0 

7250.257 

-301.739 

-334.355 

7.903 

0.955 

32.5 

7396.787 

-318.570 

-349.115 

8.521 

0.948 

35.0 

7551.996 

-334.355 

-362.878 

9.130 

0.941 

37.5 

7715.356 

-349.115 

-375.677 

9.730  ■" 

0.933 

40.0 

7886.350 

-362.878 

-387.550 

10.321 

0.926 

42.5 

8064.471 

-375.677 

-398.539 

10.903 

0.918 

45.0 

8249.228 

-387.550 

-408.686 

11.476 

0.911 

47.5" 

8440.148 

-398.539 

-418.035 

12.C39 

0.903 

50.0 

8636.778 

-408.686 

-426.630 

12.594 

0.895 

52.5 

- 883§’.687 

-418.035 

-434  .517 

13.140 

0.867 

55.0 

9045.464 

-426.630 

-441.737 

13.676 

0.879 

57.5 

9256.722 

-434.517 

-448.334 

14.204 

0.872 

60.0 

9472.094 

-441.737 

-454.347 

14.723 

0.864 

6?.  5 

9691.239 

-448.334 

-459.816 

15.233 

0.356 

65.0 

9913.832 

-454.347 

-464.778 

15.735 

0.848 

67.5 

10139.572 

'  -459.816 

-469.268 

16.228' 

0.341 

70.0 

10368.178 

-464.778 

-473.320 

16.714 

0.833 

72.5 

10599.388 

-469.268 

-476.964 

17.191 

0.826 

75.0 

10832.956 

-473.320 

-480.231 

17.660 

0.819 

77.5 

11068.656 

-476.964 

-483.148 

18.122 

0.812 

80.0 

11306.276 

-480.231 

-485.740 

18.576 

0.804 

82'.  5 

11545.620 

-483.148 

-488.032 

19.C23 

0.798 

85.0 

11786.507 

-485.740 

-490.045 

19.463 

0.791 

S7.5 

12028.767 

-488.032 

-491.802 

19.896 

0.784 

90.0 

12272.246 

-490.045 

-493.321 

20.322 

0.777 

92.5 

12516.799 

-491.802 

-494.620 

20.741 

0.771 

95.0 

12762.292 

-493.321 

-495.715 

21.155 

0.765 

_9T.5 

13008  .'601 

-494.620 

-496.622 

21.562 

0.758 

100.0 

13255.613 

-495.715 

-497.356 

21.962 

0.752 
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Appendix  B 

Derivation  of  Equations  for  Geometric 
Elements  Using  _r  and  i? 

The  following  notation  is  used  in  the  derivation: 
r  =  vehicle  position  vector 
x, y,z  =  coordinates  of  vehicle 

r-  \x\ 

r  =  radial  velocity 

r  =  radial  acceleration 

AAA 

i,j,k  =  unit  vectors  in  x,  y,  z  directions 
u  =  earth  gravitational  constant 
S  =  1 r |  =  magnitude  of  vehicle  velocity 
a  =  semimajor  axis 
E  =  eccentric  anomaly 
e  =  eccentricity 
0  =  true  anomaly 

Z  =  semilatus  rectum 
h  =  angular  momentum 
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Starting  with 

r2  =  x2  +  y2  +  z2 

and  differentiating  with  respect  to  time  obtain 

•  •  •  • 
rr  =  xx  +  yy  +  zz 

Taking  the  derivative  of  rr  leads  to 

r2  +  rr  =  (xx  +  yy  +  zz)  +  (x2  +  y2  +  z2) 

Now  x2  +  y2  +  z2  =  S2  and 

x  =  -  yr/r3  ,  y  yy/r3  ,  z  -  -  yz/i 

then 


r2  +  r r  =  S2  - 


9  2  2 

201  +  OL  + 


=  S2  -  — r  (x2  +  y2  +  z2) 
rJ 


=  S2  -u. 

r 


From  the  energy  equation 


S2  =  y 


2  _  1 
r  a 


r^+rr  =  -  &  -  s. 


2y_  _  y_  _  y_ 
r  a  r 


r2  +  rr  =  —  -  — 


(B-6) 


(B-7) 

(B-8) 


•  -f'i 


wr 


\ 


> 

I 
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Therefore 


a 


(B-9) 


a  relatively  simple  expression  for  the  semimajor  axis  a. 

Two  more  expressions  are  required  to  obtain  the  eccentricity  and 
eccentric  anomaly.  The  first  expression  is  the  equation  of  the 
conic  section 

r  =  a  (1  -  e  cos  E)  (B-10) 

from  which 


e  cos  E  =  1  -  -  (B-ll) 

a 

The  second  expression  is  derived  from  the  conic  section  equation 
involving  the  true  anomaly  0 


1  +  e  cos  0 

where  A  is  the  semilatus  rectum. 

The  first  derivative  yields  r 

r  =  1  e  (sin  Q 

(1  +  e  cos  0)2 


(B-12) 


(B-13) 
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then  multiplying  and  dividing  by  Z 


Z2  e  (sin  8)  6  _  r2  0  e  sin  6 

Z  (1  +  e  cos  0)2  ^ 


(B-14) 


However,  r2  0  is  the  angular  momentum  per  unit  mass,  i.e.,  ' 


h  =  r2  0  =  J^Z 


(B-15) 


so 


*  -  M  e  sin  0 

Z 


e  sin  0 


(B-16) 


from  which 


sin  0  = 


ft 


(B-17) 


The  orbit  plane  coordinate  system  is  shown  in  Fig.  45.  In  this 
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coordinate  system 


yw  “  r  sin  6 


•i/I 

_  r  r  V  u 


(B-18) 


and  a Is  j 


'w 


=  a  /j.  -  e2  sin  E 


(B-19) 


then 


/],  -  e2  sin  E  = 


_  r  rVn 


(B-20) 


For  an  ellipse,  £  =  a  (1  -  e2) ,  therefore 


*/T-e2  sin  E  = 


r  r  ^1  -  e2 
e  /y* 


(B-21) 


then 


e  sin  E  =  — r 


/ya” 


(B-22) 


The  two  equations  for  e  and  E  are 


e  sin  E  = 


r  r 
/ya" 


(B-23) 


e  cos  E  = 


1-i 


(B-24) 


Dividing 


tan  E  =  *-■?-  & 

/y  (a  -  r) 


(B-25) 
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from  which  E  can  be  determined.  Note  that  if  r  is  positive,  the  radial 
velocity  is  away  from  the  focus  and  E  is  less  than  180°.  Knowing  E, 
the  eccentricity  may  be  determined  from  either  Eq  (B-23)  or  (B-24) . 

The  expression  for  cos  6  for  the  parabolic  case  is  easily 
derived.  Since  the  total  velocity  can  be  determined,  i.e., 

S2  =  —  +  r2  +  r  if  (B-26) 

r 

•  • 
and  the  radial  velocity  r  is  known,  the  transversal  velocity  r0 

must  be 


(r0)2  =  S2  -  r2  =  —  +  r  if  (B-27) 

r 

The  angular  momentum  h  is  then 


h  -  r  (r6)  r +  r  if  (B-28) 


Since  h2  =  y A 


h2  =  rp  +  r3  if 


(B-29) 


Z  = 


The  general  equation  of  the  parabola 


1  +  cos  6 


(B-30) 


(B-31) 
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Appendix  C 


The  Error  Involved  in  Functions 
of  Several  Variables 

Suppose  that  U  =  f(x,y)  is  a  continuously  differentiable  function 
on  some  set  of  values  of  the  arguments  x  and  y.  Suppose  x  and  y  are 
replaced  by  their  approximate  values  a  and  b.  Then  the  approximate 
value  of  the  function  is 


u  =  f($,b) 


(C-l) 


Let  Aa  and  Ab  be  the  errors  in  the  arguments,  i.e., 


./ 


A  =  a  +  Aa  ,  B  =  b  +  Ab 


(C-2) 


where  A  and  B  are  the  exact  values  of  the  arguments.  The  exact  value 
of  the  function  is  then 


IJ  =  f(a  +  Aa  ,  b  +  Ab) 


(C-3) 


Using  Taylor's  theorem  for  functions  of  two  variables 


U  =  f  (a,b)  +  A. 


ML 


a 
y=b 


+  Ab 


3f 


3y 


Jx=a 

y=b 


2! 


A  3  ,  .  3 

~  +  Ab  97 


3x 


f(x>y)x=a  +  *  * 

y=b 


(C-4) 
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Assuming  the  errors  are  small  and  dropping  the  terms  containing  squares 
of  errors 


Au  =  U  -  f(a,b) 


3f 

A„  + 

3f 

3x 

a 

x=a 

_sy. 

*b 


y=b 


y=b 


(C-5) 


where  Au  represents  an  approximate  value  of  the  error  (since  only  the 
first  two  terms  of  the  expansion  were  retained) .  Then 


y=h  y=b 


Now  let  Eu,  Ea,  and  E^  be  the  maximum  absolute  values  of  error. 
Then 


u 


3f 

E_  + 

3f 

3x 

a 

x=a 

3y 

y=b  y=b 


(C— 7) 


Given  the  maximum  errors  of  the  arguments  of  a  function,  the  maximum 
error  in  the  value  of  the  function  can  be  obtained.  If  Ea  and  E^  are 
unknown,  the  conditions  necessary  to  minimize  the  total  error  can 
still  be  obtained  by  examining  the  partial  derivatives  3f/3x  and 
3f/3y.  The  derivation  can  be  extended  to  functions  of  more  than  two 
variables  with  analogous  results.  This  derivation  and  further  support¬ 
ing  material  may  be  found  in  Ref  20:34. 
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Appendix  D 

Derivation  of  Equations  for  True 
Anomaly  Determination 

Figure  46,  an  illustration  of  the  vehicle's  path  across  the 
celestial  sphere,  will  be  used  in  the  derivation  to  follow. 


GA/AE/69-1 


From  the  spherical  triangles  bed  and  acd 


cos  <J<2  =  cos  4*0  cos  (A0g)* 


(D-l) 


cos  =  cos  ij>0  cos  (A0  +  A0g)* 


,(D-2) 


where  A0g  =  0S  -  02  and  A0  =  02  -  0j. 


Dividing  the  above  equations  yields 


cos 

cos  ij>2 


cos  (A0S  +  A0) 
cos  A0„ 


which  reduces  to 


tan  A0g  =  cot  A0  -  cos  4^  sec  4<2  esc  A0 


(D-3) 


(D-4) 


The  angles  4^  and  4>2  can  be  obtained  by  direct  measurement  or  by  the 


vector  dot  products 


I  •  r,  =  cos 


X  •  r2  =  cos  4>2 


In  polar  coordinates,  the  dot  products  become 


4'j  =  cos  61  cos  6S  cos  (Oj  -  ag)  +  sin  6^^  sin  5£ 


cos  =  cos  <$2  cos  6S  cos  (a,  -  as)  +  sin  6,  sin  6 


2  s' 


'2  ws 


*  tp0  is  termed  "minimum  coaltitude"  in  Bef  17 : 46. 


(D-5) 


(D-6) 
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By  direct  measurement,  the  angles  iji  and  <J/2  can  be  determined  from 

$  =  180°  -  (ys  +  A/2)  (D— 7) 


as  can  be  seen  from  Fig.  26. 

The  angle  A6  represents  the  true  anomaly  difference  between 

A  A 

rj  and  r2-  The  vector  dot  product  yields 


cos  A  6 


rxl  rx2  +  ryl  ry2  +  rzl  rz2 


(D-8) 


or 


cos  A0  =  cos  61  cos  62  cos  (a2  -  +  sin  61  sin  <$2  (D-9) 

All  of  the  quantities  appearing  in  Eq  (D-4)  can  therefore  be  evaluated 
and  A0g  can  be  determined.  With  A0S  known  for  one  position,  the 
angle  ipD  can  be  calculated  from 

cos  il>9 

cos  i p  =  - T-r—  0  <  1 1>0  <  90°  (D-10) 

°  cos  A0O  -  o  — 


and  this  angle  remains  constant.  Then  for  any  subsequent  measurement 
using  the  same  star 


cos 

cos  A6S  =  - r— 

cos  4»0 


-  cos  (ys  +  A/2) 
cos  <|>0 


(D-ll) 
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Appendix  E 


Equations  for  Calculating  £,  and  k 


This  Appendix  gives  the  equations  for  calculating  E,  F,  and  k  for 
use  in  the  Hamer  method  of  range  determination. 
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The  angles  E  and  F  (see  Fig.  47),  as  well  as  the  value  k,  pertain¬ 
ing  to  a  reference  trajectory  can  be  pre-calculated  for  any  given  star 
with  known  direction  cosines  £,  m,  and  n. 


To  Star 


Fig.  47 

Illustration  Showing  Vectors 
b  and  h  (From  Ref  8:23) 

The  angle  F,  which  is  the  angle  measured  from  the  vehicle  between 
the  lines-of-sight  to  a  star  and  a  body  center  (earth,  for  example), 
is  obtained  from  the  dot  product  of  the  position  vectors  of  the  earth 
center  and  star  so  that 


£  xev  +  m  yev  +  n  zev 

cos  F  =  - 

rev 


(0°  <  F  <  180°) 


(F,-l) 


By  defining  a  vector  h  perpendicular  to  the  instantaneous  earth-moon- 
vehicle  plane,  the  angle  E,  which  is  the  angle  between  the  line-of- 
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sight  to  a  star  and  its  projection  in  this  plane,  is  obtained  from  the 
dot  product  of  the  vector  h  and  the  unit  position  vector  of  the  star 
giving 


sin  E 


%  hx  +  m  hy  +  n  hz 


(E-2) 


where  the  absolute  value  can  be  used  such  that  0°  <  F.  <  90°  and  where 

\ 

the  values  for  h  ,  h  ,  and  h  are  obtained  from  the  cross  product  \af 
x  y  z  V 

the  position  vectors  of  the  earth  and  moon  centers.  Thus,  ' V 


h  =  v  z  -  z  y 
x  Jev  mv  ev  Jmv 


z  x  -  x  z 
ev  mv  ev  mv 


x  y  y v\x 
ev  mv  ev\  mv 


* (E-3) 


h  =  (hx2  +  hy2  +  hz2)*  2 

\  ,  \  -w  ; 


A  VvV 


(E— 4) 


The  sign  of  the  quantity  k  in 

\ 

position  of  the  star  with  respect 


’*  V  \  Vu:’  \  . 

l  Eq  V8-JL33  i^.  aeteimined  by  the 
to  thV  earth-venicle  line,  r. 


The  sign  is  positive  (negative)  if,  as  viewed  frorn  the  direction  of 

h,  the  projection  of  the  l^ie  tq'  the  star  in  the  instantaneous 

^  ,  ' 

eprth-moon-vehicle  plane  is  i^o'thfe  right  (left)  of  the  earth-vehicle 


line.  As  shown  in  Fig.'  47,  the  relative  directions  of  these  two 
lines  at  4 given  time  along  a  reference  trajectory  can  be  determined 


EiVst  d  tit  in 
\  ^  1 ; 


by  fi*st  dA?i4inft  a  vector  b  which  is  perpendicular  to  the  plane 
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containing  the  vectors  rgv  and  h. 

The  vector  d  is  obtained  from  the  cross  product  of  the  vectors 
rgv  and  h  such  that  the  components  of  b  are  given  by 


b  =  y  h  -z  h 
x  '  ev  z  ev  y 


b  =  z  h  -  x  h 
y  ev  x  ev  z 


b  =  x  h  -  y  h 
z  ev  y  ev  x 


(E-5) 


where  the  values  of  h  ,  h  ,  and  h  are  given  in  Eq  (E-3) .  Next,  the 

x  y  z 

dot  product  of  the  unit  position  vector  of  the  star  and  the  vector  b 

will  yield  the  component  of  the  star  vector  along  the  vector  b  which  is 

given  by  the.  term  (£  b  +mb  +nb  )/b.  If  this  term  is  positive, 

x  y  z 

the  projection  of  the  unit  star  vector  in  the  instantaneous  earth-moon- 
vehicle  plane  is  to  the  right  of  the  earth-vehicle  line  so  that  the 
quantity  k  is  given  as 


£  b  +mb  +nb 

k  -  - * - n  <*-« 

U  bx  +  m  by  +  n  bz| 


The  quantity  k  is,  of  course,  always  either  plus  or  minus  one. 
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